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Denkt man sich die tragende Fläche aus einzelnen Hufeisenwirbeln aufgebaut, so läßt sich in an- 
schaulicher Weise eine einfache Formel für das Geschwindigkeitspotential finden, die den bekannten Formeln 
von J.M. Burgers, Th.v. KärmänundL. Prandtl gleichwertig ist. 


The author assumes that the supporting plane consists of single vortices of horse-shoe shape, and deduces 
in an obvious way a simple expression for the velocity potential that is equivalent to the known formulas of 
J.M. Burgers, Th.v. Karmanand L. Prandtl. 


Suppose que le plan portant consiste de tourbillons singles de forme de fer & cheval, d’une maniere 
claire une formule simple pour le potentiel de vitesse peut ötre trouvee, laquelle Egale les formules connues 
de J.M. Burgers, Th.v. Kärmänet L. Prandt!. 


Ilyerp Hecyımas HOBEPXHOCTB COCTOHT U3 OTACAIBHBIX IIONKOBOOOPASHEIX BUXpeü; TorAa Aa 
HOTEHUHAJIA CKOPOCTeH MOST ÖbITB HauNeHa, BecbMA HATIAAHBIH CIIOCOÖOM, IIPocraa dop- 
My.la, PABHOCHIBHAA H3BECTHBIM PopMyıaMm Byprepa (J. M. Burgers), Kapmana (Th. v. Kärmän) 
u HOpanıras (L. Prandt]). 


I. Einleitung 


Die Behandlung von Tragflügelfragen, insbesondere die Berechnung des induzierten Ge- 
schwindigkeitsfeldes, gestaltet sich besonders einfach, wenn man von der tragenden Fläche das 
Geschwindigkeitspotential kennt. 

Erstmalig wurde eine Formel für das Geschwindigkeitspotential von J.M. Burgers [1] 
angegeben. Diese Formel, deren Ableitung im einzelnen etwas langwierig ist, wird aus der all- 
gemeinen Wirbeltheorie gewonnen. Th.v.Kärmaän [2]leitet ebenfalls auf der Grundlage der 
Wirbeltheorie eine geschlossene Formel für das Geschwindigkeitspotential ab, indem er die 
Methode der Fourierschen Integrale anwendet. Diese Formel ist noch recht kompliziert. 
Sie läßt sich aber durch Ausführung gewisser Integrationen !) in die Formel von Burgers 
überführen. Einen anderen Weg ist L. Prandtl [3] gegangen. Aus der Vorstellung der mit 
Druckdipolen belegten Fläche findet er eine Formel für das Beschleunigungspotential. Hieraus 
läßt sich in einfacher Weise das Geschwindigkeitspotential berechnen ?). 

Alle bisherigen Ableitungen sind vom physikalischen und mathematischen Standpunkt aus 
z. T. nicht immer ganz leicht zu verfolgen. Dieses veranlaßte uns, nach einer einfacheren und 
anschaulicheren Lösung zu suchen, die wir dadurch finden, daß wir uns den Flügel aus einzelnen 
Hufeisenwirbeln aufgebaut denken. Wir behandeln im folgenden den stationären, inkompressib- 
len Fall; glauben aber, daß sich die hier gegebene Darstellung auch auf den Fall kompressibler 
Strömung (Unter- und Überschallströmung) erweitern läßt. Daß die Strömung außerdem noch 
homogen, reibungslos und schwerefrei sein soll, sei der Vollständigkeit halber noch erwähnt. 


II. Das Wirbelgebilde der tragenden Fläche 


In Anlehnung an eine bereits bei H. Glauert [4] angegebene Darstellung kann man sich 
den Tragflügel aus sogenannten Elementarflügeln zusammengesetzt denken. Unter 
einem Elementarflügel verstehen wir einen Flügel von der Breite des Spannweitenelementes dy’. 
Jedem Elementarflügel ist im Gesamtverband des Flügels ein ganz bestimmter durch die Geo- 
metrie des Flügelgrundrisses gegebener Platz zugeordnet. Für eine Tragfläche beliebiger 
Grundrißform gibt Bild 1 eine schematische Darstellung. Kennt man die Wirkung aller Ele- 
mentarflügel, so findet man die Gesamtwirkung durch Integration über die Spannweite (y’). 

Wie wir bereits in der Einleitung sagten, wollen wir den Flügel aus Wirbelsystemen auf- 
bauen. Um uns ein klares Bild über die zu einem Elementarflügel gehörenden Wirbel zu machen, 
gehen wir noch einmal kurz auf die Entstehung der Zirkulation an einem auftrieberzeugenden 
Flügel ein. Beim Ingangsetzen des Tragflügels führt das Umströmen der Hinterkante mit 


1) Siehe Anhang, 
1 
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theoretisch unendlich großer Geschwindigkeit zur Bildung einer Trennungsfläche, die sich zu 
einem Wirbel, dem sogenannten Anfahrwirbel, aufrollt. Da dieser Wirbel immer an dieselben 
Flüssigkeitsteilchen gebunden ist, schwimmt er mit der Strömung fort. Auf Grund der Helm- 
holtzschen Wirbelsätze muß am Flügel eine Zirkulation entstehen, die entgegengesetzt gleich 
groß ist wie die des Anfahrwirbels. Diese Zirkulation 

| Anströmrichtung wirdnunnach Kutta-Joukowski solange wachsen, 

bis an der Hinterkante glatter Abfluß herrscht. Die Ver- 

bindung zwischen dem am Flügel verbleibenden Wir- 
bel, der auch als gebundener Wirbel bezeichnet wird, 
und dem Anfahrwirbel bilden die sogenannten freien 
Wirbel. Diese Wirbelfäden folgen den Stromlinien. 
Entsprechend den hinter dem‘ Flügel herrschenden 
senkrechten und seitlichen Geschwindigkeitskompo- 
nenten beginnt sich das gesamte Gebilde der freien 
Wirbel von der Flügelhinterkante her aufzurollen. Im 
Rahmen einer linearisierten Theorie, d.h. bei Beschrän- 
kung auf kleine Profilwölbungen der einzelnen Flügel- 
schnitte und kleine Anstellwinkel, wollen wir den Auf- 
rollvorgang vernachlässigen und annehmen, daß die 
freien Wirbellinien geradlinig in der Ebene der Trag- 
fläche bis ins Unendliche reichen. Genau so wie schon 
in der zweidimensionalen Skelett-Theorie dünner 
schwachgewölbter Profile, vgl. H. Glauert [4] und 


u b ; ; 
Bild 1. Grundrißplan der tragenden Fläche. W.Birnbau “ [5], nehmen eat daß sich die 

a valte gesialeltigeie gebundenen Wirbel in der Flügelebene (Sehnenebene) 

dy,= Spannweitenelement des Elementarflügels 1 

(y’ = Tiefe des Elementarflügels befinden. 


Die um das Flügelprofil herrschende Zirkulation 
verteilt sich über die ganze Flügeltiefe. Bezeichnet k (#’, y’) die Wirbeldichte in einem Punkt 
x, y, 2’ =(, dann schreiben wir für die Wirbelstärke in diesem Punkt 


al (a; Yy) = klar, y)dr era ee 


im Fall des Elementarflügels hat dieses Wirbelelement die Länge des Spannweitenelementes dy’. 
Wie wir nun oben erläutert haben, gehen von diesem Wirbelelement zwei bis ins Unendliche nach 
hinten reichende gerade Wirbellinien konstanter aber entgegengesetzter Zirkulation der Stärke 
dI' ab. Sie haben ihren Ursprung in der Ebene 2° =0 in den Punkten #°, y’ und #, y’+dy. 
Ihr seitlicher Abstand entspricht also dem Spann- 
weitenelement dy’. Wir bezeichnen diese beiden Wir- 
bellinien als freies Wirbelpaar. Das gesamte Ge- 
bilde aus gebundenem Wirbelelement und freiem y 
Wirbelpaar nennen wir Tragflügelelement. ‚ 


(z'=0) 


‚bebundenes Wirbelelemenf 


Freies Wirbelpaar 


Bild 2. Wirbelgebilde eines Elementarflügels. Bild 3. Wirbelgebilde eines Tragflügelelementes 
u,v,w = induzierte Geschwindigkeitskomponenten 


Es ist gleichbedeutend mit dem bekannten Hufeisenwirbel. Bild 2 und 3 geben die schematische 
Darstellung des Elementarflügels und des Tragflügelelementes. 


Unseren Betrachtungen liegt das folgende Koordinatensystem zugrunde: 


x-Achse: Schnittgerade der Flügelebene mit der Strömungsebene; positiv in Anströmungs- 
richtung, 


an 
A TE; 
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y-Achse: Schnittgerade der Flügelebene mit der senkrechten Strömungsebene; positiv 
nach rechts, 


z-Achse: Senkrecht zur Flügelebene, positiv nach oben. 


III. Die induzierten Geschwindigkeitskomponenten 


Unter Anwendung des Biot-Savartschen Gesetzes wollen wir jetzt die Wirkung des 
im vorigen Abschnitt beschriebenen Wirbelsystems berechnen. 

Zunächst ermitteln wir die durch das Tragflügelelement, vgl. Bild 3, induzierten Ge- 
schwindigkeitskomponenten. Die Geschwindigkeiten rechnen wir positiv in Richtung der 
positiven Achsen. Der Abstand eines Punktes x, y, 2, an dem wir die Geschwindigkeitskom- 
ponenten suchen, vom Ort des gebundenen Wirbelelementes #’, y’, 2° = 0 beträgt 


r= Mr — vo) +(y— y) +2 en ie (2). 


Das gebundene Wirbelelement liefert: 


dr HEN) & away. DEREREISEERNEDN RT: (3) 


dr 
a. . (4), 
k(z', Y') dass x ’ ’ 
= BT EZRTE Gr Re ORERTRER) 
Um die Wirkung des freien Wirbelpaares zu bestimmen, gehen wir von der Wirbellinie aus, 


die im Punkte #°, y’, 2° = 0 ihren Ursprung hat und geradlinig nach hinten bis ins Unendliche 
reicht. Hat diese Wirbellinie die Zirkulation d!’= k.d«’, dann gilt 


dv —=0 


dw — 


a ie ie ne (6); 
k(#, y’‘) 2 | e— X h 

dv — a 1 + = ) de Be 00207) 
key) yo | T—a@\ ,, 

Be. me 14) Ge 1 (R): 


Nehmen wir jetzt die zweite im Abstand Ay’ parallel verlaufende Wirbellinie hinzu, die die 
gleiche aber entgegengesetzte Zirkulation wie die erste hat, dann finden wir die Geschwindigkeits- 
änderungen durch das freie Wirbelpaar durch die Grenzwertbildung 


= 1 a) — fly’ + Ay‘) WE af y’ 9 
En re erg ® 
wobei f jeweils die Geschwindigkeitskomponenten du, dv, dw der eben betrachteten ersten Wir- 
bellinie bedeutet. Bei der auszuführenden Differentiation gilt die Wirbelstärke k(a#’, y’) da’ 
als Konstante. 

Die Differentiation der Beziehungen (6) bis (8) nach %’ liefert dann nach Einführen des 
negativen Vorzeichens die folgenden durch das freie Wirbelpaar induzierten Geschwin- 
digkeitskomponenten: 


Ge N une. (10), 
ii rer k(a, y') | 2(y —Y)2 %—H W-Y))% Taf ’ ea 
= 47 er ur” ag era N) 
Re a —ı El, 19 
re yı2r I nee! Tor a error TE 


Die Wirkung des Tragflügelelementes erhalten wir jetzt durch Addition des 
Anteils des gebundenen Wirbelelementes und des Anteils des freien Wirbelpaares. Da das freie 
Wirbelpaar keinen Beitrag zur Geschwindigkeitskomponente « liefert, stellt (3) zugleich die Ge- 
samtwirkung in z-Richtung dar. Entsprechend liefert das gebundene Wirbelelement keinen 
Beitrag zur Geschwindigkeitskomponente v, so daß (11) die Gesamtwirkung in y-Richtung an- 
gibt. Für die Geschwindigkeitskomponente w in z-Richtung gilt dagegen 


ka,y)| y-y®’—z ( u —® 2 S—a\,, 07 
u Zul VREGITE ST NACGIF 
Be 4a yy) Fapl 1 ) Wr 
Die Wirkungeines Elementarflügels finden wir jetzt, indem wir die Gleichun- 


gen (3), (11) und (13) über die Flügeltiefe (2°) integrieren. Das Geschwindigkeitsfeld der gesam- 
13* 
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ten Tragfläche erhalten wir, wenn wir nun noch über die Spannweite (y’) integrieren. 
Auf die Wiedergabe der Integrationsformeln können wir verzichten ?). Außerhalb der tragen- 
den Fläche ist die Wirbeldichte %k («’, y’) identisch null. Sie verschwindet wegen der Kutta- 
Joukowskischen Abflußbedingung ebenfalls an der Flügelhinterkante. 


IV. Die Gleichung für das Gesehwindigkeitspotential 
Von dem Ausdruck 


DD, y2, 205 Y) = u As Al | 5 . de.dy a ı ı an (1A) 


läßt sich leicht zeigen, daß die partiellen Differentiationen nach &, %, z die Geschwindigkeits- 
komponenten 


E12 Re, UBER? 
eu=,.(#B), Se ®), d Kr (daD) au ul 


eines Tragflügelelementes darstellen. Da wir außerdem noch zeigen können, daß der Ansatz die 
Laplacesche Potentialgleichung 
A(d?®) =0 


erfüllt, stellt also (14) das Geschwindigkeitspotential eines Tragflügelelementes dar. 
Integrieren wir jetzt über die Flügelfläche, so erhalten wir das Geschwindigkeitspotential 
der gesamten tragenden Fläche zu ?) 


oo 


x — % 


| engere Er Br) ir 06) 


—o 


Wir haben die Grenzen mit — oo bis oo angegeben; merken aber an, daß außerhalb des Flügels 
die Wirbeldichte k (x, y’) = 0 ist. 
Wir können (16) noch etwas anders schreiben, wenn wir die Abkürzung 


’ te A Fr 4 a re rk De — “ . . . . 
Ga, nun) |M De (17) 


einführen. Diese Funktion, mit 'z/ [(y — y’)® + z?] multipliziert stellt den Einfluß eines Schnittes 
y’ auf einen Punkt x, y, z dar. a für den Fall, daß y = y und z = 0 werden, gilt 


—o 


G@(2, y,0,y)= für nl + : 2) ne =n2Ll2)y) (18), 
Es ist 
Ioay)= Ben KR) Y)dEE one Do ee 


die an einer Stelle z, y von der Vorderkante des Flügelschnittes y aus integrierte Zirkulation. 
Wir setzen (17) in (16) ein und erhalten 


1 2 
D (x, y, 2) = af (ee N array Re ANDO), 


Wir können der Gleichung (16) noch eine andere Form geben, wenn wir bei der Integration 
über &’ zunächst eine partielle Integration vornehmen. Es wird dann 


oo 


1:0 ; I'(z, y'‘) 
©(2, y,2) = ae > Kane 
By) 2 ee dd (21) 


Die Geschwindigkeitskomponenten erhalten wir dureh partielle Differentiationen der Formel 
für das a zu 


uw) my), a2) = By), wa) ann) ' 2). 


#) Ganz ähnliche Formeln zur Berechnung der Geige ir bereits von J.M. 
Burgers [l], S. 139 angegeben worden. 
3) Eine ähnliche Formel gibt auch wieder J.M. Burgers[l], S. 128 an. 
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V. Sonderfälle 


a)Das Geschwindigkeitspotentialbeisenkrechter Annäherung 
andie Flügelebene 
Einige Aufmerksamkeit erfordert die Untersuchung, wenn der Aufpunkt z, 9, z in die 
Nähe der Flügelebene wandert. Mathematisch bedeutet das, daß der Grenzwert lim z—0 zu 
bilden ist. 
Wir entwickeln das Geschwindigkeitspotential in eine Taylor sche Reihe nach z und be- 
trachten nur Glieder 1. Ordnung. 


Da 1,2 0)= Da, 0) (my 0)e +: ee (23). 


Um die auftretenden Grenzwerte bilden zu können, müssen wir die Voraussetzung machen, daß 
die Funktion @(z, %, 2, y') in der Umgebung des Punktes y’ = y,2 = Ostetig nach y’ differenzier- 
bar ist. 

Differenzieren wir (17) partiell nach y’ und setzen anschließend y = y und z=0, dann 
folgt: 


0G@ k 2 —ıa ol 
a9 (2, 9,0, y) = J > (wi, nl a2 a) da’ = 2 dy (U) +. (24). 
Um also die stetige Differenzierbarkeit der Funktion @(&, %, z, y’) in der Umgebung des Punktes 
y = y,2=0 zu erreichen, muß die Zirkulationsverteilung /'(&, y’) in Spannweitenrichtung (%‘) 
im Punkt y’ = y stetig nach y’ differenzierbar sein. 

Wir wollen darauf verzichten, im einzelnen auf die mathematisch nicht uninteressanten 
Grenzwertbildungen einzugehen, sondern geben sofort das Ergebnis der Gleichung (23) an. 


Ö(z, y,2—0) = 


| 


m T(z, y) .- Im E G(&, y,0, y) - Bi = en 2 ir) (25). 


Das Zeichen $ bedeutet, daß der Bereich y—e< y’ <y-+e von der Integration auszulassen 
ist. Das Integral ist also im Sinne des Gauch yschen Hauptwertes zu nehmen. 

Die Formel (25) ergibt die bekannte Tatsache, daß sich das Potential am Ort der Flügel- 
fläche oder in der hinter dem Flügel abgehenden Wirbelschicht von der Oberseite (2>0) zur 
Unterseite (2< 0) sprunghaft um den Betrag der Zirkulation /'(z, %) ändert. 

Die induzierten Geschwindigkeitskomponenten erhalten wir aus (25) zu 


iE2207; 1e% 
DEN are Dir )E 
#2’ al 
v(®, Oz, Y)» Eng une (20), 
2 Z h = 
ea m en (@ er En un 


Durch partielle Integration können wir die letzte Formel HER umschreiben und erhalten 


| Te, y | 
w(z, y,0) = an. (x, y) = Kr u az da’ ir (27). 


Für Punkte, die sich inunendlicher Entfernung hinter dem Flügel 
befinden, geht die Funktion @ nach (17), vgl. auch den nächsten Abschnitt, über in 
Bere Our + (28), 
wobei /'(y') = I'(», y’‘) die Gesamtzirkulation um den Schnitt y’ bedeutet. 
Setzen wir diesen Ausdruck in (25) ein, dann erhalten wir 


set Dur) ea. :..20 9). 


Die neue Größe 


ed = lim ya aa yı- das ln ei. (80) 
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stellt nun gerade den nach Prandt]1 bekannten induzierten Abwind dar. Meistens findet man 
diese Formel so geschrieben, daß an Stelle der Zirkulation /’die Ableitung d]/dy darin vorkommt. 
Man erhält diese Form durch partielle Integration von (30). 


1 Bra 


Auf die Wiedergabe der Formeln für die er wird an dieser Stelle 
verzichtet. Sie werden im folgenden Abschnitt, Gl. (40), mitgeteilt. 


b)Die Verhältnisse weit vor und hinter dem Flügel 


Wir wollen das Geschwindigkeitspotential in größerer Entfernung vor bzw. hinter dem 
Flügel bestimmen. Hierzu nehmen wir an, daß 


vr 7r>y—y)’+2. 
Die Funktion @(z, y, z, y’) nach (17) läßt sich dann folgendermaßen entwickeln: 


@(2] >o, y, 2, y) = en 7 ah 7) (y — y’) +21! 7). 22 
wobei 


x 


I) = 10.9), | Kiasu)de 


—» 


die Gesamtzirkulation um den Schnitt ist. 
Wir setzen jetzt (32) in (20) ein und erhalten 


Can f , 2 i % f - 
Dial > 0, y,2) = Wer ES) al ee jrwar. .(84). 


Für © = ® stellt das erste Glied das Geschwindigkeitspotential in unendlicher Entfernung 
hinter dem Flügel dar ®). 


en | 
5(®, WER eat. ET PETER 


Das zweite Glied bildet die Integration der Zirkulationsverteilung über die Spannweite. Nach 
Kutta-Joukowski gilt für den Gesamtauftrieb 


A=oV[Ty)ay . 22 ece nn). 


Es bedeutet V die Anströmungsgeschwindigkeit und og die Dichte des strömenden Mediums. 
Wir können für (34) also auch schreiben 


2.5 4 r 
ld, zit 2) oe. Vo en 
Die Geschwindigkeitskomponenten erhalten wir zu 
ul >, 2) = en 
: 4 el? oV’ 
lin = (142) ve 12), ER: 
1 x ae ne | 
w(2] —@, v2) =5(1+5,)w(, y) — gen ar, 
mit 
ae MIETE 
0(0, 4,2) = = I 


uz ber = Ä 
w(&, Yy,2) = Wy) Fap dy’ 
‚*) Die Formel (35) kann man auch ändunh herleiten, daß man sich die abgehende Wirbelschicht aus 
einer in y-Richtung verteilten Dipolbelegung hergestellt denkt. 
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Die Werte in der Wirbelschicht lassen sich aus den Ergebnissen des vorigen Abschnittes leicht 
bestimmen. 


u(&, y,0) =0 | 

dr 
v(o, ON N we (AU). 
w(@, 4,0) = — 2 w;(y) | 


Während für die w-Komponente nach (38) zu dem Wert im Unendlichen noch ein von x ab- 
hängiges Korrekturglied hinzukommt, fehlt ein ähnliches Glied bei der v-Komponente. 


ohDer Flügelunendlichen Seitenverhältnisses (Ebenes Problem) 
Für den Fall des ungepfeilten Flügels unendlichen Seitenverhältnisses erhält man bei 
konstanter Flügeltiefe auf Linien gleicher Rücklage in Spannweitenrichtung (%) eine konstante 
Wirbeldichte. 
KEN. 20 ar suay or rl). 
Setzen wir diesen Ausdruck in (16) ein und integrieren zunächst über y’, dann finden wir die be- 
kannte Beziehung 


Bi K “Var — — ß BER En > R 
z | k(x’) arc eig m ) a (42) 


In der Flügelebene (2 = O)ist vor dem Flügel das Potential null, während am Ort der Fläche und 
hinter dem Flügel der Potentialsprung 


ra, SV er ey 


besteht. 


Die Geschwindigkeitskomponenten erhalten wir durch partielle Differentiationen von (42) 
nach & und z zu 


1 i 2 , 
u(X,2) = 37 | k(x) Gern da’, 


oo . (44). 
1 ee —% h 
Bei senkrechter Annäherung an die Flügelebene wird daraus 
1 2 1 D k(«’) ‚ r 


Die letzte Beziehung können wir durch partielle Integration, ähnlich wie im allgemeinen Fall (27), 
wieder umformen. 
Ba) | 


Tue 
| ee Ar Een tn. (46), 


w(z,0) = — 


d) Der Flügel verschwindenden Seitenverhältnisses 
Der andere Grenzfall ist nun derjenige, bei dem die Spannweite im Vergleich zur Flügeltiefe 
sehr klein ist. Auf Grund der gemachten Vernachlässigungen bezüglich des Aufrollvorganges 
hinter dem Flügel gelten die nachfolgenden Betrachtungen nur für schr kleine Anstellwinkel. 
Wir wollen uns auf Punkte in unmittelbarer Nähe der Flügelebene (z— 0) beschränken. 
Für den vorliegenden Fall können wir schreiben 


(.— a >(y—yY). 


Setzen wir diese Beziehung in (17) ein, dann erhalten wir 


@(z, y,0, y) = | re. Yy) (1 AR r =, )de = 21 (m; y‘) dauadelanld7) 
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und nach (25) 
2 


1 
D(2,y,2>0)= 3% Ney) 220, (2%) 0. a (A 
Hierin bedeutet gemäß (30) 
Pie CT(a,y) , 
N ler N 


den induzierten Abwind für die Zirkulationsverteilung /'‘(z, ). 
Von den Geschwindigkeitskomponenten wollen wir nur die w-Komponente angeben. 


Sie lautet 
wen y) 2 ul) DLR Er 0 


e)Dietragende Linie 
Denkt man sich die über die Flügeltiefe verteilte Zirkulation in einem Punkt zusammen- 
gezogen, so gelangt man zum Gedankenmodell der tragenden Linie. Längs dieser Linie herrscht 
dann an jeder Stelle y’ die diesem Schnitt entsprechende Gesamtzirkulation J'(y’). Bezeichnen 
wir mit % = x(y’) die Lage der tragenden Linie in x-Richtung, dann können wir nach (17) für 
die Funktion @(z, y, 2, y') den folgenden Ausdruck herleiten: 


—% 
G(&, y, 2, y’ en TEEN 
RL. rn Nez) + ie) n 
Es ist interessant festzustellen, daß man diese Formel auch dadurch findet, daß man (17) für 
große Abstände r des Aufpunktes z, %, z vom Ort der tragenden Linie 7, y’ entwickelt. Hieraus 
ergibt sich dann auch, daß 


ge 
a 2 
[ k(&, y') di’ 


gerade der Schwerpunktsabstand der über die Tiefe verteilten Zirkulation ist. 
Für den Wert der Selbstinduktion (y’ = Y, 2 =0) folgt aus (51): 


Kay) ETW...) 


Das Geschwindigkeitspotential eines Flügels, der durch eine tragende Linie dargestellt wird, 
erhalten wir, wenn wir (51) in (20) einsetzen. 


1 | 2 x — u 

29,2) = 1 DNa)darııa Hr 
ne RER Fl eg ) wi 
Die Geschwindigkeitskomponenten erhält man wieder durch partielle Differentiationen nach 
&, Y, 2. Für alle anderen Formeln der vorigen Abschnitte, in denen die Funktion @ vorkommt, 

kann diese für den vorliegenden Fall der tragenden Linie durch (51) und (53) ersetzt werden. 
Die auf Seite 169 gemachte Voraussetzung der stetigen Differenzierbarkeit der Funktion @ 
nach y’in der Umgebung des Punktes y’ = y, z = ist erfüllt, sofern + zist. Es ist nämlich 


@ t—E\af 
rd a: 


Die Zirkulationsverteilung I(y’)in Spannweitenrichtung muß also stetig nach %’ differenzier- 
bar sein. 

Für.den ungepfeilten Flügel, bei dem wir die Lage des Koordinatenanfangspunktes i in &- 
Richtung in den Ort der tragenden Linie verlegen, ist 7 =. 

Die Formel (54) läßt sich sehr einfach anwenden, wenn die tragende Linie entweder schräg 
steht (schiebender Flügel) oder gepfeilt ist. Man gewinnt die entsprechenden Formeln durch 
einfache Koordinatentransformationen. 


| — 


VI. Zusammenfassung 
Es wird eine einfache und anschauliche Ableitung für das Geschwindigkeitspotential der 
tragenden Fläche bei inkompressibler Strömung angegeben. Die Formel ist den bisher bekannten 
Formeln von Burgers,v.Kärmän und Prandt| gleichwertig. Einige Sonderfälle, wie 
z. B. der bei senkrechter Annäherung an die tragende Fläche, werden besonders untersucht. 


R 
u: 


D- 
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VIH.:Anhang 
Im folgenden geben wir die bereits bekannten Formeln für das Geschwindigkeitspotential 
an. Die Formeln, die wir in unseren Bezeichnungen anschreiben, formen wir jeweils noch um, 
um sie mit der unsrigen vergleichen zu können. 


a) Formel.von J.M. Burgers [1], S. 128 
Aus der allgemeinen Wirbeltheorie folgt: 


2’ 


1 7 2 — 
KB * ’ E ’ ’ 
o- ||] * dz er dy'. 


Da wir annehmen, daß sich sämtliche Wirbel in der Flügelebene 2’ = 0 befinden, braucht eine 
Integration über 2’ nicht ausgeführt zu werden. ‚Wir können.also setzen k*de’ = k(«, y‘). 
Setzen wir die Beziehung nach (2) für r ein, dann können wir die Formel leicht in unsere For- 
mel (16) überführen. 


b) Föormelvon Th.v. Kärmän [2] 
Unter Zuhilfenahme der Methode der Fourierschen Integrale wird: 


0-0, +6, 
mit 
= | | | [Her mE une" aa ayı dad 
EI TS PR 4 
und 


1 f s he 
en | | | k(@, y)cosu(y—y)e "da’dy du. 
—o —o 0 
Führt man die Integrationen über A und « in den Grenzen von 0 bis oo aus, so gelangt man un- 
mittelbar zu unserer Formel (16) 5). 


c) FöormelvonL. Prandt] 3] 
Für das Beschleunigungspotential gilt 


oo 


= nn | | (Po— Pu) j da’ dy'. 
Hierin bedeutet % 
le kr 1407 —u,) 5 Vhce, y‘) 
den Potentialunterschied am Ort der Flügelfläche. 
Zwischen dem Beschleunigungspotential und dem Geschwindigkeitspotential besteht noch 
der Zusammenhang 


x 
1 
Setzen wir die beiden ersten Beziehungen in diese Formel ein und führen noch die Inte- 


T 
gration | = aus, dann folgt sofort daraus: (16). 
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Der symmetrische Kreisel mit zeitlich periodischem 


Richtmoment 
Von W. Braunbek in Tübingen 


Die Bewegung des symmetrischen Kreisels unter einem konstanten plus einem zeitlich harmonisch ver- 
änderlichen Richtmoment auf die Figurenachse wird untersucht. Denkt man sich das Richtmoment durch 
Einwirkung eines homogenen Magnetfeldes auf einen in der Figurenachse eingebauten Magnetstab erzeugt, 
so sind zwei Fälle zu unterscheiden: Wechselfeld parallelundsenkrechtzum konstanten Feld. Im 
ersten Fall ergibt sich Resonanz, wenn die Frequenz des Wechselfeldes mit der Frequenz der freien N utation 
des Kreisels, im zweiten, wenn sie mit der freien Präzession übereinstimmt. Außerdem treten auch subhar- 
monische Resonanzen auf. 


The motion of a symmetric gyroscope is investigated, if @ constant and a harmonically variable directing 
moment are acting atihe awis. Ifthe directing moment is generated in such a way that a homogenous magnetic 
field acts on a bar magnet that is situated in the awis, two cases are to be distinguished, namely that the alter- 
native field is situated parallel and normal to the constant field. There is resonance in the first case, 
if the frequency of the alternative field is equal to that of the free nutation, in the second case, if that frequency 
is equal to that of the free precession of the gyroscope. Besides, subharmonic resonances are possible. 


Le mouvement de la toupie symmetrique & un moment directeur constant et un tel harmoniquement variable 
al’egard du temps, agissant sur l’aze est examine. Suppos£ que le moment directeur est produit par l’influence 
d’un champs magnetique homogene sur une barre aimentee montee dans l’aze, il faut distinguer deux cas: un 
champs alternatif parallele et verticalau champs constant. Au premier cas il y a une resonance, si 
la frequence du champs alternatif correspond & la frequence de la nutation libre de la toupie, au second cas, si 
elle correspond a la precession libre. De plus, des resonances subharmoniques peuvent se produire. 


Nccuenyerca ABUKeHUe CHMMETPHYeCcKOTO TUPOCKONA, OHTYPHAA OCh KOTOPOTO HOABEPTAETCH 
BJIHAIHHIO IIOCTOAHHOTO INIIOC TAPMOHHYEeCKH BO BPEMEHH IIEPeMEHHOTO HAIIPABJIAIOLILHX MOMEHTOB,. 
Ecım HampaBıaamımHÜü MOMEHT CO3NAeTCH ONHOPOAHBIM MATHHTHEIM HOJIEM, BIIHSIOIIUM Ha MAT- 
HUTHBIH CTeP>KeHb B OCH OUTYPEI, TO HA3INHYAMT ABA CIYy4YaA: TePeMEHHOE NONE TA PAMIENBHO 
U TepHeHAHUKYIAPHOK IIOCTOAHHOMY 1010. Pe30HaHc NOAYyYaeTcH IIPH COBIAACHHH YAYTOTBL 
TIePeMEHHOTO NO.IA, B IIEPBOM CHYYAe, — C YACTOTOU CBOÖOAHOH HYTAIHH THPOCKOIA, a BO BTOPOM, 
— c cBObonHoü ımpemeccueii. KpoMe TOTO BOSHHKAWT TAKIKE CYOTAPMOHHYeCKHE Pe3OHAHCHI. 


Einleitung 


Kreisel mit einem richtenden Drehmoment auf die Figurenachse, das von deren Lage und 
außerdem explizit perodisch von der Zeit abhängt, sind auf zwei Gebieten der Physik schon ein- 
gehend behandelt worden: Einmal im Rahmen der Störungen der Erdrotation durch Kräfte 
aus dem Planetensystem (Präzession und Nutation der Erdachse), zum anderen bei der Unter- 
suchung des Verhaltens von Elektronen oder Atomkernen, die infolge ihres mechanischen und 
magnetischen Momentes Kreisel mit Richtkräften darstellen, im magnetischen Hochfrequenz- 
feld. Auf beiden Gebieten sind aber die Untersuchungen stets speziell auf das vorliegende be- 
sondere Problem ausgerichtet gewesen. Die atomaren Kreisel bilden zudem ein quanten- 
mechanisches System, das nur in bezug auf seine Erwartungswerte klassisch behandelt werden 
kann, und weichen auch sonst mehr oder weniger von dem einfachen, makroskopischen Krei- 
sel ab. 


Es soll daher hier einmal in einer mehr systematischen Weise die Bewegung des klassi- 
schen, makroskopischen, symmetrischen Kreisels unter zeitlich periodischen Richtkräften be- 
handelt werden. Das besondere Interesse soll dabei der Resonanz zwischen der aufgeprägten 
Frequenz und charakteristischen freien Frequenzen des .Kreisels, wie etwa seiner Präzessions- 
frequenz oder seiner Nutationsfrequenz, gelten. 


Um die Verhältnisse nicht zu sehr zu komplizieren wirdvonjeder Dämpfungab- 
gesehen. Eine ‚„Resonanzkatastrophe‘‘ kann — im Gegensatz zur linearen Schwingung — 
beim Kreisel auch für verschwindende Dämpfung nicht eintreten; sie wird durch die Nicht- 
linearität der Gleichungen verhindert, oder anschaulich einfach dadurch, daß alle Ampli- 
tuden als Winkelamplituden auf endliche Werte beschränkt sind. 


Den symmetrischen Kreisel mit Trägheitsmoment J um die Figurenachse und Trägheits- 
moment J’ um eine Schwerpunktachse senkrecht dazu denken wir uns im Schwerpunkt unter- 
stützt (Bild 1). Das Richtmoment läßt sich am einfachsten gedanklich so konstruieren, daß 
man in den Kreisel starr einen Magnetstab mit magnetischem Moment u in Richtung der 
Figurenachse einbaut und den Kreisel in ein Magnetfeld 9 bringt. Dieses soll sich zusammen- 
setzen aus einem zeitlichkonstanten Anteil $, und einem Anteil 9, derharmonisch 
mit der Kreisfrequenz w von der Zeit abhängt. Beide Anteile seien räumlich homogen. Stellt a 
einen Einheitsvektor in Richtung der Figurenachse dar, so wirkt auf die Figurenachse das 
Drehmoment: 


M = la, HF 91 2 ER 


® 
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In dem konstanten Feld 9, allein würde sich der Kreisel wie ein symmetrischer, schwerer 
Kreisel verhalten. Für das Wechselfeld 9, sind der Richtung nach zwei typische Fälle zu unter- 
scheiden: 


Fall I: 9, parallel zu 9 
Fall II: 9, senkrecht zu 9. 


Beide Fälle geben zu verschiedenen Bewegungstypen An- 
laß und sollen im folgenden getrennt behandelt werden. Im 
Fall II kann 9, eine elliptische Polarisation in einer Ebene L 9, 
besitzen. Hier interessieren besonders die Spezialfälle der 
linearen und der zirkularen Polarisation. In Fall I gibt es 
natürlich keine Polarisation. Hier schwingt der Vektor 9, 
stets in der festen Richtung 9.- 


Die Lage des Kreisels soll durch seine Euler schen 
Winkel 9, y, ®in einem ortsfesten Koordinatensystem %, %, z 
charakterisiert werden (® Winkel (z, a); y Azimuth der 
z-a-Ebene um die z-Achse; @ Azimuth des Kreiselkörpers 
um die Figurenachse). $, hat stets die Richtung der z-Achse. 
Außer der aufgeprägten Kreisfrequenz w des Feldes 9, spielen 
noch die Kreisfrequenz ®, der Rotation des Kreisels um die 
Figurenachse, die Kreisfregenz w, der freien, regulären Präzes- 
sion im Feld 9, und die Kreisfrequenz w; der freien Nutation / 
in diesem Feld eine Rolle. Bild 1. Magnet-Kreisel im Magnetfeld 


Fall I. Wechselfeld parallel Grundfeld 


'$, und 9, haben die Richtung der 2-Achse. Also gilt: 
H,=H,+ H,coswt=H,(1+.&cos wit), 


wenn a als Abkürzung für das Verhältnis 4,/H, eingeführt wird. Auf die Figurenachse wirkt 
um eine Achse senkrecht zur z-a-Ebene das Drehmoment (positiv, wenn es # zu vergrößern 
trachtet): 

= — uH,(l-+0cos wi) 


und für die Euler schen Winkel gelten die bekannten Differentialgleichungen: 


otYcosd—ydsind—0.. en ee 2 tage), 
and N IN B), 
Jd+ Jor sin d+(J — J')ö2sin dcos#® = —uH,(l+acoswti)sind ... (4) 


Diese sind von zweiter Ordnung und nicht linear. Für & =0 stellen sie die Bewegungsglei- 
chungen des symmetrischen, schweren Kreisels dar. Das für diesen Fall bekannte allgemeine 
Lösungsverfahren, das auf eine Darstellung von 9, y und Yin elliptischen Funktionen führt, 
wird aber durch das zeitabhängige Glied mit a auf der rechten Seite von (4) unmöglich ge- 
macht. Dagegen lassen sich folgende Verfahren aus der Theorie des schweren Kreisels auch 
auf unseren Fall übertragen: 


a)Überlagerung von kleinen Schwingungen über die reguläre 
Prazessionnch, 
Die partikuläre Lösung von (2)— (4) mit & = 0 welche die reguläre Präzession in H, dar- 
stellt, ist: 
Pp=uwil, y—=@sl, EI re (+) 
mit: 
Van OHREN 17: FE :) F 
(Bestimmungsgleichung für &,; die beiden Wurzeln der quadratischen Gleichung entsprechen 
der „‚langsamen‘‘ und der ‚‚schnellen‘‘ Präzession). 
Man macht nun für eine Nachbarbewegung von (5) mit «=# 0, jedoch a& 1, den Ansatz: 
9 =w,t+ aeimt+ aeie! 
y=Wgt+ beint+ b/eiet huhu Aral): 
dB=d +cei@t + c’elet 
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Die zweiten Glieder der rechten Seite dieses Gleichungssystems bedeuten eine ‚‚freie Nu- 


tation‘‘ mit der Frequenz @,, die dritten Glieder eine ‚‚erzwungene Nutation‘‘ mit der Zwangs- 
frequenz w. Dabei soll 


a<]1 
WAZ ou < w, 
Wb< Ws ob’< ws 
c<d SE, 


sein. 

Einsetzen von (7) in (2) bis (4) und Wegstreichen aller in den kleinen Größen quadratischen 
oder höheren Glieder liefert 3 homogene, lineare Gleichungen für a, b und c, sowie 2 homogene 
und eine inhomogene lineare Gleichung für @’, b’ und c’, Aus der Bedingung für die Lösbarkeit 
des homogenen Gleichungssystems kommt die Frequenz der freien Nutation: 


1 . 
= Wo +(J—2J") w; cos DEN WE SED ern er, 
und als Lösung des homogenen Gleichungssystems: 
DRU 7 © 2.J° cos d, 
OE j’ 9, sin Bar | (9) 
@ „7 21608 8% + J @, cos? d,—J'@; (1 +.cos? u) ler Ah, 
Ögr J’ @3; sind, 
weiter als Lösung des inhomogenen Gleichungssystems: 
ui ad [J w, cos d,-+-J wg cos? du — JS’; (1+cos? d,)] 
Jo (0? — 2) 1 0" 2 0 2 0 
’ . uH, « ’ 5 
b=—1i Ta Yo-+(J—2J') w; cos d,]| mel 
j KH H, sin d, 


Ir al) 


In der hier durchgeführten Näherung stellt sich die Bewegung also als Überlagerung von 
regulärer Präzession, freien Nutationen und erzwungenen Nutationen dar. Die Formeln für die 
freie Nutation sind bereits vom symmetrischen, schweren Kreisel her bekannt. Die Amplitude 
der erzwungenen Nutation geht formal mit Annäherung von w an w; wie 1/(w® — w3) gegen 0; 
es tritt also Resonanzfürw=w, ein (in Wirklichkeit versagt hier natürlich die ange- 
wandte Näherung). Beim Durchgang der Frequenz w durch w; findet ein Phasensprung der 
erzwungenen Nutation um 180° statt. Die 9- und y-Schwingung hat sowohl bei der freien wie 
bei der erzwungenen Nutation 90° Phasenverschiebung gegen die 9-Schwingung. Das Ampli- 
tudenverhältnis der 9- und y-Schwingung ist für die freie und für die erzwungene Nutation das- 
selbe. Die Amplitude der 9-Schwingung richtet sich nach dem Frequenzverhältnis w/w;; im 
Resonanzfall = w; würde auch die d-Schwingung gegen die beiden anderen für freie und 
erzwungene Nutation dasselbe Verhältnis aufweisen. 

Für die Annäherung an den Grenzfall des ‚‚schnellen ee und Ausgang von der 
„langsamen Präzession‘‘, also für 


<wı 
wird: 

Jet | 

a j a O3 

— rzictgd; — er — clgd; ee a a IL 
Fred ri Ad (11) 
5b, TH ee 

ce. sind c’ © sind 


Von den Verhältnissen bei der erzwungenen Nutation der Erdachse unterscheiden sich 
die hier betrachteten Fälle wesentlich, einmal dadurch, daß auf die Erdachse außer einem Dreh- 
moment M, auch noch ein Drehmoment M, gleicher Größenordnung wirkt, zum andern da- 
durch, daß man bei der Erdnutation so weit von der Resonanzstelle entfernt ist (832 : 1 Tag, 
o=27n:19 Jahre), daß von vorn herein mit anderen Vernachlässigungen gerechnet wird. 
Daher rührt es, daß bei der Erdnutation entgegen den Gl. (11) 5’ und c’ von gleicher Größen- 
ordnung sein können. 


u 


. 
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ß) Aufstellung einer Differentialgleichung für #allein 

Da die Diff.gl. (2) und (3) kein zeitabhängiges Glied auf der rechten Seite enthalten, 
können sie genau wie in der Theorie des schweren Kreisels ein mal integriert werden mit dem 
Ergebnis: 

or Vs =- en Fendt Wmcdh - 2.2.0102), 

J'y sin?9+ JO, cos d = (0, = Jw, cos dy+ wa (J cos? 9, + J’ sin? d,)\ (13 

— J 9,08 Bd +Y(J cos? Bd, + J’ sind) I) ) 

Es ist hier absichtlich, auch in der Bezeichnung, zwischen den allgemeinen Anfangswerten %, 

und y, und den zur regulären Präzession gehörigen Werten w, und @,, die alsspezielle 
und y, auftreten können und durch Gl. (6) zusammenhängen, unterschieden. 

Die Integrationskonstanten C/, und (%, haben bekanntlich eine sehr einfache anschauliche 
Bedeutung: C\, ist die Projektion des Drehvektors w, J'C, diejenige des Drehimpulsvektors ® 
auf die Figurenachse, (', die Projektion von ® auf die raumfeste z-Achse. Die vom schweren 
Kreisel her bekannte zeitliche Konstanz dieser beiden Projektionen wird durch das Wechsel- 
feld H, nicht gestört, da dieses weder um die z-Achse noch um die Figurenachse ein Dreh- 
moment ausübt. 

Aus den Gl. (12) und (13) läßt sich @ und » in d und den Integrationskonstanten C, und 
C, ausdrücken: 

. O4 (J cos? 9 + J’ sin? 9) — (, cos d 

er —— I a nee a (14 
? J’ sin? 9 (14), 
ki J' sin? d un 


Speziell für den Kugelkreisel (J”—= J) erhält (14) die zu (15) bemerkenswert symmetrische 
Form: 
. J C re 65 cos ) 
— N. ERLLACN 
u J sin? 9 er 
Setzt man (14) und (15) in (4) ein, so findet man für 9 allein die nichtlineare Diff.gl. 
zweiter Ordnung: 


SÖLRrl@)=—uHll+atoswoi)sind ..:..-..-. . (16) 
mit: 
 J0,,&—JC,cosd 1 .(05— J0ycos d)? cos d 
a Fi sin d nich sin? 9 ws; 


Die allgemeine Lösung von (16) enthält zwei weitere Integrationskonstanten (, und (,. 
Aus ihr ergibt sich mit (14) und (15) auch $ und y und damit die ganze Bewegung des Kreisels. 
Die 4 Integrationskonstanten sind aus den Anfangsbedingungen, nämlich aus den willkürlich 
wählbaren Werten 9, ®, » und $ zur Zeit {—=0 zu bestimmen. Diese erfüllen natürlich im 
allgemeinen nicht die Bedingung (6) der regulären Präzession. 


Mittels (12), (13) und (16) rechnet man übrigens leicht nach, daß 
F (9,) S— — uH, sin d, 
ist, daß also 9 — #,, die reguläre Präzession, eine partikuläre Lösung von (16) für den Fall 
& = 0 darstellt, wie es sein’ muß. 
Von hier aus kann man wieder eine Näherungslösung für «+ 0, jedoch « & 1, gewinnen 
mit dem Ansatz: 
= Ute; Ed, 
Fo) =F(d,)+ er (9), 
womit aus (16) wird: 
J’&+ e[F’(9,) + uH, cos 9] = — a u H,sin 9, cos wt. 
Die Lösung dieser inhomogenen, linearen Diff.gl. führt wieder auf die schon früher angegebene 
Lösung (7) mit (8), (9) und (10). 
Außer dieser Näherung läßt sich noch eine Näherung für 
v<1 
durchführen, welche also alle Bewegungen umfaßt, bei denen die Figurenachse dauernd nur 
wenig von der £-Achse ($,-Feldrichtung) abweicht. 
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y)Die achsennahenBewegungen 
Für den Fall #< 1, Figurenachse wenig von der Feldachse abweichend, lassen sich die 
Funktionen von ® in (16) in eine Potenzreihe nach steigenden Potenzen von 9 entwickeln, die 
wegen dem sin? d im Nenner des zweiten Gliedes von (17) mit 93 beginnt. Da nur ungerade 
Funktionen von # vorkommen, enthält die Potenzreihe nur Glieder mit ungeraden Exponenten. 
Das Glied mit 9-1 verschwindet überdies 'exakt. Die Diff.gl. (16) geht also, falls man Glieder 
bis 9 anschreibt, über in: 


A 


= +(B +alcoswt)d +(D+aE cos wi) d° . a er I 
mit: 
1 
A= u (0 JC}) , 
H, 1 
De ar JO Ce r0) 
H, 
=, ern (19). 
ri ge ) 
Be by) 3094 37907 0 02 + : 
_ uH, 
E Or 


Aus der Form von (18) ist zu ersehen, daß ö nur dann dauernd sehr klein bleiben kann, 
wenn der Koeffizient A entweder auch sehr klein ist oder aber exakt verschwindet. Diese 
beiden Möglichkeiten führen zu verschiedenen Gleichungs- und Bewegungstypen und sollen 
getrennt behandelt werden. In beiden Fällen kann, falls Bnicht sehr klein gegen Dist, was 
für das folgende stets vorausgesetzt werden soll (C und E sind sowieso von gleicher Größen- 
ordnung), das Glied mit 0° vernachlässigt werden. 

Im Fall A=0 wird damit aus (18): 


A n 
= atr(B+oalcoswN)d. a a he Fr. 


Da immer A> 0 ist, wird mit 9— 0 die Beschleunigung 9 positiv sehr groß; die Figurenachse 
wird also sozusagen bei Annäherung an die z-Achse sehr stark von dieser abgestoßen. Ein 
Durchgang durch 9 = 0 ist unmöglich, und falls ® klein bleibt — wofür allein die hier ver- 
suchte Näherung gilt — muß es sich zwischen zwei Librationsgrenzen 9; und 9,,ur bewegen. 
Dies ist offensichtlich nur möglich, wenn B negativ ist, da bei positivem B das ® unbegrenzt 
zeitlich anwachsen würde. Bei negativem B dagegen gibt es für a = 0 eine Gleichgewichtslage 


4 
A b 
4 : U Vs 5 I — 0) Be te rs (21), 


für welche 9 verschwindet. Die Bewegung mit diesem konstanten d, ist wieder die reguläre 
Präzession, diesmal im Spezialfall kleinen 9. Um die Gleichgewichtslage (21) sind nach (20) 
für « = 0 anharmonische, freie Schwingungen möglich; bei «+0 treten durch das cos wit — 
Glied erzwungene Schwingungen hinzu, die sich aber natürlich nicht einfach linear überlagern. 


Die notwendige Bedingung B< 0 läßt sich nach (19) leicht diskutieren. Die Klammer 
im Ausdruck für B ist stets positiv, da ihr Wert zwischen <7 - 5 (Ga — JC,)’und 2 (+ JCh 


liegt. B wird also negativ für alle positiven und einen ER Bereich negativer Werte von 
u H,, d.h. für alle Werte, die der Ungleichung 
7 
u >— (+4 70,04 Pc) 

genügen. Anschaulich bedeutet dies, daß das Drehmoment des Grundfeldes 9, die Figuren- 
achse entweder zur z-Achsehinodernursch wach vonihr weg drängen muß, wenn ® klein 
bleiben soll, wobei die Begrenzung durch die in ©, und (0, steckenden Anfangsbedingungen der 
Bewegung gegeben ist. 

Die anharmonischen, erzwungenen Schwingungen von ® (die durch Hiovaabres von % 
nach (15) zur räumlichen Bewegung der Figurenachse ergänzt werden) kann man wenigstens 


S 
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qualitativ überblicken, wenn man in (20) das d im zeitabhängigen Glied durch einen festen 
Mittelwert Ö ersetzt. Dann entsteht die nichtlineare Schwingungsgleichung mit Zwangskraft: 


ö+ (12105) IR OTEORE BEN. >, un Ars Yhaı" (22), 
Ohne die rechte Seite würde die Lösung dieser Gleichung, die dann auf elliptische Funktionen 
für Ö führt, eine Schwingung ergeben, deren Frequenz mit steigender Amplitude ansteigt. Aus 
der Theorie der erzwungenen, nichtlinearen Schwin- 

gungen weiß man, daß in einem solchen Fall die | Amplitude 

stationäre, erzwungene Schwingung eine nach grö- 
Beren Frequenzen hin geneigte Resonanzkurve nach 
Art des Bild 2 hat. 


Die Resonanzfrequenz in der Grenze unendlich 
kleiner Amplitude (also mit &« — 0) erhält man als: 
34 
ZB nr —=4|B|. 
Ü) 0 


9=% 


Setzt man in dem Ausdruck (19) für B die Werte u H,, 
0, und 0, nach (6), (12) und (13) ein und geht zum Bild 2. Schiefe Resonanzkurve des nichtlinearen 
Limes do —( über, so wird: Schwingungs-Systems 


6renz - Resonanz Frequenz 


1 
WW — z Yo +(J—2J')w;]. 


Dasselbe erhält man übrigens auch aus (8) mit 9,— 0, da für unendlich kleine Amplitude 
auch die achsennahe Bewegung wieder in den Rahmen der früheren Näherung der kleinen 
Schwingungen um die reguläre Präzession fällt. 

Was nun endlich noch die Bedingung A<&1 anlangt, die für die achsennahe Bewegung 
erforderlich ist, so ist sie von selbst erfüllt, falls nur % zu Anfang genügend klein ist. Denn aus 
(12) und (13) folgt für u <1: 


Or (ot) 8, 


Be: ypi on 
(,7.J (9 + Yo) — > Pr (F—J’) Yo d% - 
Daraus nach (19): 
1 2: 
ABIT 


Wenn also nur der Anfangswert 9, von 9 genügend klein ist, so ist für beliebige Anfangs- 
werte 9, und Y, die Bedingung A<1 erfüllt. Den Anfangswerten 9, und y, wird hierdurch 
keine Bedingung auferlegt. Nur muß 9, um so kleiner gewählt werden, je größer 9, und Y, sind. 


Dune peazialfallider Machıeuschen’Diffierentialgleichung 
Die im Abschnitt y) behandelte Bewegung ist von der des Abschnitts a) nicht allzu 
verschieden. Sie stellt ihr gegenüber eine Erweiterung auf nichtlineare Schwingungen von ® 
um eine Gleichgewichtslage dar (allerdings nur näherungsweise, insoweit im zeitabhängigen 
Glied der Differentialgleichung 9 durch # ersetzt werden darf) unter gleichzeitiger Spezialisie- 
rung auf kleine 9. Einen gänzlich anderen Bewegungstyp erhält man jedoch, wenn A exakt 
gleich Null gemacht wird. 
Zunächst soll untersucht werden, was die Bedingung 


A=0 
anschaulich bedeutet. Da , nicht Null sein soll, erhält man aus (23) mit A = 0 sofort: 


He ae ng ve) Me EEE 27 

Yo J 
Es muß also, damit A verschwindet, die Drehgeschwindigkeit des Kreisels um die Figurenachse 
zur Drehgeschwindigkeit der Figurenachse um die 2-Achse zu Anfang der Bewegung in einem 
ganz bestimmten, durch (24) gegebenen Verhältnis stehen. Dieses Verhältnis bleibt, da ja ©, 
und (, Konstanten der allgemeinen Bewegung sind, auch während der ganzen Bewegung 
erhalten, solange die Bedingung 9< 1 bestehen bleibt. Da J’/J je nach der Form des Kreisels 
alle Werte zwischen den Grenzen 1, (unendlich dünne Kreisplatte senkrecht zur Figuren- 
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achse) und & (unendlich dünner Stab in der Figurenachse) annehmen kann, können für das 
Verhältnis (24) alle Werte zwischen O0 und ® vorkommen, und zwar 0 bis 1 beim abgeplatteten, 
1 bis oo beim gestreckten Kreisel; es muß jedoch stets @, dasselbe Vorzeichen haben wie %,. 

‚Setzt man nun A=0in die Diff.gl. (18) für d ein und vernachlässigt wieder das Glied 
mit 0% und höhere, so folgt diesmal die Mathieu sche Differentialgleichung: 


d=(B+olcosot)#d ..... ME . (25) 
wo Bund C durch (19) gegeben sind, hier wegen A = 0 aber (0, — 10, Hedsor werdäi kann 
und Bssich dadurch zu: 2 

NET eh 0? (26) 
y ae N a az 


vereinfacht. 

Die Lösung von (25) unterscheidet sich grundlegend von derjenigen von (20). Mit A=0 
ist die „‚Abstoßung‘‘ der Figurenachse durch die z-Achse verschwunden, und 9 kann ohne 
weiteres durch Null gehen und sein Vorzeichen wechseln. Die Zeitabhängigkeit von ® läßt sich 
damit auch nicht mehr näherungsweise als Schwingung um einen 
Gleichgewichtswert d, auffassen. 

Betrachtet man zunächst die Gl. (25) für «=0, also die 
achsennahen Bewegungen des symmetrischen, schweren Kreisels, 
so folgt für 9 einfach die lineare Schwingungsgleichung 

9 Be 
die für negatives B Schwingungen mit der Frequenz 
ergibt. 
Die Hinzunahme von y nach (15) liefert die verschiedenen 


Typen a)—.e) (Bild 3) der Bahn der Kreiselspitze je nach der Größe 
von B, d.h. von u. A,. 


Der Fall b) entspricht der regulären Präzession des kräfte- 
> freien, der Fall d) der re- 
N .. .. ” 
c gulären Präzession des 
D c B schweren Kreisels in der 
20 | 21a $ d Grenze = 0 (für 90 


Zn 


würde d konstant bleiben 
und ® daher, ähnlich wie 
im Fall e), den Bereich 
ISIS kleiner Werte verlassen), 
Si T A der Fall e) ergibt positives 
Na B und damit Entfernung 
N SS der Figurenachse von der 
\ x z-Achse auf einer Spiral- 
bahn. 


er 


S 
| 
j 


7 
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Bei #0 ist die 
Mathieu sche Gl. (25) 


/ N \ N 
Ei 0) \ NN zu lösen und die Bahnen 
-75 / NNEAN der Kreiselspitze werden 
7 1/ \ schon sehr kompliziert. 
-20 Man kann jedoch unter- 
- U Er, 10 ge suchen, unter welchen Be- 


Bild 3. Achsennahe Bahnen der Bild4. Darstellung inderA— -EbenederMathieu- dingungen ® ständig im 

Kreiselspitze füra = 0 schen Gleichung Bereich kleiner Werte ver- 
bleibt. Dies ist offenbar der Fall in den sogenannten ‚‚stabilen Bereichen‘ der Konstanten 
der Mathieuschen Gleichung, während in den ‚instabilen Bereichen‘ # nach GI. (25) im 
Laufe der Zeit über alle Grenzen wächst, in Wirklichkeit also aus dem Bereich der hier durch- 
geführten Näherung entweicht. Man kann daher die instabilen Bereiche als ‚‚Resonanzbe- 
reiche‘‘ bezeichnen, wenn man unter Resonanz versteht, daß bei noch so kleinem (aber von 
Null verschiedenem) Anfangswert 9, das d durch die Einwirkung der Zwangsfrequenz & aus 
dem Bereich kleiner Werte hinausgeführt wird. 

Schreibt man die Mathieusche Gleichung in der üblichen Form: 


FYyy adhreos2e)y=0 28 
Fe cos 2.2) y= 0 uni AS ER, 


Ü 
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die mit (25) durch die Beziehungen: 


De ae Brauer a0 2 

Val, =yzut, ar But eg la aa: (29) 
verknüpft ist, so stellen sich die stabilen und instabilen Bereiche in einer h2-A-Ebene ent- 
sprechend Bild 4 dar, wo die stabilen Bereiche schraffiert angegeben sind !). 

Läßt man nun bei vorgegebenem B und «C' in Gl. (25) die Zwangsfrequenz w alle Werte 
von O bis oo durchlaufen, so entspricht dies wegen (29) im Bild 4 einem Punkt, der auf einer 
Geraden AO, bzw. BO u. s. w. aus dem Unendlichen bis zum Ursprung die h2-A-Ebene durch- 
läuft. Bringt man auf dieser Geraden nach den Beziehungen (29) eine w-Skala an, so läßt sich 
unmittelbar aus der Figur ablesen, für welche w-Werte ‚‚Resonanz‘‘ auftritt, nämlich für alle 
diejenigen, für die der Punkt in der h?-A-Ebene in ein instabiles, nicht schraffiertes Gebiet fällt. 
Es gibt also hier keine Resonanz stellen, sondern Resonanz bänder, die sich in unend- 
licher Zahl bei niedrigen Frequenzen häufen. Dabei sind je nach der Lage der Geraden, die 
durch das Verhältnis & C/B bestimmt ist, folgende Fälle zu unterscheiden: 


a) B negativ; |«C|< |B| schwaches Wechselfeld; Gerade OA 

-b) B negativ; |«C| und |B| von gleicher Größenordnung; stärkeres Wechselfeld; 
Gerade OB 

c) B=0; uH,= — J?0?/4 J’ Gerade OC 

d) B positiv; |«C] >|B| stärkeres Wechselfeld; Gerade OD 

e) B positiv; Ja C|<|B| schwaches Wechselfeld; Gerade OE. 


Im Fall a) befindet sich der Bildpunkt fast überall in stabilen Bereichen mit Ausnahme 
ganz schmaler Bänder bei etwa 
Ne 1A, 8, nr 
also bei 


yı ege,.e 


@®—4|B|, |Bl, S |B 
oder bei 


2 2, 
= * „er * — * N ee a IT De Te GT 
w—=2nt, 0%, 2 0%, ...— w*... (30) 


Es gibt also hier Resonanz bei der doppelten d-Frequenz der Bewegung ohne Wechsel- 
feld, sowie bei allen ganzzahligen Teilen davon. Es treten sogenannte „subharmonische 
Resonanzen‘ auf. Jede Resonanz ist nicht völlig scharf, sondern ein Resonanz band. 
Die Breite der Bänder ist jedoch bei |«C|< |B| äußerst gering, insbesondere die der höheren, 
da die Grenzkurven in Bild 4 sich bei größerem / in immer höherer Ordnung berühren. 

Im Fall b) wechseln breitere Resonanzbänder mit breiteren Lücken ab, wobei die 
genauen Verhältnisse sich nach der Steigung der Geraden, also nach dem Verhältnis der Größe 
«C und Brichten. 

Im Fall ec), der in bezug auf die Größe u H, des konstanten Richtmomentes dem Fall d) 
des Bildes 3 entspricht, läuft die Gerade OC, hier die h2-Achse, ganz überwiegend durch nicht 
stabile Gebiete. Hier sind also, umgekehrt wie bei a) breite Resonanzbänder durch schmale 
Lücken getrennt. Nur in engen Frequenzbereichen bleibt # klein. Die Fälle a) bis c) haben 
jedoch gemeinsam, daß die Umgebung des Ursprungs, also die ganz gro ßen Frequenzen, zu 
diesen stabilen Frequenzbereichen gehören. Es gibt also bei negativem B und beliebigem « C 
stets eine Grenzfrequenz, oberhalb deren keine Resonanzstellen mehr liegen. 

Anders ist es im Fall d), bei positivem B. Da verläuft die Gerade in der Nähe des Ur- 
sprungs im instabilen Gebiet. Der ganze Frequenzbereich von O bis oo besteht hier aus Reson- 
nanzbändern, die nur durch ganz schmale stabile Zwischenräume getrennt sind. Nur in der 
Nähe ganz bestimmter Frequenzen tritt also hier keine Resonanz ein, bleibt 9 klein. Daß es 
überhaupt solche Frequenzen gibt, ist besonders bemerkenswert. Denn ohne Wechselfeld war 
ja der Fall positiven B (Fall e) des Bild 3 durchweg instabil. Hier stabilisiert also das Wechsel- 
feld eine an sich instabile Bewegung, ein paradox anmutendes Ergebnis, das aber eine gewisse 
Analogie in dem Pendel mit vertikal schwingendem Drehpunkt hat, das bekanntlich unter 
gewissen Bedingungen auch mit dem Schwerpunkt über dem Drehpunkt (also in einer an 
sich instabilen Lage), stabil wird. 

Fall e) schließlich unterscheidet sich von Fall d) nur durch die geringere Steigung der 
Geraden, die jetzt überhaupt kein stabiles Gebiet mehr durchschneidet, weil dieGrenzkurven 

2 
in Bild 4 in ihren nach links oben auslaufenden Ästen asymptotisch die Steigung En > ns 
1) Siehe z.B. M. J.O.Strutt in „Erg. d. Math.“ 1. Bd., Heft 3, S. 229 ff., 1932. 
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annehmen. Es gibt also bei positivem B und zu schwachem Wechselfeld überhaupt keine im 
Bereich kleiner # stabile Bewegung. Dies erscheint unmittelbar plausibel, da ja bei positivem B 
erst das Wechselfeld die Bewegung stabilisiert und dies offenbar unterhalb einer gewissen 
Stärke nicht mehr vermag. 


Fall II. Wechselfeld senkrecht zum Grundfeld 


Der Fall, in dem das Wechselfeld 9, senkrecht zum Grundfeld 9, steht, soll von vorn- 
herein auf den ‚‚schnellen‘‘ Kreisel spezialisiert werden, für den y < p ist und damit der Unter- 
schied zwischen Impulsachse und Figurenachse nicht berücksichtigt zu werden braucht. Das 


Drehmoment 
M— ulad] 


des Gesamtfeldes 9 auf die Figurenachse (a Einheitsvektor in deren Richtung) kann dann 
geschrieben werden: 
@z 


ne a | 
wo ® das Impulsmoment und P sein Betrag ist, da für , < 


P=Jo, 
und nach (6): 


ist. 
Damit wird die Kreisel-Grundgleichung R- M zu der linearen, homogenen Vektor- 
Differentialgleichung erster Ordnung für ®: 
Ri 


Ba all ee ent 


Da B stets senkrecht zu ® steht, bleibt der Betrag P des Impulsmomentes konstant. @; ist 
die Präzessionsfrequenz im konstanten Feld $,. Das Gesamtfeld 9 besteht jedoch jetzt aus 
den beiden zueinander senkrechten Komponenten 9, und $,, von denen 9, eine harmonische 
Zeitfunktion mit der Frequenz @ ist. 

Im allgemeinsten Fall kann hierbei die Spitze von 9, eine Ellipse durchlaufen (ellip- 
tischpolarisiertes Wechselfeld). Legt man die x- und y-Achse in die beiden 
Achsenrichtungen der Ellipse (2-Achse in Richtung von $,), so stellt sich das Feld 9 dar als: 

H,=H,cos wt, Hahn | 
H,=Hisin ot, Hr=add 
wen | 
und die Differentialgleichung (31) wird in Komponenten: 
P,= w(— P,+ o P,sin wi) | 
P,=w(P,—«aP, cos wt) | a NG (33). 


P, = w,(—a' P,sin wt + P, cos wt) 


. (32) 


” 


Dazu kommt noch wegen der Konstanz des Betrags P: 
P% = Ps + Pr P=constul, „ia Aigen 


was auch formal aus (33) durch Multiplikation der drei Gleichungen mit P,, P,, P, und Ad- 
dieren-hervorgeht. 

Am einfachsten zu behandeln ist das Gleichungssystem (33) im zirkular polarisierten Fall 
&=.a. Außerdem soll dann noch der linear polarisierte Fall a —=0 betrachtet werden, der 
schon alle Besonderheiten des allgemeinen, elliptisch polarisierten Falles zeigt. 


o) Zirkular polarisiertes Wechselfelda=a 
Hier läßt sich die allgemeine Lösung von (33) mit 3 willkürlichen Konstanten angeben, 
und zwar am einfachsten, indem man das Gleichungssystem auf ein mit & um die z-Achse 
rotierendes Koordinatensystem transformiert, das für £ = 0 mit dem ruhenden übereinstimmt. 
Dann fällt in den neuen Koordinaten «', y', 2’ die Zeitabhängigkeit auf der rechten Seite weg: 


Pr=(w—o)Py,: : Pr=—lw— wo) Pr—0amPr,.. Pr=omPy : (85) 
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und man erhält die allgemeine Lösung dieses linearen Diff.gl.-Systems mit konstanten Koeffi- 
zienten: 


W—Wy . w—W 
Pr = 0, — sind — 0, cos At+ Rn 
Py = C,cosAt+ (0, sin At BE MEERES 90) 
or 70) 0 — 
Pr = 01T sin dt — 0, z wsit—Q; n = 
mit: 
EN er 0 0 EI ee 2a 9 (37): 


Die 3 Integrationskonstanten (C\/, C,, ©, können leicht in den Anfangswerten Pzo, Pyo, Pzo 
ausgedrückt werden. Durch Rücktransformieren von (36) in das ruhende System x yz erhält 
man: 


P,= > cos wtsinAt — sin wteosäi) 
Sees wt cos At-+ sin otsinAr)+ Gr cos wt, 

1 0, (cos @tcosit+ = sin otsinät) (38). 
+0,(eos wisin At — sin wteosär)+ 0, sin wt, 

P,= 0, sinit— 6, — hl, 


Hieraus ergeben sich folgende spezielle Formen: 
Mit«—=0 und damiti=w— wo, diefreie, reguläre Präzession im Feld $,, da dann : 
cos wtsin At — sin wt cos At= — sin w,t wird u. s. w.: 
P, = — (sin &t — (%, cos @st, P, = Cı c0s wgt — (5 sin Wet, P,=-—-(; (39) 
Mit «<1 (Glieder mit &® und höhere werden vernachlässigt; sehr schwaches Wechsel- 
feld 9): 
@®; 


— (,cos wt, 
—=.(s - 


P, = — Cı sin wgt — (0%, cos Wt+& 
2) 


[OP F 
mu C, sin wt, er (AN), 


2 


P,= (1 608 wt — (, sin Wet 2, 


wg 
w— Wg 


P,=— (,+a [— 02 005 ( — @) t+ Cı sin (o— @,)i] 


Setzt man hierin C, und (, =(, so kann man die verbleibende Bewegung als reine, 
erzwungene Präzession betrachten: 


7) m) 
P,= (eos WE: P,= 2 El sm @E, Bar (Fun rlal), 
© an 3 Y no, 3 2 3 (41) 


Dies ist übrigens, wie man leicht durch Einsetzen sieht, sogar eine strenge Lösung von (33) 
für =. Die allein hierin noch vorkommende Konstante ©’, bestimmt sich mittels (34) aus 
dem vorgegebenen Drehimpuls P. 

Die Bewegung (41), läßt sich, was P, und P, anlangt, als lineare Überlagerung der freien 
Präzession (39) mit der Frequenz w, und der erzwungenen Präzession (41) mit der Frequenz 
auffassen. Die Verhältnisse sind weitgehend analog denen bei der Überlagerung freier und 
erzwungener Schwingungen. Für ® = w, geht die Amplitude der erzwungenen Präzession 
gegen oo (die hier durchgeführte Näherung versagt dann natürlich). EstrittalıoResonanz 
mit der Frequenz der freien Präzession auf, eine Erscheinung, die von der Theorie der Atom- 
kerne im magnetischen Hochfrequenzfeld her bekannt und dort gründlich untersucht ist ?). 
Der Resonanznenner erscheint allerdings in der etwas ungewöhnlichen Form (w — w,) statt 
(© — 2); der Grund hierfür wird später deutlich werden. 

Zur Veranschaulichung der allgemeinen Lösung (38) verwendet man zweckmäßig die 
Form (36), die im rotierenden K.S. #’ y’ z’ gilt. Führt man hierin die komplexe Größe 

2)z2.B.F. Bloch, Phys. Rev. 70, 460 (1946). 

14* 
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NY = Es ein, so geht (36) über in: 


A & 3 


Perl eit) 1. Pe=— 7 Ueirtt — O, 


sl 


ren — — 2 ei 0, 


oder vektoriell: 


(42) 


P= + Pı- 

Aus den Komponenten sieht man, daß ®, stets senkrecht zu ®, steht, sowie daß der Betrag von 
P, konstant und gleich |A| ist (wegen (37)), und daß die Amplituden von ®, in der Ebene x’ y’ 
und senkrecht zu ihr gleich groß und um 90° gegeneinander phasenverschoben sind. Der 
Endpunkt von ®B läuft also ineiner Ebene, die zudem system- 
festen Vektor B.ssenkrecht; stehts und’ zwar beschreihtucre, 
dieser Ebene mit der .Winkelge&chwindiıakeıt/i einensakreıs 
vom Radius |\umden Endpunktvon %, (Bild 5a). Der Vektor ® selbst, in 
unserer Näherung des schnellen 
Kreisels also die Figurenachse, 
beschreibt einen entsprechen- 
den Kreiskegel. Dies alles be- 
zieht sich aber auf das rotie- 
rende System &’ y’ 2’ und man 
muß also das ganze Vektor- 
gerüst noch zusätzlich mit der 
Winkelgeschwindigkeit ® um 
die 2’-Achse rotieren lassen um 
die Bewegung im raumfesten 
K.-S. 2 yz zu erhalten. 

Dabei schwankt offenbar 
der Winkel 9 von ® gegen die 
z-Achse zwischen zwei Libra- 
tionsgrenzen 9, und d,, die sich 
aus der Lage des B-Kegels ge- 
gen die 2-Achse ergeben. Der 
„Resonanzfall‘“‘, ® = w,, zeich- 
net sich jetzt anschaulich allein 
dadurch aus, daß der — in 
& y' 2’ — systemfeste Vektor ®, 
in die ©-Achse fällt, der Vektor 
%, in einer Ebene parallel der 
y' z'-Ebene rotiert, und daß da- 
mit 9, = nr — d, wird (Bild 5b). 
Fällt im Resonanzfall außerdem 
®, die Figurenachse, zu irgend 
einer Zeit (etwa zu Anfang, 
t = O)in die y’ z’-Ebene, so muß 
C; —=0 sein und der P-Kegel 
artet in die y’ z’-Ebene selbst 
aus. (Bild 5c). Dies ist der 

DE ' i einzige Fall, in dem der 

Bild 5. EEE ED et ME 5 im rotierenden Sklınlaal ” alle Werte zwischen 

0 und z erreicht, die Kreisel- 

spitze im raumfesten System & yz also eine Spiralbahn durchläuft, die sich über die ganze 
Einheitskugel erstreckt. 

Der einfachste Fall andererseits ist der, wenn in (42) die komplexe Amplitude Y=0, 
also ®, = 0 wird. Dann liegt im #’ y’ 2’-System alles fest: 90-Richtung (z-Achse), Figuren- 
achse (Richtung B— %,) und 9,-Richtung (#’-Achse) liegen in einer Ebene (x z’-Ebene), 
die im Raum mit der Winkelgeschwindigkeit um die z-Achse rotiert. Dies ist die von 
Bloch?) besonders herausgestellte einfachste Lösung des Gleichungssystems (33) für &’ = a. 
Sie stellt gleichzeitig die anläßlich unserer Gl. (41) als ‚‚reine, erzwungene Präzession‘‘ bezeich- 
nete Bewegung des Kreisels dar. 


ß) Linear polarisiertes Wechselfeld „=0 


Der Fall des linear polarisierten Wechselfeldes ist wesentlich schwieriger zu behandeln als 
der zirkular polarisierte, da sich hier kein Mittel findet, die explizite Zeitabhängigkeit aus den 


() 
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Differentialgleichungen zu entfernen. Aus diesem Grund wird bei der magnetischen Kern- 
resonanz fast stets mit dem zirkular polarisierten Feld gerechnet, obwohl experimentell ein 
linear polarisiertes angewandt wird. Bloch und Siegert?) haben übrigens auf wellen- 
mechanischer Grundlage die Gleichwertigkeit beider Felder im Resonanzfall nach- 
gewiesen, ein Ergebnis, das mit den nachfolgenden, mit klassischen Methoden erhaltenen 
Resultaten in Übereinstimmung steht. 


Schwingt das Feld 9, in der @-Achse, ist also a’ =0, so kommt an Stelle von (35) das 
Gleichungssystem: 


P,=—wP,: P,=w(P,—aP,cos ot), P,—= w,aP, cos wt E22 (49): 


Man kann zunächst versuchen, durch Weiterdifferenzieren und Eliminieren von P, und 
P, und deren Ableitungen eine Differentialgleichung für P, allein zu erhalten. Man muß 
hierbei die erste der Gl. (43)zweimal, die beiden anderen je einmal weiter differenzieren und erhält 
die lineare, homogene Differentialgleichung 3. Ordnung für P,: 


P,+ wtgwt P,+ @%(1+ a2 cos? ot) P,+ wiotgwtP,—0. eh, (AA), 


Verwendet man bei der Elimination die Beziehung (34), die bereits eine erste Integration 
mit der Integrationskonstanten P darstellt, so kommt man mit der 2. Ordnung aus; die Diffe- 
rentialgleichung wird aber dafür nicht linear. Sie lautet 


EBEN lila) ©a5); 


9 
© 


P,+ 2 P,= aw} cos our Sep 


In beiden Fällen kommt man nicht weiter. Weder für Gl. (44) noch für Gl. (45) läßt sich 
eine allgemeine Lösung angeben. Auch gehört keine von beiden Gleichungen einem bereits gut 
durchdiskutierten Typ an, aus dem sich allgemeine Ergebnisse ableiten ließen. 


Es bleibt also nichts übrig, als wenigstens eine Näherungslösung für «<1 zu suchen, 
etwain Form einer Potenzreihe nach Potenzen von a. Man kann hierfür etwa mit dem Ansatz: 


PR=mtemt®pt--, P,=9ptaeıt®g+t'', P,=rbartort (46) 


in das Gleichungssystem (43) eingehen und erhält zu den einzelnen Potenzen von «a folgende 
Gleichungssysteme: 


M= — %% 

0°: g=WgPo 2 EEE EEE TER U (47), 
Be 

BER ED Et COSDE 2 0. 2 20 (48) 


| 11 = W9 C0S wit 
usw. 

Das erste System (47) ist homogen und hat als Lösung die ‚‚freie Präzession‘‘ mit 3 will- 
kürlichen Konstanten, genau entsprechend Lösung (39). 

Die weiteren Systeme sind inhomogen und ergeben nach Einsetzen der je vorhergehenden 
Stufe als Lösungen eine Lösung der homogenen Gleichung, die identisch mit der von (47) ist 
und mit dieser zusammengezogen, d. h. weggelassen werden kann, plus einem ‚‚erzwungenen‘‘ 
Glied ohne weitere willkürliche Konstante. Die Lösung von (48) wird z. B.: (nur erzwungenes 
Glied): 


2 


92 
—= ——— (5 cos wt 
pı w: — 3 3 ’ 
27) 
= -— (5sin wt, 
Eurer Ne (A9), 
le h j 
rn = (1 ——— (wsin @t cos @gt — @, cos wLsin @zt) 
wo? — w3 
@s E B 
— (5, —— (w sin wt sin @gt + @, cos @E COS @zt) 
02 — 


2 


®)F,BlochundA.Siegert, Phys. Rev. 57, 522 (1940). 
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Die C, ©, C, sind dabei die Konstanten der Lösung von (47), d. h. der Gl. (39). Fügt man 
aus (39) und (49) eine „‚„Lösung erster Ordnung‘‘ zusammen, d. h. bis zur Potenz «! unter Ver- 
nachlässigung höherer Glieder, so wird diese 


2 


10) 
P,= — 0, sin @t — (0, cos wt-+ a isn O, cos wt, 
X) 5 
P,= 01 cos wt — (,sin &t+& 0 C, sin wt nn. (0). 
w* — w2 
P,= — (+ [4 Glieder nach (49)] 


Dies ist für das linear polarisierte Feld das Entsprechende zu den Gl. (40), die für das 
zirkular polarisierte Feld gelten (be ide nur näherungsweise für « < 1). Die Ähnlichkeit zwi- 
schen (40) und (50) ist sehr groß. Nur tritt jetzt in (50) der ‚‚richtige‘‘ Resonanznenner (u —3) 
auf, und außerdem sind die Amplituden der erzwungenen Glieder von P,und P,verschie- 
den, während sie beim zirkular polarisierten Feld gleich waren. Während also beim zir- 
kular polarisierten Feld die reine erzwungene Präzession (C, = C, —0) auch zirkular um die 
z-Achse ‚‚polarisiert‘‘ ist, ist sie beim linear polarisierten Feld „elliptisch polari- 
siert‘ Dieser Unterschied verschwindet aber mit Annäherung an den Resonanzfall, da für 
© — &, zwar beide Koeffizienten oo werden, der Grenzwert ihres Verhältnisses aber 1. 


Um die Entwicklung bis zu höheren Potenzen von « zu erhalten, kann man nach (48) und 
den entsprechenden folgenden Gleichungssystemen allgemein schreiben: 


nach der jeweils ersten Gleichung: „= —— Dn 
@® 
nach der jeweils dritten: Yn+1ı = W@ COS WE = — P, Cos wit, 
Ya+ı = — [ pn cos wit dt. 


Dieses in die zweite Gleichung, und sie wiederum in die weiterdifferenzierte erste eingesetzt, 
ergibt: . 
Pn+2+ © Pn+2 = WErn+ı COS wi 
— — 2 cos wt / m cos wtd=fu(t) 
mit der Lösung: 


Pr = OHMCOs N | DR cos wid ‚1e Fre al 
woO einen Integraloperator bedeutet, der definiert ist durch: 
Oy= wzlcos @t [ ysin wt dt — sin wgt | y cos 0 Klarheit. er 5, 


(51) stellt zusammen mit (52) eine Rekursionsformel für die p, dar, aufbauend auf p, und p.. 
Das Verfahren läßt sich bis zu beliebig hohen p, durchführen, da unter den Integralen stets nur 
Potenzsummen von harmonischen Funktionen vorkommen. Die Formeln werden aber rasch 
sehr umfangreich. Aus den p, lassen sich nach (48) usw. die 9, und r„—ı durch Differentiation 
gewinnen. Ob das Verfahren konvergiert, erscheint fraglich. Mindestens lassen sich aber die so 
erhaltenen Reihen der Form (46) für genügend kleines a als semikonvergente Reihen verwenden. 


Etwas einfacher wird die Lösung, wenn man von vornherein nur die reine erzwungene 
Präzession berechnet, also (, = (0, —=0 setzt. Dann verschwinden nach (51) wegen p, = 0 
alle 9, und gq, mit geradem n und alle r, mit ungeradem n. Nach Umrechnung der cos- und 
sin-Potenzen in Winkelfunktionen vielfacher Argumente sieht die Lösung so aus: 

1 = c{ cos wt [Aua+ A,303+ A,50® + gt -)+ cos 3 wt[A,0°+ A;0®+ ir -] 
| +cos5wt[A,0d + -]+ +}, 
Den 3 
Pi 12 sin wt[A1ı1a+ Aysad+ Asad-+--]+ 2 sin 3@t[A30°+ Asa +] a 
Du ae (53) 
- 9 ins ot [A,0d+- S + u he 
©, J 


P,= C{1+cos2wt[A3a0?+ Ayal+-- -]+ cos 4wtlAuat+:)+:} 
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mit folgenden Werten der Koeffizienten 


ie [D} wo; 
eg ER Rn 
a of 1. 0 
8.(0? — w3)? ' #78 (02 — w2)2 (9. — w2) ’ 
nee) Bw —w) 
k Be a ee nenne 
nen MN N EN (54). 
re [07 7,19 3 8 
8 (02 — 3) (9 0 — 2) “732 (02 — w2) (9.0 — w?) ’ 
TEE in Eis END Ban BER 22 
64 (0? — 3) (9 w® — w})? ae 308 Fr: 
BR ANA bau nl 64 (0? — w2) (I. m — w2) (25 w® — @?) 


Der Ausdruck (53) mit den Koeffizientenwerten (54) ist allerdings nicht die einzige mögliche 
Form der Lösung. Man könnte nämlich die einzige willkürliche Konstante ( selbst als Potenz- 
reihe in a schreiben: 


ee) nd), 


in (53) ausmultiplizieren und neu nach Potenzen von a ordnen. Die k; wären (beliebig) vorzu- 
geben; C', wäre jetzt die willkürliche Konstante der Lösung. Selbstverständlich würden dabei 
alle Koeffizienten A,, anders, je nach Wahl der %,. Diese Willkür der Darstellung rührt daher, 
daß man bei der Berechnung der einzelnen Stufen nach (51) in jeder Stufe eine neue Integra- 
tionskonstante erhält. Die Lösung hat aber natürlich nicht unendlich viele willkürliche Kon- 
stanten, sondern (bei der ‚‚reinen erzwungenen‘‘ Präzession) eine, nämlich in der allge- 
meineren Darstellung C,, die etwa durch den Betrag P des Impulsmomentes festgelegt ist. 
Bildet man nämlich aus (53) P? = P};-+ P,-+ Pi, so verschwindet dabei tatsächlich, 
wie es sein muß, inje der «-Potenz die Zeitabhängigkeit, und es bleibt für P? eine Potenzreihe 
in « mit zeitunabhängigen Koeffizienten. Wollte man P selbst als willkürliche Konstante der 
Lösung einführen, so müßte man die %k, in (55) so bestimmen, daß alle Koeffizienten dieser 
Potenzreihe für P Null werden mit Ausnahme des Absolutgliedes. In der Form (53) der Lösung 
hat dagegen die willkürliche Konstante C die Bedeutung des zeitlichen Mittelwertes von P,. 
Besonders interessant ist nun, daß in den Koeffizienten A;, außer dem Resonanznenner 


(w® — 2) auch die weiteren Resonanznenner (90? — 2), (250 — 3) .... auftreten. Es gibt 
also, und zwar in immer höherer Ordnung in «, auch für 
= @/3, W/5 '**- Wwf2k+1 


unendliches Anwachsen der Amplitude (in dieser Näherung); d.h. es treten auch hier sub - 
harmonische Resonanzenauf. Rechnet man mit von Null verschiedenen Werten 
von C, und (,, so erscheinen auch die Resonanzen bei geraden Bruchteilen von @;. 

Bricht man die Entwicklung (53) nach den Gliedern mit « ab, so kommt, wie es sein muß, 
wieder die Näherungslösung (50) mit C, = (, —0. 


y) Elliptisch polarisiertes Wechselfeld a+#«=+0 
Die Verhältnisse liegen hier grundsätzlich ähnlich wie beim linear polarisierten Wechsel- 
feld. Es soll daher nur die Näherungslösung für a< 1, a’< 1 bis zu den in a und a’ linearen 
Gliedern angeschrieben werden. Man erhält sie einfach durch Superposition (die in der linearen 
Näherung noch erlaubt ist) zweier Lösungen der Art (50) je mit a und a’ geschrieben, wobei aber 
bei a’ P, mit P, vertauscht und eine Phasenverschiebung von 90° eingeführt wird. Das Er- 
gebnis ist, wenn man die Konstanten ebenso benennt wie in (50): 


E 12) 
P, = — C, sin ot — (©, cos wt + CO, (ww; + & w) Pr cos ot | 


. (56). 


P, = Ci c0s &t — (, sin wt-+ 0,(aw-+ a’ @,) ne sin wt,. | 
il! 
P, = (;+ a[4 Glieder] + a’[4 Glieder] 


In dieser Lösung muß sowohl diejenige für zirkulare wie auch die für lineare Polarisation 
enthalten sein. Tatsächlich kommt man mit a’ = 0 auf (50) und mit a’— « auf (40) zurück. 
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Von hier aus erscheint auch der Resonanznenner (o — @,) in (40) verständlich, da aus dem 
ursprünglichen (w® — 2) vom Zähler her (w + @,) fortgekürzt wird. 

Die Entwicklung zu höheren Potenzen von & kann genau wie beim linear polarisierten 
Wechselfeld erfolgen, wobei wie dort in den höheren Gliedern subharmonische Resonanzen auf- 
treten. Nur werden die Formeln im einzelnen noch wesentlich umständlicher und unüber- 
sichtlicher. 


Eingegangen am 7. Juni 1952. 


Durchbiegung von geschichteten Plattenkörpern 
mit leichter Füllung 
Von V. Vodicka in Pilsen 


Durchbiegung von geschichteten Platten ist bisher nur für unendliche Platten berechnet worden. Im 
Falle endlicher Abmessungen betrachtete Prusakov das Problem ohne die einschränkenden Voraus- 
setzungen, daß die innere Zentralschicht undeformiert bleibt und die Spannungen gleichförmig über die 
Höhe der äußeren Schichten verteilt sind, jedoch immer nur für den einfachsten Spezialfall eines aus 
3 Schichten bestehenden Plattenkörpers. 

Dieser Aufsatz behält die Prusakovschen Verallgemeinerungen und fügt darüber hinaus noch 
andere hinzu. 


The deflection of sandwich plates was, till now, thoroughly investigated for plates of infinite dimensions 
only. For finite plates, Prusakow considered the problem without the usual restrictions, namely thal 
the inner layer is not deformed and that the stresses are equally distributed in the outer layers, but for the 
special case of a plate consisting of three layers. T'he present paper goes far beyond Prusakow’s 
generalizations. 


La deflexion de plaques placees en couches les unes sur les autres jusqu’&a present n’a EiE calculee que 
pour des plaques infinies. Au cas de plaques finies Prusakov aconsidere le probleme sans les suppo- 
sitions restrictives que la couche centrale interieure ne soit pas deformee et les tensions soient distribuees 
egalement dans les couches exterieures, mais seulement pour le cas special le plus simple d’un corps de 
plague compose de trois couches. 

Cet article garde les geneeralisations de Prusakov et au delü ajoute encore d’autres. 


Pacyer nmporuöa pacclIoeHHbIX INACTHH ÖBLI AO CHX HOP BO3MO3KEH JIHIIB IA ÖECKOHEYHBIX 
INacTuH. B CIybae KOHEYHEIX Pa3MepoB NFOÖJIEMA paccMaTpuBanach [IpycakoBbIM bes OTPaHH- 
YHBAINHX NPeANOJOKeHHha O TOM, YTO BHYTPeHHHH HeHTPaAbHBIa CIOH OcTaeTca ea Nedop- 
MAalHH, a HAIpASKeHHA pacıpenenamTca PABHOMEePHO HA BEICOTE BHEIIHUX CIIOeB, — BCETNA, 
OAHAKO, B IIPOCTeÄIIeM YACTHOM CAYY&e KOPNYCA INIACTHH, COCTOAMEeTO M3 TpexX coeB. B 
ES paöore, coxpansoımei o6boömenna Ilpycakosa, mpensaraertca pa AaubHeinux 
oDoÖImeHMi. 


In Prikladnaja matematika i mechanika, Tom XV, 1951, S. 26—35, befaßt sich Herr 
Prusakov mit der Durchbiegung einer aus 3 Schichten bestehenden ebenen Platte und 
zwar unter sehr allgemeinen Bedingungen. Jedoch sind seine Ausführungen noch mancherlei 
interessanter und weitgehender Verallgemeinerun- 
gen fähig, wovon hier eine ausgeführt werden soll. 

Ein Plattenkörper setzt sich aus 2r +1 
planparallelen Schichten zusammen. Die innerste 
Zentralschicht ist aus einem leichten und nicht zu 
festen Material (wir sagen weiter kurz: Füllung), 
die übrigen von verhältnismäßig schweren und 
festen Materialen. Ihre symmetrische Anordnung 
um die Füllung, sowie die gleiche (im Vergleich mit 
der Zentralschicht kleine) Dicke von je zwei sym- 
metrisch gelegten Schichten ist aus Bild 1 zu er- 
sehen. Weiter sind immer je zwei um die Füllung 
symmetrisch liegenden Schichten von gleichem 
Material vorausgesetzt. 

Die Füllung befindet sich zwischen den Ebenen 
2= + 2%, die #-te Schicht oben (die Numerierung 
erfolgt von der Zentralplatte aus) zwischen 2=2;_1, 
z2—=2; und die -te unten zwischen den Ebenen 
2=2_i41, 2=2_,. Dabei ist ganz allgemein 


u 22h NE 2er 


Bild 1 A a Pe 


(1), 


Ü 
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wodurch die z-Dimensionen des Körpers festgelegt sind. Das &-y-Profil ist ein Rechteck mit 
den Seiten a, b. Die auf die Füllung sich beziehenden Größen werden im folgenden ohne Index 
geschrieben, für die 2-te Schicht oberhalb wird der Index 7, für die i-te unterhalb -% benutzt, 
wobei alle Indizes rechts oben an die nämlichen Größe in Klammer beigefügt werden. 


Die Grundgleichungen des Problems (es wird nur der statische Fall betrachtet), nämlich 
die Gleichgewichtsbedingungen: 


Or DT, ER QISSn 3 dry. |. 00, GE ROTEN dr,, | 00% h 
== - ö — — - =: = . 2a 
0% 7 0% er 02 0% 0Y 7 02 0% | 0% r 02 Si En 
das Hookesche Gesetz: 
ran eh Kiss A et Ber Ms Mas Te. 
En Era Er EATR | 
2.2 
Tey Tyz Taa 
A el m: iZen, 
mil 
Hs Me Ms Me Han Ms (2.3) 
EB Ban Basis Ei : 
und die Gleichungen 
ou v ow ou, mw ov : dw ow gu 
=, „>75 a = n5 u = Fa) 2” "=. RE I me a: 2.4 
a  Yay RT N le’ Ur a 


vereinfachen sich grundsätzlich infolge zusätzlicher Voraussetzungen, worüber jetzt Näheres 
gesagt werden soll. 

Was zuerst die Füllung betrifft, so kommen hier gewöhnlich Materialen mit den Elasti- 
zitätsmoduln von der Ordnung etwa nur 2.10° kg/cm? in Betracht, während man bei den anderen 
Schichten mit Hunderttausenden von kg/cm? zu tun hat. Deshalb setzen wir für das Gebiet 
2, <2< 2, näherungsweise 


BE NE 0, 


Die übrigen elastischen Konstanten E,, @3; @s3, Hgı, Hg; und auch die Größen &,, &y, Yzy 
sollen, der Sache angemessen, endlich sein. 
Nach (2.2) und 2.3) ist nun ganz allgemein 


1 1 1 
& — 7 (0, — Ua Oy — His 0,), &, — (2 Haı Ix ze Gy —— Uaz 0,) ( Uz10x Usz2ly Fr 0,); 
E, BR E, 


so daß die Voraussetzung E, = E, = O und die verlangte Endlichkeit von &,, &, notwendig 


O4 — Ua Oy — Hz 92 = — Yaı 95T Oy — Hay 0, — 0 
zur Folge hat. Aus (2.3) folgt aber in unserem Falle weiter 143 = sg; = 0 und die eben auf- 
gestellten Gleichungen ergeben 0, =0,=0. Für g, folgt dann &, = 5 und endlich ist 
nach (2.2) [wegen G@,5 = 0 und der Endlichkeit von yzy] ty = 0. 
Zusammenfassend hat man also für 2,<z2<2: 
Bee nein een ,=0 |. .:....6(8) 


und die Grundgleichungen (2.1), (2.2), (2.4) gestalten sich für diese Zentralschicht wie folgt: 
On ya 20T, © 


Toy, 09% 9. 
Om Me. - ü N Te 
sl Betr wu far N 4 
rhnterg, mug ( (4) 
_Ww , _w 0 „ _dw, du 
rn ne 


Die Integration ist hier einfach, und die unbekannten Größen werden in der Reihenfolge 
Ta Tyz Oz, W, u, v berechnet. Die dabei auftretenden beliebigen Integrationsfunktionen (man 
hat ja mit partiellen Gleichungen zu tun) hängen nur von &, yab und werden 9, - . ., 9 ge- 
nannt. Führt man die Rechnung wirklich durch, so erscheint die allgemeine Lösung, unter 
Beachtung von früheren Ergebnissen, fürdie Füllung z2,<z<z, in der Form (statt E, 
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schreibt man, um später Mißverständnis zu vermeiden, lieber #,): 


op, MP 
A 2 ); Tey = 0, Tyz = Pa Tzx = Pı5 


oT 0y 
ME [9 en) RR. le a] 
ei. ed PR: 2 E, 9% 9 3% ai 0% 
a 9), 
en 4 (2 zen 2ER | N, (9) 
A ET he RN) Sl 
en ln 9pı )| 
Ver % Ir - A % 0 


Kleine Dicke und große Steifigkeit legt es nahe, alle anderen Schichten als Platten zu be- 
trachten. Man setzt hierfür 
EB = de = 00, red, = Esendlich, 


so daß auch 1 = Hy = u und 6, =G — Sn) endlich ausfallen. 
Unter diesen Umständen sind die beiden letzten der Gleichungen (2.3) nur dann möglich, 
Wenn il = U —= 0 ist. Daraus folgt aber nach der dritten Beziehung (2.2) weiter e, — 0, 
während zwei letzte Gleichungen (2.2) [wegen 63 = 63 = %]Yyz =Yza =0 ergeben, Zu- 
sammenfassend gilt also für jede die Füllung umgebende Schicht () 
E, = Bs—=:E:endlich, £, = Gy = 60 00, a el endlich, 

= Hg = WUnd ug; endlich Kıs = Uas = 0; & = Yyz = Yzu — 0 und die Grundgleichungen 

en) 1), (2.2), (2.4) gestalten sich folgendermaßen: 


RL er 00, |, ya __ Mn | Tzy 00 
— —— m z —e?. — 8 
0% r 0Y - 02 r 0Y r 02 02 ' 9y 02 
1 en Te 2(1-+ ! 
sen — un), genannte), ae ttn, (6). 
ou ER ow er ou, Mm v ow ow ou IR 


Bi ZA Duni 7 aa 7 a7 Tai 7 
Auch hier ist die Integration einfach und die gesuchten Größen werden in der Reihen- 


folge w, u, ®, 0%, Oys Teys Tags Tyz, 0, berechnet. Mit beliebigen Integrationsfunktionen 
®,,..., D,nur von x, y bekommt man 


oD o® 
v—=o®, u 02, ee mt 
ee dike: 08, 0b, 0° 0° | 
ae I (& +u|=; A let Ey een or. =) Bi 
E | o®d 1 2Ps ( fo E, | 
R — a |f a ET GIF = 0y? 2 G 


e E 
0y = 1-8 (HE; + „=; 
= : E © | >) E Ex 4 08 9% _ 
(1+u)\oy 9% 2(1+ u)\0y 0% 0% 0% 


er ae Zend RK Le = 
"| Ps (+ 04)” =. Te oy2 P, (7); 
I LP ) | 
S,: — 25,0% 2 EG 
5% OT, ee) 1 zE le 0? 2) 
<a P; IA Ti oY E m Te 2 002 | Ay? D; 
EL N] 
+5 2 Va AL 
> ne ee zu) |= (2 = 
.|-®, |& we 
. z2E ( 09,1 +2 
Fa) a ae 
! 92 92 
der Operator A ist durch A = Fr + dp definiert. 
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Zu der für die Füllung geltenden Lösung (5) tritt also hinzu noch die allgemeine Lösung 
für die <-te Schicht oben: 


E; [od oo») (2 0° | | ) 

“ln 2 RM I (i) 

9% ale a dt) Pi 
E; | oa AD. | 0° es 2 | 
(?) = 2 - 2 3 “ m (i) 
Ze 1 wi” 0% 627) re 0? 2 
DU) =] 2?E, ( BU 005W 2 ) 

W— DU — = 4 2 ea ei @). 
nr = re) WIE re 
ne (eh 

»v» 2(1-+) Ay: 02 0Y (8.1), 

zE, [f1— )® 0 .... 2 ) 

BI Di % () Se 2 ur (i) 

4: Ps al 2 n | aD 00. Wi 
zE, [[ & 1— u; 2 = ( +n08% z ) 
) — DW i i (i) ade 1. Fi a))l. 
Na Ser, (= N Day |; 
(i) (i) { > 
) a ni us vu) — DB) — z Be DR ud: 
} 
PETE EHIR, 1) 23 PR ARE (' 
sowie die allgemeine Lösung für die»-te Schicht unten 
E, [a8 oD(-") ( 0 0 ) 3 ] 1 
(—i) — ® a ; 3 TR (—i) 
Put | a aa th) Pi : 
BE. EX Ye Ya) ( 0° 02 ) et | 
Pe Re v ‚ 2 3 ra u et 
9 1— uw BR 0% 0 Eat 9y? 2 
e od en) z2E, I 98) 2 =) 
r it). I 5 | > A (—%) R 
P, 2 00 or STIEFEL Trek 
ae TA ( U RL RL ) 
vv» 2 dtm)\ 9 0% 0% 0y (8.2 
=) 
zE; “. RR: ) RUE 5 2 >) 
1, Dr 12 2 002 ' 92 ®, Piayx 2 0% 2 a7: 
er, —i :E |(® | 1m # )® ie tmoB") 2 o) 
m —®. + 2 02 BrReR . fra 
r od) ; ; „DE 
N N Ne u Dam 6 =) — B-; 
ı) U = D\ 2 I r V® &D\ dy 1 


ER w—=1,2,...,n 


Der 


Nach der anfangs erwähnten Bezeichnungsweise sollte eigentlich überall in den Formeln (8.1) 
E®, u® und in den Beziehungen (8.2) E®, u —®stehen. Jedoch setzen wir schon vom Anfang 
voraus, daß die ö-te Schicht oben vom gleichen Material besteht wie die ?-te unten; darum ist 
EÜ — E, u) — u”, wofür man einfach E;, u; geschrieben hat. 


Die Lösung n% (8. Ri 4 2) enthält a + 1) unbekannte Funktionen 9, ..., Ps; 
DU,..., BO; DB, ,.., d@—=1,2,...,n) von @, y. Sie sind im Einklang mit den für 
die Ehenen == 2, BEE — ne {es =0,1 2... n und für die Seitenflächen = (0, 2 —=«, y=(, 
y=b geltenden Bedingungen zu bestimmen. Jetzt wollen wir uns mit den erstgenannten 
eingehender befassen. 


Unter Voraussetzung, die Platte sei an ihrer oberen Stirnfläche z2—=z, nur durch den 
Normaldruck q,(z, y), an der unteren durch g,(z, y) belastet, hat man für 2—= 2,: 


mm Mn dy) 
für 3 2, (= — 2): 


er ne). 
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Nach (8.1) und (8.2) drücken sich diese Ban Be aus wie folgt: 
2 Ben 1 (n 
DW = En ( 0 — 1 Un 0 )® (n) ur ai Er Un od" -7 400) = 0 


— un L\0x°? 27.07 2 m 
ee le u ee on: 
FRE ar (= Eu : P a. = 1 Pr er ( _ u. % 3 20) u. 
er) an (er Zon)| =0 02). 


An den Berührungsflächen zwischen je zwei benachbarten Schichten wird natürlich 
stetige Änderung vod Spannungen und Verschiebungen verlangt. Mathematisch heißt das 


fürs a al, 2 le 


Turin, Wert), Mod; en ud, ae — url 
und für @=r (= —- 2) Lt = la nn): 
(—) — 2-1) = (— i—U), (—) gel): 
ee 27 ? T, Ir 0, 0, 3 
ya = ui» ; ve ee vi), va) = wii), 


Unter Benutzung von (8.1) und (8.2) bekommt man daraus erstens 
z; E; ( ® 1— =) ak; ar 2; ) 
Di ee pe (ü) 4 377 TA DU 
®, 1— uw et 2oy: NE 02 Ve 20% 2 2 
1-+ wir oDü+n 


- a Dr I 97 2; “ 
— Büi+N) 3 le 5 Ir DÜ 5 aV >49; Hl) 
,E |[1—- u % 0? ,)® in od | 
@)__ ı ı ae DU) RR Tas . Ua ee 2 (i) 
®, 1— u ( 2 5000 z 07 2 702 2 AP; 
; 24 Bi+ı ( — Heel a Ya 1+w40BdfH) * 
ee er 2% EEE ae 
" — Wir 2 00a) ® oy\ - 2 0% (10.1) 
BD) &i E; ES oa ) 
End | () 
5 =; ( 0% 0y/ 21 21 — A) 0% 0 3 au 
KR Lin Sy 2 By, Benz Bi+D Ar ) 
— pH) 2 DIE 
D, ( 0% aa — uf4ı) 0% 0y 3 u. 
od“) opü+n ’ od) oDü+n 
KaE I 1 DEMENZ, 1 KaRat ae 1: — Hut): 1 } 
| % &\ pr D\ 2; dy D\ 2; day 
DÜ — Dü+); En Fe Te 
und analog 
BE (® 1— uw ei 1+ u, 007 2; >) a 
(—i) we, = (—i) % —))| — Bi) 
en, = wi oa+ 2. 99% er 0% Al: 2n 00% + 2 By P, 
2 Eizı 0 ln 2) ee ( tm 0DTT | 2 | 
au 1— us (+ 2 ay? ®; r 0% 2 oy ii 2 rt: 
;E |[1—- uw 9% 0° nn (- +m9dT) 2% ) Br 
(-i) Bed u (—i) Zi AB-)| — B-i—ı) 
er, ( tz) A wa ne mat ®: Ba 
% Bi+i (@ — BI, =) ® en re Hi u % a) rt 
u — ur 2 dat a)? Ar 2 0% ar 2 49; 


O 
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} od HD a EB eo aaa 
Di = 4 5 ) ii | 2 3 Er >) ii But 
R 0% 0Y = 2(1 — ui) 0% 22 oy Ar 3 or Pen 
od -i—) ad-i) 22 BE: adl-i—V Di) 2; 
+ 2; ( Ei 5 { i+1 | 2 3 % ne 
0% | oy 2(1 — ui4ı) 0% [67] 3 ag: 
Fi aD) od ie a Dar) 
od i) 2; — BA) ; eg j DB) 2: = Bl-i—1) : 1 
in 0% s > 0% a oy Bor 7% Day K: 
DI) — BI, 3 RT 


Dieselben Bedingungen für die Stetigkeit von Spannungen und Verschiebungen gelten 
freilich auch an den Trennflächen zwischen der Füllung und den beiden sie umgebenden 
Schichten. Nach (5), (8.1) und (8.2) bekommt man folgende Bedingungsgleichungen: 


Beam 1m 0 Eon; 
a Te Ar Hı 1) Hı OD; zo 
ERIK le Fa ir) Ya | 2. 3 = 100) 


E, I1— o/l+modV 7 
— BU _ 2% #ı |( Hı © )® je | U 2 %o »]) 
a 1—ı2l\ 2 a, 2 Ara 2m dr at 
_, [9 2) — a „(er > 2 Zn Br D» dd 5 ») 
% (de Rı 9% De 0% 04 1 12) 0% 0 ee 
Ialal))s 
ee, ala Blwiop,0P Sr 
P5 0 — o| 
9 2 E, 0% 0% 
2 ni 2 | E23 we | - a) 98; 
Re Io: 04) 2m,” 3 eu are) 
Kl er = — GW 
+2 R 9) Per =6d; 
2,E, [[ % 1— u, 9‘ 1+m008V9 2% 
u per (5 iı ee) a Hi 9 ev ») 
m ln ah? I as Kae) 
E f — tt 02 =) 8 ud 2 ) 
EINEN 21 ee AR Hı 0#a 0 = 
an ae Ft 
op, on) — Er Be N 
RA: dy B"+% 9% 4 dy 
22 E, Se 002.2 >) 
B 2(1— u) ey 163 he (11.2). 
2b = a = a (M e)| — Te 
2 a (E 02) 28,0% ont oy/ı er 
als er) 2 ae | 0,0 od” 
Pe n(d8 oy/ 28,04 er N 0 ee ay 
FR. 20 (9 | — 8-0 
Pı 2 |m+ 5 (2 a 
Den aufgestellten Bedingungen (9.1)—(11.2) entnimmt man leicht, daß sich alle die 
6(2% + 1)unbekannten Funktionen 9,,...,965 Di, ..., Did; DI=Ü,... o-) (wel) 2,2...) 


durch 8, 85”, DI), &,"), 85", DI” ausdrücken lassen. Das ganze len reduziert sich 
also auf die Bestimmung der eben geschriebenen.6 Funktionen, oder — was dasselbe ist — auf 
die Bestimmung der (übrigens noch nicht voneinander unabhängigen) Größen 


4, , am, A”, 2", B2 i Ba) BE EC, (12.1) 
wenn man ganz allgemein (also auch für die anderen Indizes) 
er BD=z0) - B"i=l2..,n,k=12,...6 (12.2) 


setzt. 
Aus den Beziehungen (9.1) und (9.2) folgert man dann mit den Abkürzungen 
E, I 1 | 
= £: a » Bea) =]; 2, 18.06% 13 
31m) &s 5 Ks 5 Bi; ? n (15) 
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leicht 
er | 
BT Ar Esel 


und die Bedingungen (10.1), (10.2) ergeben 


5 0° 0 oBÜ | 
fü) _ 9% 5 Pe le dr BETTER i) 
As 2 (04 Pi a2 a) > 1 le % 2 Fe: ) 


i oBli+n 2; a 
— Ali) _ 2847 Ale: A PiH15 a) BED + lan —& 140%) 


0y 2 
2 ;r R 
a a a en) 
A 


(i) (i 14) (i) 
a@ „(Bi “) aaltd "2a Be 


Br) r 
a 


(+1) ee [3 (i-+1) 
Ba ur 0% a ) (151): 


0% [oR7] 0% oy 
‚(eBi> Bü+n i+1) AU+D F 
e Bü) j (#1) (i) i Bü+n 
AND — 2, = an ne AD — 2; — eure a R 
AD — Ali+n, ER RL 


ei: 0° n) rd ( 0AU a | 
le je WE = | ae SERFREREE (i) 
B\ 2,2% (+8 ap AD+ 55 \ dy 5 AB! 


i 2 0° 9 0) Ay 2 
= BÜ+V — 28412 (= Bat 5) a va FE 3 bar yore 


| Be+) 
% 

2 2 
BO 2u (By 2) aw+ 2 (m —& 42%) 


0Y 2 


0Y° 0% 
— Bit) 941% hei: an nn) an ln — & aBeHn)) Ah, 
BO— 2 ZN —_ Bi+V nn BO — 4, 2 = Bin Rate) 


Bü — Bit); era, 1 


) 


| 
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EEE 


Endlich hat man nach (11.1) und (11.2) 


PEN EI APR Bew, 

A at - © Er | — 5 an). 
en Hesshifogen 

2 (et + er le | Fe a 1an) (16.2). 
al Den 
tet ml 


Die Gleichungen (14.1)— (16.2) ermöglichen jetzt alle in Betracht kommenden Größen, 
nämlich 


AP, BP d=12...n k=12..,6); 9, 90::.,08 
durch A, Aw, A, BW), Bi), BW) auszudrücken. Damit wollen wir uns nun befassen. 
Zuerst ergeben die letzten der Beziehungen (15.1) und (15.2) 
N RE TE ET 


und daraus folgert man nach den zwei vorletzten der Gleichungen (15.1) und (15.2) sofort 


Bean Be Bu = Be el, 2... 
Im ganzen hat man also 
BAT Ep Br Ber. =l2,3) ©: 0.04 (171) 


und es bleiben noch die AW, BO6—=1,2,..,n;k=4, 5, 6), 9, Pa - - -» Ps zu berechnen. 


Für Aw, Aw, B@, B{» ist das schon durch die Formeln (14.1) und (14.2) getan. Um 
auch A und BW) zu finden, hat man nach (14.1) und (14.2) 


Am 9A (Bar oBm z 
4 ER; a TEEN (n) 
dw dy TA TTIREE Ta, en 
7.2). 
(n) (n) (n) (n) 
oB! oBi _ REN er 04! An Ay) 
0% ay 0% 0 2 
Damit erhält man aus den letzten Gleichungen (14.1) und (14.2) leicht 
(n) (n) 
An a A) ++] 
| Br 2 > . (17.3). 


oBm oBw 


2 
Bm + u rap) + | 
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Weiter werden wir uns der Bestimmung von AN Bi e@=12, REIS) 
zu, was mit Hilfe von (15.1) und (15.2) gelingt. Danach ist mit Rücksicht auf (17.1) zuerst 


92 02 ) 
Am = AMT) 22 ß - ir) th & — Pi+1&i+ı) En Bi) 
Ir 0: Bw ß) 
+2 2,(0; & — &i41 &i+1) ar 2y zi(e &i41) 5, AAr 
- > es (18.1) 
A — PK) 1) 2 2; Be mei ßi 11 &i+1) dm? | (e; &i+1)7 ko | B 
+22, — K+ıe sr il — Ei )3 „a4 ee 1 
il & — Mir &itı) Sn, le i+1) = wa bar un 
Hı &itı) 7 dy +1 
92 
er} H+n+ 2 alle 4) d% Be a ee ar Pi +1 2a 


02 Alm) ö 
+ 223(0; 8 — +1 &41) EEE CB &41)5,ABr" 
(18.2). 


0° a 
Be = BET) + 2%, R & — Piz &i+ı) at (& — &+1) a A) 
92 Am 


2 — Hl; — gr ABW; a —A; ER 


+22(% & — &i+ı ag 


Die dritten der Gleichungen (15.1) und (15.2) ergeben unter Benutzung von (18.1) und 
(18.2) nach kleiner Rechnung 


An) Am) 2 2; , 
AD = AED + a — au) Al + A) | 


(18.3), 


5 . oB\ n) 0B") 2 2; j . 
BO — Bit) + Z( — &+ı) 1 =: Sr TE: aan); a | 
womit alles zur wirklichen Auffindung von au, Br G=123,;,n—1; k=4,9,,6)vor- 
bereitet ist. 
Mit den Abkürzungen 


n—1 n—l n—1 


& (& — &+ı1) 2 = Sm » SM 8 — &41&4)4a = 6: &— Brr&H)s—=Z (19) 


bekommt man durch is über ©=1,1-+1,...n— 1 aus (18.1)—(18.3) leicht 


02 Bin) d 
AU — Aw+ 2(, 0 1 a“ I; =) Doaa2 Yıes'r TR 87,44 
ANE Am 214 Es Sı, 1 Bm + 2Y, Tu By EB S7 ea AA") 
oy? 0% oy a 
04\ n) = 9) 
ee) 3 Ba (n).» > u 
AD — At Al ta) nA; I=1%.. nm 
BO — Bm 2(s,,2 tz =) Aw + RAR N: 25, ZB 
4 4 wu oy2 x x 0y 
BU — Bm + 2(Z Am Hay _ 9,212 
= Eh 1 + en Ey 1a 


oBm oBW 
0% 0Y 
Führt man hierin nach (14.1), (14.2) und (17. u die Aw, Bw; k—4, 5, 6 ein, so hat man die 
gewünschte Darstellung von AawB % G=1l 43.210 ? k = 4,5, 6) durch die Am. A, 
Aw, Bm, Bm, Bm 
02 Bm 


0° 
Aau=2 + En 2m) 2; da? + (Zi + Pn En en) 5 a BRACH) 02 0y 


BD — Bin 5, 4( \-2 8,444 1=12,...n-1. 


(20,1) 
—e (1,2 1 En 2.) AA 
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02 B® 


0% 0Y (20.1) 
Mrz (81, 2 fr En %n) >, „dan 


02 0° 
Anz 2 1 + Pr En 2) 92? a (S1,1 un En Fe op? 1 B®+2 2 (Y, + On En m) —— 


0A Am 2 
ta) ta AAN ++ Mm 


Deel,2,..3n—1 


AP—=(S,a +) A | 


0? Am) 
0% 0y 


ER (S1,2 5 2) 3 => = uBi 


0° je 
Bi —2 [& 8% Zn) I2 = (Z, Fr Pr En DE a am +2 (Yı FIR OnEn N 


02 Am) 


Bo — 2l(z 90; 
end! 3e29 (202). 


En En) Sr (5 Et a En DER | Am Pa @4 ier An En I 
Pe (S1,2 Er En ) AB 


0 Bi) oB" 
0% 0% 


2 
) 3 Sat m 2) AHA + N— 5 
Mr 2: ‚n—|1 


Zum Schluß haben wir noch die g,, ge durch die Größen (12. 1) auszudrücken, wozu 
wir uns der Beziehungen (16.1), (16.2), (17. 1) und (20.1) bedienen. Zuerst hat man nach 
(16.1), (17.1) und (20. Er 


92 
n=|(8, 1 &2%+ En u) + (Zı — Ba Zo + Ba En Zn) ER Bi) 


2m 1 f 
9% dad 3 Sun EN. AA 


BD— (S1.+ 6,24) A( 


) 


+ (FL — 98204 On En Zn) 
22 22 (21.1). 
Ps, @ — Pı&%+ PnEn?n) Fe) 2 A EN DE sr By” 
02. Bi" ie 
0% °y Da! 

Daraus erhält man weiter nach schöner Rechnung — insbesondere beachte man, daß im 
Sinne von (19) 


Ir En RAR — dt en) 2, AA 


z Y, +4=%9ı1=12...,n—1 
ist — 


og, , 2 OB oe 1 2 
Dr Br Alan +) Sata] (2) 


Die Berechnung von 93, 94, 95, Ps gestaltet sich nun ganz einfach. Nach der letzten 
Gleichung (16.2) und nach (17.1) ist 


E 
= > no... (21,2) 
und die drei letzten Beziehungen (16.1) ergeben dann 


(N) (n) 
DE EB SCHE LER EN EEE DE ICE Zr) 
X 3 | (21.3) 


ei a -,( Be 
Ei ae 5 5 2, 


Zusammenfassend drücken also die Formeln (14.1), (14.2), (17.1) (17.3), (20.1), (20.2), 
(21.1), (21.2), (21.3) alle unser Problem charakterisierenden Größen 


IR By (i ne 2 6), Yası- - -» Pe 
durch Aw, Bw (k—1, 2, 3) aus. Unser letztes Ziel ist die Aufstellung der Differential- 


‚gleichungen für diese 6 Grundsteine der ganzen Aufgabe, sowie ein paar Worte über die zu- 
gehörigen Randbedingungen. 
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Den Gleichungen (16.2) und (20.2) entnimmt man unter Berücksichtigung der Be- 
ziehungen (16.1), (17. 2 leicht 


6 18:20 + En 2) 5 a (Zı — Pre + Br En en) 5 a Ay 


02427; Al 
+, —- n.%+ %n En Em) 5 2 0y 5 (S,. — &% ni) ABM —.0 
RB — Pı&ıo + Pn En ap Ep vn (81,1 — &120 + En %n) | Hi (23). 
z 02 A d 
+, —- 9n:%+ ann u) ad 5 Su2 — 42%3+ na) AB” =.) 


9 Bm) Bm 9 
da aY ) 3 81,5 = &20 + &n2n) AAAP”+H—-p—=0 


Vergleicht man den aus (16.1), (16.2) und (21.2) berechneten Wert von A{® mit dem, 
welcher sich nach (20.1) ergibt, so kommt man zur Gleichung 


94) >42" 2, 
0% 0Y 3 


Differenziert man nun die erste der Beziehungen (23) nach x, die zweite nen y und addiert, 
so folgt 


(Sa. tm) A | 


(S1,2 — 04 5%) A | S1,3 294%) AAB” — > se AtR—). 


0A) 0A) 1 
(Sı,1 tn %) A ( er 3 ) — D7 (S1,2 + SRH 27) AAB”. 


Zusammenfassend erhält man so die Gleichung 


481,1 zart) (S1,3 — +4) — 3(81,2 — 2 En 
6(8,],1 — &20+ En 2m) 


m)” N E; n 
AAB® + MP —at% 


(24). 
Endlich folgt nach (16.2) mit ee auf (17.1), (21.1), (21.3), (22) und (23) 
2 
Be 2 a (81, mat 2o+ En ae = (Zi — —D, & 20+ Bu En Zn) = B & 
| 0 a 
+(Yı — &E120+ An En in) —— 2 ET zZ (Sa — at 52%) 3, di ) 

22% f) a 
— & 2 + En) (Sı3 — at m 
 3(8,8— 123-4 &n 22) B, 0% “4 (5 3.17 o+ IR 1,3 Bi) ) 

— 3(8,8 — 81254 © n)?] AAAI® — (81,1 — & 204 En 2m) (Qı — 1) 
(25)® 


- 


Bo RZ, — Arzt Ban), +8 — 2%+E N Be) 
B Gl a) n nn dn2 1,1 10 nn) dy2 3 
u 1 
+(Z, ee Eee $ 


>: a 
2 — & 204 En 2m) (Sı.3 — E12 + &u 28 
3(81,2—&%+ en) B, öyl6 [#191,1 — erto + ent) (Di ) 


— 3(81,8 — &1264 &n n)?] AAAS” — (81,1 — &120+ En 2m) (dı — 2) 
Die Beziehungen (23)—(25) stellen das gesuchte Differentialgleichungssystem für die 


unbekannten Funktionen A”, Bw; k—1, 2, 3 dar. Durch Einführung neuer unbekannter 
Funktionen 


0) 
— (81,3 ie 25 LH & ) dy aan) 


(n) (n) in) 
mA 2 = 1; 2 Pr 


Saat 

2 w—= Bm, 2) = Am) u) n en 

a T ns » 2(8,1— &20 + &n %n) 
und der Ausdrücke 


2 =. tg 2m = (26.2) 
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zerfällt das obgenannte Gleichungssystem in zwei Gruppen: die beiden ersten der Beziehun- 
gen (23) und die aaa (24) ergeben 


92 
C# — gt %%) e + (4 — Pı&ı2o + PrEn 2n) 23 N 


is 
+(Y, — mt Er 
, 0° 
Zi ham + &1% HF EnZn) & (27.1), 
0° 
ahkuns Bee DE 


4(S1,1— E120 + En m) (S1,3 — &1 25 + En 23) — 3(S1,8 — 123 + &n 22)? 
’ ’ a A RL FE EEE ER I er Da In ER Be 1 Ed ME AAw u — 
6(S,,1 — E20 + En) 1 SEEch 
die Be (25) und die letzte der Gleichungen (23) gehen über in 


29 t+ m Ow, 
se 1 &1%+ En?n) 9% 


2 
== all, 1 &2% 4 &n a em (Ze Pi &ı20 + Pr En ar 


)) G13 


ne+ 


& 
+(Yı — NE > 


22% ae 
— —i—-[4(S,ı — 3% 4 u) (Iı.3 — 122 + ©, 28 
aa le amt nm) A at nm) 


3, + m] Alm — (a ante) Pi) 
mg zo tn 2a Ows 
a 27.2) 
22 zZ 
Se I@ —Ppı &ı ?0 in Pn En a) + (8, ı 7, &% a En ar 2|e 


+ 


Y 2 
+( 1 E20 4 An En u 92 0y 


22 an 
sr — 7. _—_1—[4(S —— BR 82% Ss — 22-4 .,28 
euere) at en) 


— 3(5,8 — &2 + En A)?] AAwg — (Sı,1 — Ei2o + En m) pr 
4(S1,1— &1%9 + En) (8, & 29 + En a) — 3(S1,2 204 En we Arts 


@) 6(S,, I got) 
On 2) 
E De ee 2 2 ZA SD 
(S1,2 at en) [et Toay > 2° ) 


Im einfachen Spezialfalle a = 1 gehen unsere allgemeinen Gleichungen (27.1), (27.2) 
natürlich in die Prusakovschen Beziehungen (1.34), (1.35) über — s. seinen am Anfang 
erwähnten Aufsatz. Der Leser möge sich davon selbst überzeugen, nur muß er dabei beachten, 
daß die Größen S,, 7» Y,, Zı für n= 1 als ‚‚leere‘‘ Summen gleich Null zu setzen sind. 

Die Randbedingungen hängen freilich von der Art ab, wie die Platte an ihren Rändern 
befestigt ist. Näheres darüber findet der Leser auch in der Prusako vschen Arbeit und es 
wird ihm möglich sein die dort aufgeführten Schlüsse und Erwägungen sinngemäß auf unseren 
allgemeineren Fall zu übertragen. 


Eingegangen am 13. Juni 1952. 
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Zusammenfassender Bericht. 


Fortschritte der Turbulenzforschung”) 
Von W. Tollmien in Göttingen 


Ein zusammenfassender Bericht über den Stand der Erforschung des Turbulenzproblems 
ist seit einer Reihe von Jahren nicht gegeben worden, was bei der zentralen Bedeutung dieses 
Problems für die Mechanik und ihre Nachbarwissenschaften einigermaßen verwunderlich ist. 
Auch die einschlägigen Berichte auf den Internationalen Kongressen für Angewandte Mechanik, 
die 1946 in Paris und 1948 in London veranstaltet wurden, bemühen sich kaum darum, die ge- 
samte Problemlage zu umreißen, sondern behandeln mehr oder minder Ausschnitte, wie sie der 
Arbeitsrichtung der Vortragenden entsprechen. Da diese Berichte zudem in manchen Ländern 
schwer zugänglich sind, habe ich mich entschlossen, der Anregung einiger wissenschaftlicher 
Freunde zu folgen, auf dieser GaMM-Tagung über die Ergebnisse und die neueren Tendenzen der 
heutigen Turbulenzforschung einen Überblick zu geben, soweit es eben im Rahmen eines einzigen 
Vortrages möglich ist. Ich wähle dabei als Ausgangspunkt etwa das Jahr 1945, da seit diesem 
Jahr wieder ungehindert Grundlagenforschung möglich wurde und allmählich auch die während 
des Krieges gewonnenen Ergebnisse bekanntgegeben wurden. Bei den weiten Auswirkungen 
dieses Gebietes sowie dem überaus raschen Fluß der Spekulationen daselbst ist freilich wie- 
derum eine persönliche Färbung auch dieses Berichtes unvermeidlich, wenn das Übersichtsbild 
nicht zu verschwommen sein soll. 


1. Entstehung der Turbulenz und Stabilität der Laminarströmungen 
Am befriedigendsten ist die Lage, soweit es sich um die alte Frage nach der Entstehung der 
Turbulenz infolge einer Instabilität der Laminarströmungen handelt. Zur theoretischen Unter- 
suchung werden dabei einer laminaren Geschwindig- 
a keitsverteilung U(y) fortschreitend Wellen überlagert, 
IS Ua IRR die in Partialschwingungen aufgelöst gedacht werden, 
Den, von denen eine durch die Stromfunktion 


BAR ply) nee 


dargestellt ist, wobei i die Zeit bedeutet. Die unge- 
ER störte Strömung U(y) soll im wesentlichen parallel 
zur Wand, deren Richtung in die x-Achse falle, ver- 
laufen und merklich nur von der Koordinate ysenk- 
recht zur Wand abhängen (siehe Bild 1). Grenzschicht- 
profile sind die praktisch wichtigsten, welche von diesen 
Stabilitätsuntersuchungen erfaßt werden. a ist bis 
auf den Faktor 2x die reziproke Wellenlänge der 
Bild1. Grenzschichtprofile an einem umströmten Partikularschwingung, der Realteil von c gibt die 
ie Rue LEERE, der vr gegen die Phasengeschwindigkeit, während der Imaginärteil von 
ac je nach dem Vorzeichen ein Maß für die Anfachung 
oder Dämpfung der Schwingung ist. Der Realteil von &c werde mit ß, bezeichnet, er ist die 
Kreisfrequenz der Partialschwingung. Die Differentialgleichung für die Amplitudenverteilung 
p(y) der Störung heißt 


rt 
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Rist eine Reynoldssche Zahl, welche aus einer charakteristischen Geschwindigkeit (etwa 
U,) und einer charakteristischen Breite (etwa der Verdrängungsdicke ö*) des auf Stabilität 
untersuchten Geschwindigkeitsprofiles U(y) gebildet ist. Es wurde nun eine „asympto- 
tische Stabilitätstheorie‘ entwickelt, d.h. die obige Differentialgleichung wurde 
bei sehr großen Werten des Parameters R unter Beachtung der Randbedingungen integriert. 
Diese asymptotische Stabilitätstheorie führte nach fördernden ArbeitenvonLordRayleigh, 
L.Prandt1l,0.Tietjens, W.HeisenbergundF.Noetherinder von W.Toll- 
mien und H. Schlichting gegebenen Form zu einem vollen Erfolg (für die historische 
Entwicklung dieser Forschungsrichtung vgl. [1] und [2] des Literaturverzeichnisses). Einmal 
wurden die theoretischen Ergebnisse von Tollmienund Schlichting für die längs an- 
geströmte Platte von G.B.SchubauerundH.K.Skramstad[3], die unter der Leitung 
von H.L.Drydenim National Bureau of Standards in Washington arbeiteten, in glänzender 


*) Vortrag auf der GaMM-Tagung in Freiburg i. Br. am 29. März 1951. Nur wenige Zusätze, welche 
das während des Jahres 1951 noch erschienene Schrifttum zum Teil berücksichtigen, wurden bei der Vor- 
bereitung zur Drucklegung dieses Vortrages in den Text neu aufgenommen, 
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Weise bestätigt. Bild2 zeigt das Wachstum gewisser natürlicher Störungen mit größer werdendem 
Abstand von der Vorderkante der Platte. Diese angefachten Störungen sind anfangs sehr regel- 
mäßig, und das merkwürdige Ergebnis der Theorie, das mehr als ein Jahrzehnt angezweifelt 
wurde, daß nämlich Störungen ganz bestimmter Frequenz und Wellenlänge angefacht werden, 
wurde durch diese Oszillogramme vom _. Pr 

August 1940 zum erstenmal bestätigt. rl a ' iR 
Noch deutlicher wird dies im folgenden 422 ______. 
Bilde 3, wo die theoretische Kurve, 
welche nach meinen Rechnungen vom 737 
Jahre 1929 die instabilisierenden Störun- 
gen umschließt, mit Experimenten aus 153 
dem Jahre 1941 verglichen werden. Nach 
den obwaltenden Versuchsbedingungen 168 
konnte nur der obere Zweig dieser Kurve 
gemessen werden; dieser wurde aber 
durch die Beobachtungen mit einer weit 


größeren Genauigkeit wiedergegeben, als reise a ER } 
ich sie je bei diesen Stabilitätsunter- 198 Y N www j 314 10 


suchungen erhofft hatte. Durch Einfüh- 
rung künstlicher Störungen konnten 2/ | rl da el Br: 338-10° 


Schubauer und Skramstaddie 
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ganze Indifferenzkurve, welche die an- 229 en SR .© 58 ging 362 10° 

gefachten von den gedämpften Schwin- Rare res j 

gungen scheidet, der Messung zugänglich 

machen. Auch hier war die Überein- U, = 24,4 m 

stimmung zwischen Theorie und Experi- 220 2. Oszillogramm natürlicher Längsschwankungen der Geschwindig- 
eit in einer Strömung längs einer ebenen Platte nach Messungen von 


ment überraschend gut, wie aus Bild 4 c.k.Se h ubauerundH.K.Skramstad[3]im Abstand 0,58 mm 
zu ersehen ist. Das folgende Bild5 zeigt vonder Wand. Us aa he NE EINER. x Abstand 
eine weitere frappante Bestätigung der nee Saar 
Theorie, diesmal der Amplitudenvertei- 
lung p zweier neutraler oder indifferenter Störungen, die von Schlichting berechnet wur- 
den, während bisher nur von den Frequenzen der Störungen die Rede war. Auch bei der 
Amplitudenverteilung ist die Übereinstimmung zwischen Rechnung und Experiment ganz vor- 
züglich. Die Bekannt- 
gabe dieser bewunderns- 
würdigen Experimente 
erfolgte durch H. L. 
Dryden [4] auf dem 
Pariser Internationalen 
Kongreß für Ange- 
wandte Mechanik 1946; 
ausführliche Veröffent- 
lichungen von Schu- 
bauerundSkram- 
stad kamen im Jahre 
1947 heraus [3]. 

Fast gleichzeitig 


mit der experimentellen 00 700 800 1200 1600 2000 2400 2800 3200 3600, 

Sicherung der asympto- —n RR. Und 

tischen Stabilitätstheo- Bild 3. Meßpunkte aus natürlichen Schwankungen für den oberen Zweig der Indifferenzkurve 

rie gelang ihre endgül- der Plattenströmung nach G.B. Schubauer und H.K. Skramstad [3] und 
theoretische Kurve, 


tige theoretische Fun- 
dierung, wodurch schwerwiegende Einwände von F. Noether entscheidend widerlegt wurden. 
Im Jahr 1944 konnte ich nämlich die asymptotische Fehlerabschätzung für die benutzten 
Näherungslösungen der Störungsdifferentialgleichung erbringen. Der Ansatz für diese Nähe- 
rungslösungen ist verhältnismäßig einfach zu gewinnen, dagegen ist ihre Rechtfertigung ein 
tiefliegendes Problem. Diese Arbeit wurde 1947 veröffentlicht [5]. 
Bei dieser Sachlage gehört die lange angezweifelte asymptotische Stabilitätstheorie der 
laminaren Strömungen heute zu den am besten gesicherten Teilen der modernen Strömungslehre. 
Die zahlreichen erfolgreichen Anwendungen der asymptotischen Stabilitätstheorie auf 
laminare Grenzschichten mit Druckgradienten, der ja bei der längs angeströmten Platte ver- 
schwindet, kann ich hier nur nachdrücklich erwähnen. Hierzu liegen besonders Arbeiten von 
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H. Schlicehtingund J. Pretsch vor, wozu man als Literaturquelle etwa [2] vergleichen 
möge. Die Tendenz geht heute dahin, Sondereinflüsse wie die Kompressibilität und die Schich- 
tung zu berücksichtigen. L.Leesund C.C. Lin [6] haben in zwei Arbeiten die asymptotische 
Stabilitätstheorie auf kompressible Flüssigkeiten zu übertragen gesucht. An diese Arbeiten haben 
sich ausgedehnte Erörterungen angeschlossen, welche hauptsächlich den stabilisierenden Einfluß 
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0 400 800 1200 1600 2000 300 2800 3200 U%3 
—— R= Eu ——- 
Bild 4. Messung der gesamten Indifferenzkurve der Plattenströmung 
mittels künstlicher Störungen nach G. B. Schubauer und H.K. 
Skramstad [3] und theoretische Kurve. 


Bild5. Amplitude Ver der Längsschwipgungen in 
einer Plattenströmung für zwei neutrale Schwingungen 
I und II, diein Bild 3 gekennzeichnet sind. Ausge- 
zogene Kurven nach Rechnungen von H. Schlich- 


einer gekühlten Wand betreffen. Ich möchte an dieser tin, gestrichelte nach Messungen von G.B.Schu- 
R N E aaa; = d t bauerund H.K.Skramstad. Beim Verschwin- 
Stelle nicht in die noch laufende useinandersetzung den der Amplitude ist das Vorzeichen von Y[p/]® ge- 
[7] eingreifen, die wertvolle Aufschlüsse verspricht. ee un a 2: 
Eine Schichtung, wie sie oben erwähnt wurde, kann experimentell an den Oszillogrammen sehr deutlich 
; . uf; ; ' festzustellen war. Die Grenzschichtdicke 8 ist gleiel 
durch eine ungleichmäßigeDichteverteilung,d.h.durch 37535 5* genommen, wo ö* die Verdrängungsdicke ist, 
eine ungleichmäßige Verteilung der Schwerkraft für 
die Masseneinheit oder auch die ungleichmäßige Verteilung der Zentrifugalkraft hervorgerufen 
werden. Der zweite Schichtungseinfluß tritt in stärker gekrümmten Strömungen auf. 

Neuerdings hat W. Wasow [8] die asymptotische Integration der Störungsdifferential- 
gleichung auf einem anderen als dem von mir in [5] eingeschlagenen Wege in sehr eleganter Weise 
diskutiert. Seine wertvollen Ergebnisse beziehen sich nicht ausschließlich auf neutrale Schwingun- 
gen, wie sie von mir seinerzeit in [5] betrachtet wurden. 

In einer bislang unveröffentlichten Arbeit hat D. Grohne das Stabilitätsproblem des 
„asymptotischen Absaugeprofils“ U=1—-e mit Mitteln aus der Theorie der Laplace- 
Transformation, die auf diesen Spezialfall anwendbar ist, exakt gelöst und asymptotisch dis- 
kutiert. Die Übereinstimmung mit der allgemeinen asymptotischen Stabilitätstheorie ist voll- 


ständig. 


2. Weiteres über die Eigensehwingung laminarer Strömungen. 


Auf ein andersartiges Problem hat kürzlich C.C. Lin [9] die asymptotische Stabilitäts- 
theorie angewandt. In einer gemeinsam mit J.R. Foote veröffentlichten Arbeit untersucht 
Lindie Instabilität zonaler Winde über der rotierenden Erde, wie sie sich etwa im Mittel entlang 
den Breitenkreisen ausbilden. Diese Winde werden nur von der geographischen Breite abhängig 
angenommen, die Störungswellen der Stabilitätsuntersuchung schreiten längs der Breitenkreise 
fort. Abhängigkeiten und Störungen in der vertikalen Richtung werden nicht berücksichtigt. 
Das so skizzierte Problem in einer Luftschicht auf einer rotierenden Kugel wird durch eine 
Transformation auf ein ebenes Problem abgebildet. Die Bedingungen für die neutralen und an- 
gefachten Schwingungen werden von Lin diskutiert, indem Methoden, wie sie von mir bei der 
Diskussion der reibungslosen Instabilität von Geschwindigkeitsprofilen mit Wendepunkten [10] 
entwickelt wurden, weiter ausgebaut werden. Es kommt bei der Existenz neutraler und an- 
gefachter Schwingungen auf den Vorzeichenwechsel der absoluten vertikalen Wirbelkomponente 
der zonalen Grundströmung an. (Dieser Gradient ist hier nicht mehr proportional der 2. Ab- 
leitung der Geschwindigkeitsverteilung.) Durch die Störungswellen wird bei westlichen zonalen 
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Winden ein Impulsmoment längs der Meridiane nach Norden transportiert. Die Störungen 
prägen sich nach C.-G. Rossby [11], der dies Problem überhaupt aufgeworfen hat, in einem 
Schachbrettmuster langer Wellen in der Westwindströmung der oberen Troposphäre aus. Es 
ist bemerkenswert und erfreulich, daß auch bei diesem ganz anders gelagerten Problem die 
Methoden der asymptotischen Stabilitätstheorie nach sinngemäßer Erweiterung sich durchaus 
bewährt haben. 

Des weiteren ist man in jüngster Zeit dazu übergegangen, nicht nur die Indifferenzkurve, 
also die Grenze zwischen angefachten und gedämpften Störschwingungen, und allenfalls noch die 
angefachten Störungen zu erörtern, sondern das volle durch die Störungsdifferentialgleichung 
und die Randbedingungen gegebene Schwingungsproblem zu untersuchen. Zu einer bestimmten 


2 
Wellenlänge } = = der Störschwingungen gehören bei einer bestimmten Reynoldsschen 


Zahl R vermutlich unendlich viele Schwingungen mit verschiedener Zeitabhängigkeit. Eine 
davon kann unter Umständen angefacht werden, während unendlich viele weitere zugehörende 
Eigenschwingungen gedämpft sein dürften. Zu diesem Fragenkreis haben C.C. Lin [12], 
H. Holstein [13] und W. Wasow [8] Bemerkungen gemacht, ohne daß indessen die ge- 
dämpften Eigenschwingungen in einem konkreten Fall diskutiert worden wären. Für die gerad- 
linige Couette-Strömung zwischen zwei parallelen Wänden hat bereits 1914 L. Hopf [14] das 
System sämtlicher Eigenschwingungen diskutiert, was vergessen zu sein scheint. Allerdings 
treten gerade in diesem Falle überhaupt keine angefachten Eigenschwingungen auf. Auf meine 
Anregung hat sich D. Grohne mit diesem Problem beschäftigt und neben notwendigen 
Korrekturen der Ho pfschen Ergebnisse bedeutende neue Einsichten in das System der Eigen- 
schwingungen bei der ebenen Poiseuille-Strömung und einer ebenen Strömung in einem diver- 
genten Kanal gewonnen. 

Trotz aller Erfolge der Methode der kleinen Schwingungen darf man nicht vergessen, daß 
damit wegen der Linearisierung nur das erste Stadium der Entstehung der Turbulenz erfaßt 
werden kann. Für die Beeinflussung einer Rohr- oder Kanalströmung durch endliche Störungen 
am Einlauf möchte ich auf die interessanten Untersuchungen von L. Schiller [15] nach- 
drücklich hinweisen. 


3. Über neuere Messungen der mittleren turbulenten Strömung in Wandnähe und im freien 
Strahl 


Wenn ich nun kurz auf die Untersuchungen über die ausgebildete Turbulenz zu sprechen 
komme, so darf ich zunächst auf die Experimente von H.Ludwiegund W. Tillmann[16] 
einerseits und von H. Reichardt [17] andererseits hinweisen. Diese Versuche haben zum 
Ziel, durch eine verfeinerte experimentelle Technik die wandnahe turbulente Strömung genauer 
zu erfassen, die für den turbulenten Hautreibungswiderstand und für den Wärmeübergang ent- 
scheidende Bedeutung haben. Für den letzten Zusammenhang vergleiche man die neue Arbeit 
vonH. Reichardt [18]. Was die freie Turbulenz anlangt, so liegt das Reichardtsche 
Forschungsheft [19] in zweiter Auflage vor, ferner ist eine nützliche Zusammenstellung aller 
einschlägigen Untersuchungen durch W. Forstalljr. and A.H. Shapiro [20] gegeben. 
In jedem Falle ist man bei den Experimenten über mittlere turbulente Strömungen bestrebt, 
zeitliche Mittelwerte der rasch schwankenden turbulenten Geschwindigkeiten und 
Drucke an einem ganz bestimmten Ort so gut, wie es bei der Größe der in die Strö- 
mung eingeführten Instrumente geht, zu messen. 


4. Turbulenzen verschiedener Stufe 


Dagegen werden in der letzthin mächtig aufblühenden sogenannten statistischen Theorie 
der Turbulenz räumlich-zeitliche Mittelwerte benutzt. Man kommt auf diese 
Weise zu dem fundamentalen Begriff von Turbulenzen verschiedener Stufe. 
Diese statistische Turbulenztheorie knüpft an Untersuchungen von G. I. Taylor [21]aus den 
Jahren 1935 bis 1938 an, die von Th.v. KärmänundL. Ho warth [22] weiter entwickelt 
wurden. Kräftige Impulse erhielt diese Forschungsrichtung durch Arbeiten von A.N. Kol- 
mogoroff [23], L. Onsager [24], L.Prandt1 [25], C.F.v. Weizsäcker [26] und 
W. Heisenberg [27]. Unter den vielen, welche das neu erschlossene Gebiet mit Erfolg 
weiter angebaut haben, sind wohl in erster Linie G.K.Batchelorund A.A.Townsend 
zu nennen, wobei ich die Verdienste anderer Forscher, die ich an dieser Stelle noch nicht nenne, 
keineswegs schmälern möchte. 

Tatsächlich kann man die Konzeption von Turbulenzen verschiedener Stufe viel weiter 
zurückverfolgen, bis in den Anfang der zwanziger Jahre, als Meteorologen wie A. Defant [28] 
undL.F.Richardson [29]sich mit Erfolg ähnlicher Vorstellungen bedienten, 
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Wenn man die Geschwindigkeiten und sonstigen Strömungsgrößen nicht nur über ein ge- 
nügend großes Zeitintervall, sondern auch über ein gewisses Raumgebiet mittelt, so wird die er- 
haltene mittlere Bewegung um so ausgeglichener, je größer man den räumlichen Mittelungsbereich 
nimmt. Bei Mittelungen über Raumgebiete verschiedener linearer Abmessungen / erhält man 
Turbulenzen verschiedener Stufe, wobei durch den Übergang zu größeren Mittelungsbereichen 
die Turbulenzen niedrigerer Stufe weggemittelt werden. Selbstverständlich bewirken aber die 
niedrigeren Turbulenzstufen einen Austausch der Eigenschaften, die bei den höheren 
Turbulenzstufen noch beobachtet werden; insbesondere interessiert man sich für den Austausch 
der kinetischen Energie zwischen den einzelnen Turbulenzstufen. L. Onsager spricht von 
einem Kaskadenprozeß, durch den Energie von den höheren auf die niederen Turbulenzstufen 
übertragen wird. L.F. Richardson hat diese Vorstellung in einen netten Reim gefaßt: 

Big whirls have little whirls 
Which feed on their velocity; 
Little whirls have smaller whirls, 
And so on to viscosity. 


Die Viskosität oder Zähigkeit mißt den Austausch, der durch die molekularen Bewegungen er- 
zeugt wird; sie bezeichnet nicht mehr eine Turbulenzstufe, sondern eben den Übergang von 
mechanischer in thermische Energie. A.Defant hatschon im Jahre 1921 den kühnen Versuch 
unternommen, den allgemeinen West-Ost-Kreislauf der Atmosphäre in mittleren geographischen 
Breiten als Turbulenzerscheinung großen Stils aufzufassen. Die Zyklonen und Antizyklonen 
walten als Elemente dieser Großturbulenz. 

Wenn man auch mit dem Begriff eines Austauschkoeffizienten, der in Analogie zu den 
molekularen Transporterscheinungen eingeführt wird, die komplizierten Verhältnisse bei den 
turbulenten Vermischungsvorgängen nur recht summarisch erfassen kann, so haben diese Vor- 
stellungen zu einem zweifellos großartigen Gedankenbild geführt, das eine Zusammenfassung sehr 
verschiedenartiger Erscheinungen ermöglicht. Es zeigt sich, daß der Austausch, der von Turbu- 
lenzen niedrigerer Stufe auf die nächsthöheren ausgeübt wird, in einer leidlich bestimmten Weise 
mit der Vergrößerung des Mittelungsbereiches zunimmt. Der kinematische turbulente Austausch- 
koeffizient sei mit e bezeichnet, wobei e wie die kinematische Zähigkeit » etwa in cm?/sec gemessen 
sei. Während sich für » in Luft unter normalen Bedingungen etwa 0,15 cm?/sec ergibt, hat man 
für den natürlichen Wind in den unteren 1000 m ein e von etwa 5 - 10° cm?/sec, wobei I etwa 
5.10*cmist. Für den durch die hochgetürmten Cumulonimben bewirkten Austausch, welchen 
uns die Wolkenbrüche bescheren, ergibt sich e zu ungefähr 10° cm?/sec, während / hier etwa 
3-10°cm ist. Für den Austausch durch Zyklonen und ähnliche Druckgebilde innerhalb des 
West-Ost-Kreislaufes erhält man ein e von etwa 101! cm?/sec, wobei man l angenähert mit 10° cm 
zu nehmen hat. Daß der Austausch mit l#/3 wächst, ist schon 1926 vonL.F.Richardson.c. 
festgestellt worden. Theoretisch wurde dieses Gesetz von A.N. Kolmogoroff,L. On- 
sagerund C.F.v. Weizsäcker zwischen 1941 und 1946 unabhängig voneinander herge- 
leitet, worauf ich noch kurz zu sprechen kommen werde. Beiläufig seibemerkt, daß C.F.v. Weiz- 
säcker l.c. diese Betrachtungen auf astrophysikalische Probleme angewandt hat und einen 
ganz gut stimmenden Schluß von dem Milchstraßensystem, dem unsere Sonne angehört, auf 
den Orion-Nebel gezogen hat. !hat dort die Größenordnung von 10!° cm. 


” 


5. Mathematische Hilfsmittel zur Erfassung der homogenen Turbulenz 


Ich komme nun auf die mathematischen Hilfsmittel zur Erfassung dieser Erscheinungen zu 
sprechen. C.F.v. Weizsäcker hat direkt Mittelungen über räumliche Bereiche verschiedener 
Abmessungen benutzt, A.N. Kolmogoroff gewisse Korrelationen zwischen Geschwindig- 
keitsdifferenzen. Am wirksamsten hat sich aber schließlich die Heranziehung der Fourier- 
Analyse erwiesen, die schon von G. I. Taylor [21] vorbereitet wurde. Physikalisch läuft das 
auf eine Spektralzerlegung der Turbulenz hinaus. 

Bisher hat man sich in diesem Forschungsbereich vorwiegend mit der sogenannten homo- 
genen Turbulenz beschäftigt. Man stelle sich etwa vor, daß die Flüssigkeit in einem großen Bottich- 
kräftig durchgerührt wird, so daß erhebliche turbulente Geschwindigkeitsschwankungen ent- 
stehen, ohne daß eine Strömung im Mittel sich ausbildet, d.h. 

db=0. 
b soll der Geschwindigkeitsvektor sein, durch Überstreichen werden Mittelwerte angedeutet. 
Der Einfluß der Wände des Bottichs sei dabei vernachlässigt. DiestatistischenEigen- 
schaften dieser turbulenten Strömung werden dann von den Ortskoordinaten unabhängig 
sein, was man eben gerade unter Homogenität versteht. Es bleibt eine Abhängigkeit von der Zeit 
übrig und allenfalls auch von der Richtung, falls man, um in unserem Bilde zu bleiben, beim 
Umrühren in dem Bottich eine Richtung bevorzugt hat. 
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Durch die Untersuchungen von G. I. Taylor [21] angeregt, führten Th.v. Kärmän 
undL.Howarth[22]JdenKorrelationstensorein,d.h.den Tensor der Korrelations- 
momente gleichzeitiger Geschwindigkeitsschwankungen an verschiedenen Orten (siehe Bild 6). 
Mit v, werde die Geschwindigkeitsschwankung in der Richtung e; an einer Stelle bezeichnet, die 
durch den Ortsvektor r gekennzeichnet ist, 


v; sei die Geschwindigkeitsschwankung in 
der Richtung e, an einer Stelle mit dem 
Ortsvektor r’. Mit r sei der Distanz- 

v.! 

d 
für den Korrelationstensor eine Beziehung 
einfacher Form erhalten. Schließt man 220g r3 
Richtungsabhängigkeit der statistischen G E-X 1) tXy 12 tXa 3 

Eine wesentliche Vereinfachung des Formalismus erreicht man durch den Übergang vom 

Korrelationstensor zum Spektraltensor, wiees J. Kampe& de Fe&riet [30] und 
G.K. Batchelor [31] taten: 


vektor 
=r—tr=&e+ &%+ &ses 
bezeichnet. Dann gilt für die Komponen- 
ten des Korrelationstensors wegen der 
Homogenität 

RW)=R(—t). 
Aus der Kontinuitätsgleichung kann man 
Eigenschaften aus, nimmt man also mit Ae = 7 2r ii. ı 
anderen Worten zusätzlich Isotropie ; Sa a2 23% 
an, so wurde diese Beziehung schon von R; (#) = vs 37 
Th. v. Kärmän aufgestellt. Unter der- ER VEN A EN 
selben Annahme konnte Th. v. Kärmän ar rn 
und L.Howarth [22]Jausden Navier-Stokesschen Differentialgleichungen eine dyna- 
mische Gleichung für den Korrelationstensor herleiten. 


ai = m I Rile) el Ir de, de.de,. 


[PX 
f ist die vektorielle Wellenzahl. Wegen der Homogenität ist der Spektraltensor hermitesch, d.h. 
n@)—1;0). 


wo die Tilde den Übergang zur konjugiert komplexen Größe bedeutet. Nach J. Kampe6de 
Fe&rietergibt sich aus der Homogenität und der Kontinuitätsgleichung für den Spektraltensor 
die einfache Form: 

PL 


= +). 


© = fhethet he, ist ein komplexer Vektor, der wegen der Kontinuitätsgleichung der Be- 


dingung Di 0 genügen muß; öist das Kroneckersche Symbol; g(f) ist eine reelle positive 
Funktion, % ist gleich |f. 
Führt man direkt die Fourier-transformierten der Geschwindigkeitsschwankungen ein: 


+ 
= as i ‚Jen ı UR, (iR, 
sm A 


wie es G. Darrieus [32] und W. Heisenberg [27] taten, so lassen sich, wie ich zeigen 
konnte, die Komponenten des Spektraltensors in folgender Weise ausdrücken: 


(Am)3!? 
H 
Dabei ist H das Raumgebiet, wo die Turbulenz als homogen angesehen werden kann. Voraus- 
gesetzt ist dabei, daß man ein solches Raumgebiet herausschneiden kann, während die Ge- 


schwindigkeitsschwankungen außerhalb H identisch Null gesetzt werden können, so daß also der 
Außenraum auf H nicht einwirken möge. 


Di = For: 
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Der Vorteil der Einführung des Spektraltensors liegt in der einfachen Form der kine- 
matischen und dynamischen Aussagen, wie sie aus der Kontinuitätsgleichung einerseits und aus 
den Navier-Stokesschen Differentialgleichungen gewonnen werden können, wozu man 
als reichhaltigste Quelle die Arbeit von G.K. Batchelor [31] vergleichen möge. 


6. Das Energiespektrum der homogenen isotropen Turbulenz 
Unter der zusätzlichen Voraussetzung der Isotropie kann man eine besonders einfache 
Spektralfunktion F(k) einführen, welche das Energiespektrum der Turbulenz charakterisiert: 
wobei k = |f| und 


v? N 
a BR 


sind (siehe Bild 7). H.L. Dryden [33] hat schon 1938 auf Grund von Messungen den Schluß 
gezogen, daß die Abnahme der kinetischen Energie in den verschiedenen Wellenzahlen der iso- 
tropen Turbulenz es erforderlich macht, daß die höheren Wellenzahlen Energie auf Kosten der 
niedrigeren gewinnen. Dies entspricht durchaus dem in Abschnitt 4 entwickelten Bilde von 
Turbulenzen verschiedener Stufe. Die Übertragung der kinetischen Energie der größeren 
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Bild 7. Schematisches Bild des Energiespektrums der homogenen isotropen Turbulenz bei quasi- 
stationärem Verhalten. Angedeutet sind theoretische eh von A.N. Kolmogoroff [23] 
bei mittleren und größeren k und eine Aussage von G@.K. Batchelor [31] bei kleinen k. 


Wirbel, deren Wellenzahl kleiner als ein gewisses % ist, auf die kleineren Wirbel muß bei 
quasistationären Verhältnissen von % unabhängig und schließlich gleich der 
Dissipation sein. Wenn wir diese Energieübertragung von den größeren auf die kleineren Wirbel 
noch auf die Zeit- und Masseneinheit beziehen und sie p, nennen, während die Dissipation mit 
bezeichnet wird, so drückt sich der obige Sachverhalt für quasistationäres Verhalten durch die 
Gleichung 9, = aus. Aus Dimensionsbetrachtungen folgt in einem mittleren Bereich des 
‚Energiespektrums, wo weder die Versuchsanordnung, welche die größeren Wirbel schafft, noch 
die Zähigkeit, welche die kleinsten Wirbel in Wärmebewegung überführt, eine Rolle spielt: 
F(k) m 92/3 K—513; 
denn die Dimensionen sind 
M=-PT=, [M-IL, W-L7S, 

wo T eine Zeiteinheit und Z eine Längeneinheit ist. Das obige Spektralgesetz ist von A.N. 
Kolmogoroff,L. Onsager und C.F.v. Weizsäcker mit verschiedenen Begrün- 
dungen aufgestellt worden. Aus ihm folgt, daß der Austausch mit der 4/3-Potenz der linearen 
Abmessung des räumlichen Mittelungsbereiches zunimmt, wie in Abschnitt 4 erwähnt wurde, 

Für abklingendehomogeneisotrope Turbulenz erhält man die Inte- 
grodifferentialgleichung: 


n k 
Mi [ ra. de (u tn) [2 7, Ddk, 
ö ö 


wenn man sich schon des summarischen Austauschkoeffizienten eg, für die Energieübertragung 
von den größeren auf die kleineren Wirbel bedienen will. Die Dissipation ist 


v2 KeF(k, 1) dk. 
N) 
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W.Heisenberg[27]setzt nun für den partiellen Austauschkoeffizienten: 
Ep -/[ (F, k) dk, 


wobei der Integrand, wie die Schreibweise andeutet, nur von F und k abhängen soll, so daß leicht 


aus Dimensionsgründen 
F 
Ey | Y Edi 
k 


mit einer Konstanten x folgt. (Für andere Ansätze, wie sie Th.v. Kärmän,A.Obukhov 
und L. Koväsznay für die Energieübertragung gemacht haben, vergleiche man etwa die 
Abhandlung von S. Goldstein [34]. Der Ansatz von Heisenberg empfiehlt sich u.a. 
durch seine Einfachheit.) Man hat nach „ähnlichen Lösungen“ der Heisenber gschen 
Integrodifferentialgleichung gesucht: 


k oo FB k 
3 F 
2 [va ) dk — 2 Vs dk +) kerdk, 
J ' 


d. h. nach Lösungen F, welche ihre Abhängigkeit bezüglich kim Laufe der Zeit t bis auf Maßstabs- 
abänderungen behalten. Es ergibt sich mühelos der Ansatz: 


Fk, ) = g(kye). 


Die dann resultierende Integralgleichung für g(x), wo kYt— x gesetzt wurde, istvonHeisen- 
berg [27] diskutiert und von J. Bass [35] sowie von S.Chandrasekhar [36] exakt ge- 
löst worden. Es ergibt sich aus diesem Ähnlichkeitsansatz, daß die kinetische Energie der Massen- 
einheit 


oo 


u [ra ) dk = 1 (0042 


wird. Eine Abnahme der Turbulenzenergie nach diesem Gesetz ist im Anfangsstadium 
tatsächlich [37] beobachtet worden, wenn sich auch im übrigen aus den Experimenten wie aus 
der Theorie manche Bedenken gegen die durchgängige Gültigkeit des Ähnlichkeitsansatzes für 
das Turbulenzspektrum ergeben haben. Im Endstadium, wenn die turbulente Energie- 
übertragung gegenüber der molekularen ganz zurücktritt, wird nach G.K. Batchelor und 
A.A. Townsend [38] 


F(k,t)z Ame —2r Rt), 


A ist dabei die Loitsianskysche Invariante [39], d.h. eine Größe, die bei der homogenen 
isotropen Turbulenz im Verlauf des zeitlichen Abklingens einen konstanten Wert behält. t, ist 


4 ; 2 
eine passend zu wählende Zeitkonstante. Es folgt dann, daß 5 im Endstadium proportional zu 


(t —t,)®'? wird, wie es in [38] auch experimentell bestätigt wurde. 

In Verfeinerung dieser Energiebetrachtungen kann man auch die Abnahme der Energie 
einer einzigen Geschwindigkeitskomponente ins Auge fassen. Unter ausschließlicher Berück- 
sichtigung der Kontinuitätsgleichung und der Tatsache, daß die mittlere Strömung im Sinne 
von OÖ. Reynolds verschwindet, erhält man 


Zeaeton  1[2 De | 5 
a a | Frude. 


Mit p werden hier die Druckschwankungen bezeichnet. G.K. Batchelor [31] hat unter 
Hinzunahme der Homogenität das Äquivalent dieser Gleichung in Fouriertransformierten 
angegeben: 


yuH+ 


TO _ nr. 


(In der Bezeichnungsweise RE wirnicht der irreführenden Konvention, daß doppeltes Auf- 
treten eines Index eine Summation über diesen nach sich zieht.) Führt man nun auch in der vor- 
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hergehenden Gleichung die Homogenität ein, so verschwinden die Trägheitsglieder auf der 
rechten Seite, mit anderen Worten, die turbulente Energieübertragung 


ae 0 —— 0) 

2 | 8 N 
- v; v ne |: 
> 0x, 1 2. ı 72 2 ı °3 


so daß für die entsprechende Fo uriertransformierte gilt: 
© „, 
II] && dr, dr, dk, = 0. 


Durch diese turbulente Energieübertragung als Wirkung der Trägheitsglieder erfolgt also ein 
Energiefluß zwischen verschiedenen Wellenzahlen für dieselbe Ge- 
schwindigkeitskomponente. Andererseits sieht man durch Heranziehung der 
Kontinuitätsgleichung, daß die Druckglieder bei Summation über alle drei Komponenten ver- 
schwinden: 


und 


ee 


Diese Glieder erzeugen daher einen Energiefluß zwischen verschiedenen Ge- 
schwindigkeitskomponentenbeiderselben Wellenzahl. Dieser zweite 
Energiefluß dürfte im allgemeinen auf gleichmäßigere Verteilung der Energie auf die Geschwindig- 
keitskomponenten, also auf Isotropie hinwirken, die ja hier nicht vorausgesetzt war. Am einfach- 
sten zu beschreiben ist die Wirkung der Zähigkeit, die durch das Glied 


vude; bzw. —2vKk2Tilt) 


gekennzeichnet ist. Die Zähigkeit beeinflußt also jede Geschwindigkeitskomponente und jede 

Wellenzahl unabhängig von den übrigen, und zwar umso mehr, je größer die Wellenzahl % ist. 

Die Summenbeziehung für die Druckglieder und die Integralbeziehung für die Trägheits- 

glieder sind an die Homogenität der Turbulenz ‚gebunden. Bei inhomogener Turbulenz würden 
n2 


diese Glieder eine zusätzliche Änderung von 2 in der Größe 


2 
— div 2 + nn ) 
0 2 


hervorbringen. 


7. Experimentelle Prüfung der Aussagen über die homogene Turbulenz an der 
Windkanalturbulenz 


- Die Aussagen über die homogene Turbulenz werden an der Windkanalturbulenz geprüft, 
wo die erstrebte gleichmäßige Hauptgeschwindigkeit in der z-Richtung U sei. Man kann an- 
nehmen, daß das zeitliche Schwankungsglied über nicht allzu große Distanzen nahezu unver- 
ändert mit U stromabwärts gespült wird. Die Voraussetzungen für die Anwendbarkeit der 
theoretischen Ergebnisse für die homogene Turbulenz sind in der Literatur nicht einwandfrei 
formuliert worden. Es scheint mir auf folgende Bedingung anzukommen: Wenn 7 die Mindest- 
zeit ist, die zur Gewinnung zuverlässiger statistischer Aussagen gerade genügt, so dürfen sich die 
zeitlichen Mittelwerte der Windkanalturbulenz auf einer Strecke Ur nicht merklich ändern, was 


man als lokale Homogenität bezeichnen mag. Dann kann der Operator — in der homogenen 


% 
Turbulenz durch den Operator un in der Windkanalturbulenz bei der Anwendung auf sta- 


tistische Eigenschaften ersetzt werden und die homogene Turbulenz kann in einem Windkanal 
hinter einem ‚‚Turbulenzgitter“, einem aus Stäben aufgebauten Rost, der etwa quadratische 
Maschen bildet, studiert werden. Der Chronist mag mit ironischer Befriedigung verzeichnen, 
daß die Windkanalturbulenz, mit der die experimentierenden Aerodynamiker bei ihren Modell- 
versuchen so große Not hatten, jetzt”zu einem wertvollen Substrat grundlegender Versuche 
in der Strömungslehre geworden ist. 

Durch die neueren Arbeiten von I. Proudman [40] und R.W. Stewart and A. A. 
Townsend [41] ist der Vergleich der Experimente mit den bisherigen theoretischen Ergeb- 
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nissen in breiter Weise und mit kritischer Vorsicht durchgeführt worden. Es zeigt sich, daß 
manche Züge der Experimente richtig wiedergegeben werden, während andererseits systematische 
Abweichungen von grundsätzlicher Bedeutung nicht zu verkennen sind. Besonders die Er- 
haltung der Form des Turbulenzspektrums, die in Abschnitt 6 als Annahme eingeführt wurde, 
um eine Lösung der Heisenbergschen Integrodifferentialgleichung zu ermöglichen, scheint 
nur sehr bedingt zu stimmen. 


8. Lokale Isotropie und Ähnlichkeit im Kleinen 


Diese homogene isotrope Turbulenz scheint auf den ersten Blick eine sehr eingeschränkte 
Bedeutung zu haben, ja geradezu einen ausgearteten Fall von Turbulenz darzustellen. Eine 
universale Bedeutung wird diesen Turbulenzerscheinungen aber durch Hypothesen zugesprochen, 
die von A.N. Kolmogoroff herrühren. Danach sollen die Bewegungen der Wirbel, deren 
Abmessungen kleiner als eine gewisse Länge sind, universalen Charakter haben, gleichgültig, wie 
die Turbulenz im großen beschaffen ist. Die statistischen Eigenschaften der Turbulenz im 
Kleinen sollen isotrop sein und durch 9 (die Energiedissipation für die Zeit- und Masseneinheit) 
und die kinematische Zähigkeit » vollständig und allgemein bestimmt sein (siehe Bild 7). 


Diese großartige Konzeption von Kolmogoroff ist bis zu einem gewissen Grade 
durch Messungen von A.A. Townsend [42] bestätigt worden. Von der Zufallsverteilung 


2 die wohl hauptsächlich durch die Feinstruktur der Turbulenz bestimmt ist, wurden die 
Ti 


ersten vier Monate gemessen und daraus dimensionslose Größen gebildet. Es ergab sich 


2 (=) (2) 3/2 
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Die angegebenen Zahlenwerte wurden sowohl hinter einem Turbulenzgitter mit quadratischen 
Maschen als auch im Windschatten hinter einem einzigen Draht mit leidlicher Übereinstimmung 
festgestellt. S ist die ‚‚Schiefe‘‘ (skewness) der Zufallsverteilung (siehe Bild 8), 7, ist der „‚Ab- 
flachungsfaktor“ (flattening factor), wobei dieser letzte Name wenig glücklich gewählt ist. Je 
steiler eine Verteilung nämlich verläuft, um so größer wird gemeinhin T,. 7, —3 ist der „Exzeß‘“ 
(measure of kurtosis), da 7, für eine 4 
Gaußsche Normalverteilung 3 ist. ° 
Die in dem Bild 8 aufgezeichnete 
Zufallsverteilung von = ist in der 

ı 
Mitte steiler als eine Gaußsche 
Normalverteilung, wie sie etwa v; be- 
folgt; an beiden Enden läuft sie mit 
längeren Schleppen aus. 


Die Schiefe 8 (3) tritt in der 


Wirbelgleichung auf. Wird der Wir- -5 


belvektor ur ung Sei 
rotd = W= WC + Wat WyEz 
gesetzt und der Wirbelbetrag Bild 8. Statistische Verteilung von v, (gestrichelt) = (ausgezogen) nach 
|| 00: Messungen von A.A. Townsend [42], [2]. Die Maßstäbe sind willkürlich. 


so gilt in homogener Turbulenz 


a nr) 2 3 >) 
Ir = 2m(mVb) —2» 2) 2, im‘ 


i=1j-1l 
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Unter Hinzunahme der Annahme der Isotropie kann man diese Gleichungnach G.K.Batche- 
lorund A.A.Townsend[43]so ausdrücken: 


9w? 7 Teo\ 82 (=) 14 ‚— 3/2 @ 
— = — —— (02)”” 8\3,.) — —— (02)"" —, 
ot 3y5\ 02; ay5\ R, 
wo 
Kar 0°», 2 0v, 2]? 
ul Ey & 
Tanss 
Ri= ih \ (2) 
sind. 


Das Glied mit der Schiefe S mißt den von Helmholtz eingehend diskutierten Effekt, 
daß durch die Dehnung der Wirbelfäden eine Vermehrung des Wirbelbetrages eintritt. Das Glied 
mit @ mißt die Diffusion der Wirbelung durch Reibung, die in den Experimenten [43] stets größer 
ausfällt als das Helmholtzsche Glied und zwar um 35% bis 100%. 

Um auf die Kolmogoroffsche Hypothese der lokalen Ähnlichkeit zurückzukommen, 
muß ich erwähnen, daß gerade in letzter Zeit auf Grund experimenteller Befunde [41] ein- 
schneidende Beschränkungen für den Gültigkeitsbereich dieser Hypothese angemeldet sind. 
Doch halte ich die Diskussion darüber nicht für abgeschlossen. 


Schlußbemerkung 


Wenn wir dies in großen Linien entworfene Bild der heutigen statistischen Turbulenz- 
theorie überblicken, so müssen wir feststellen, daß bislang wenig quantitative Aussagen gewonnen 
wurden, wenn man von reinen Dimensionsüberlegungen absieht. Doch haben diese theoretischen 
Vorstellungen, wie gerade die zuletzt erwähnten Untersuchungen zeigen, in Verbindung mit Ex- 
perimenten zu einer tieferen Einsicht in den Mechanismus turbulenter Strömungen, in ihre 
Feinstruktur, geführt. Es hat nicht den Anschein, als ob die Kraft dieser Ideengänge schon er- 
schöpft wäre. Freilich haben diese Anschauungen bislang noch keinen Beitrag zu dem Problem 
der mittleren Geschwindigkeitsverteilung im Sinne von OÖ. Reynolds geliefert. Selbst der 
Abfall der Turbulenzenergie hinter einem Gitter in einem Windkanal konnte bisher nicht bündig 
theoretisch erfaßt werden. Es offenbart sich hierin eine eigentümliche Schwäche dieser For- 
schungsrichtung, so daß man bis auf weiteres lieber von einer statistischen Analyse als von einer 
statistischen Theorie sprechen sollte. 
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KLEINE MITTEILUNGEN 


Über das Ausgleichen von Kurven in Para- 
meterdarstellung. 

1. Einführung. In der x,y-Ebene sei eine 
Punktreihe P, (z,, y,),v =1,2,...n, durch Messung 
der Koordinaten (2,, y,) festgelegt. Die fehlerfreien 
Punkte P* (z*, y*) sollten auf Grund einer Hypo- 
these oder auf Grund von physikalischen oder tech- 
nischen Überlegungen auf einer Kurve 


EXERYTTAREBINEN 0 SE) 


liegen, von der man zwar die Form F aber nicht die 
Parameter A*, B*, ... kennt. Man pflegt dann mit 
Hilfe der Ausgleichsrechnung geeignete Parameter 
A, B,...so zu bestimmen, daß die mit A, B,... be- 
rechnete Kurve den Verlauf der Punktreihe ‚‚mög- 
lichst gut“ wiedergibt. Die allgemeine Lösung dieser 
Aufgabe scheint für den Fall, daß der Zusammenhang 
zwischen x und y nicht in der Gestalt (1.1) sondern in 
Parameterform 


2214, Bee) unde u — ut AwB, sel.) 
gegeben ist, noch nicht bekannt zu sein. Kürzlich 
wurde diese Aufgabe von Jaeckei behandelt }), 
jedoch mit der Einschränkung, daß in dem gemessenen 
Wertepaar (z, y) eine Koordinate fehlerfrei bestimmt 
wordenist. Da nun die Ergebnisse des Ausgleichens, 
nämlich die Zahlenwerte der Parameter A, B,... und 
die geglätteten Koordinaten (%,, 9,) durchaus davon 
abhängig sind, ob nur die x-Werte, nur die y-Werte 
oder beide fehlerhaft sind, so ist es nicht gleichgültig, 
unter welcher Annahme man die Rechnung durch- 
führt. Man ordnet zwar häufig die Meßfehler einseitig 
einer der gemessenen Veränderlichen zu, mußsich aber 
dann darüber klar sein, daß man die Ergebnisse damit 
von vornherein in nicht übersehbarer Weise fälscht 
und oft den eigentlichen Zweck des Glättens — die 
Herabsetzung des Einflusses der Meßfehler auf das 
Ergebnis — gar nicht erreicht. Wir wollen deshalb im 
folgenden die Ausgleichaufgabe unter sehr allgemeinen 
Voraussetzungen lösen, wenn der funktionale Zusam- 
menhang zwischen x und y in Parameterform gegeben 
ist. 

2. Aufgabenstellung. 
etwas verallgemeinerter Fassung) 


G(z,1|4,B,...)=0 und H(y,1]4,B,...)=0 (2.l) 


Es sei also (in 


die Parameterdarstellung einer ebenen Kurve 
Fis yasıBas)=D,. 


Die Funktionen G=0 und H=0 sind der 
Form nach bekannt. Von allen drei Kurven F, @ 
und H setzen wir voraus, daß sie frei von singu- 
lären Punkten sind. Die Elimination des ‚‚veränder- 
lichen‘‘ Parameters t aus G = 0 und H = 0 sei nicht 
(oder nur mit erheblichem Rechenaufwand) möglich, 
sodaß der sonst übliche Weg des Ausgleichens über 
F = 0 nicht gangbar ist. Neben dem veränderlichen 
Parameter i enthalten € = O0 und H = O.noch die un- 
bekannten Parameter A, B,..., die wir als ‚feste‘ 
Parameter bezeichnen wollen. Gegeben ist eine aus 
genügend vielen Punkten bestehende Meßreihe 
P,(&, y,), so daß bei der Berechnung von A, B,.. 
wirklich ein Ausgleichsproblem vorliegt. Die Genauig- 
keit bei der Messung der Punkte P sei durch die 
Streumatrix 


Oyz Oyy 0% 0y0 c, 


gegeben. Ihre Elemente o%,,... sind normalerweise 
Funktionen von (z, y)also 045 =0gr (2; y),.... Nur 


1) K. Jaeckel: Fehlerausgleichung bei Funktionen in Para- 
meterdarstellung. Z. angew. Math. Mech. 31 (1951) S. 185. 


wenn die Punkte P in der ganzen Ebene mit gleicher 
Genauigkeit gemessen werden können, ist 8 — konst. 
Sind die Meßfehler von x und y unabhängig vonein- 
ander, so ist der Korrelationsfaktor 0 = (0, und die 
Matrix reduziert sich auf ihre Spur. 


3. Lösung der Umkehraufgabe. Wir 
lösen zunächst die einfachere Aufgabe, bei bekannten 
Parametern A, B (deren Zahl wir auf zwei beschränken 
wollen) und gegebenem Meßpunkt P,(2,, y,) den aus- 


geglichenen Punkt P,(%,, y,) zu finden. Der Kürze 
wegen wollen wir den Index» in der Zwischenrechnung 
weglassen, wenn kein Mißverständnis möglich ist. 
Für die an den Meßwerten (z, y) anzubringenden Ver- 
besserungen Aa =z2— x und Ay=y-—y gelten 
(wie an anderer Stelle?) gezeigt wurde) die Glei- 


nd OxFr + oyoFy _F (31) 
O% (0 F2)? + 2020, 0 FzxFy + (oyFy)? l 


o%  (onFr)? +2050y0 FzFy + (oyFy)? : 
wobei #F =F(x, y| A, B)=F, den an der Meß- 
stelle P, berechneten normalerweise kleinen Funk- 
tionswert darstellt. In dieser Gestalt sind die Glei- 
chungen natürlich noch nicht verwendbar, da wir 
nach Voraussetzung die Gestalt der Funktion F nicht 
kennen. Wir entnehmen aber aus (3.1) und (3.2) 
die Beziehungen 


Az 
„ir = (0,F5 + oyoFy) | 
c 
und r . (3.3) 
Y 
mit vorläufig unbekanntem A, 
RL F ala 3.4), 


(Or Fz)’+ 205040 Fa Fy + (oyFy)? 
Für die ausgeglichenen Werte (7, y) soll gelten 


G(2,1|4,B)=0 und H(y,1|4,B)=0 (3.5). 


Setzen wir in diesen Gleichungen die Meßwerte x, y 
und einen passend dazu bestimmten Näherungswert ! 
ein, so werden @ (z,1|4, B)=@ und H (y,t|4,B)=H 
im allgemeinen nicht verschwinden. @=G, und 
H = H, sind normalerweise kleine nur wenig von 0 
abweichende Werte. Entwickelt man die Funktionen 
G = 0 und H = Oin Potenzreihen, so erhält man aus 
(3.5) 


G(2,1)A, B) = @(z,t|A, B) 
+E- Et Ya+ı-=0 
und 
H(%, 1A, B) = Hy, t|A, B) 
+@# -WH,+t-DH+t=0. 
Mit der beim Ausgleichen üblichen linearen Näherung 
gilt also 
Ar +GA+eE=0, | 
H,4y+Hı A+H=0 J 
Setzt man in diese beiden Gleichungen Ax und Ay 


aus (3.3) ein, so erhält man zwei lineare Gleichungen 
zur Bestimmung der Unbekannten A und 4t, 


x G2(0&Fz + oyoF,)A+HAt+G=0, Ion 
2) K. Stange, Über das Ausgleichen einer fehlerhaften li- 


nearen Punktreihe bei korrelativer Verknüpfung der Meßifehler. 
Mitt.-Bl. für Math. Statistik. 1052 S.48. - - - 
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Sofern die Determinate dieses Systems nicht ver- 
schwindet, wird insbesondere 
Ze 
GH, — HG, 
(@% Hr 62) (oa Fz+oyeFy) — (HyGroy) (oyFy-+ 0% 0F%) 
(3.8). 
Setzt man diesen Ausdruck für A in die erste Glei- 
chung (3.3) ein, so kommt zunächst 


Az 
rim 
(oFz% + oyoFy) (@ HH — HG,) 

(03 Ht0z)\ 05 Fz + 0y0Fy) — (HyGroy) (oyFy + 020 Fr) 

: (3.9). 


Hier eliminieren wir mit Hilfe von 
et _#: . a —_ GaH: 
Fy x Hy Fi Gx H, G 


die partiellen Ableitungen nach F. Zweckmäßig setzt 
man abkürzend 


G, Hr; 0% = & und H,G;: Oy= w bes (3.10). 
Dann kommt schließlich 
da _ ®—oYV 


€ (GH, — HG,) (3.11) 


% M—-20209Yy + PR 
und entsprechend 

are ri reD 
oy DB —-209Y +? 


Damit sind die an den Meßwerten z,, y, anzubringen- 
den Verbesserungen Ax,, Ay, bestimmt. Sie sind ab- 
hängig von der Streumatrix $, an der Meßstelle und 
von den für %,, 4», tv berechneten Werten der Funk- 
tionen G@ und H und ihren partiellen Ableitungen. 

Ist insbesondere bereits @(z, {!4,B)=0 und 
H(y,t|4,B) =0, so liefern die Gleichungen richtig 
Az =0Ound Ay =0. 

Es bedarf noch einer kurzen Erörterung, unter wel- 
chen Bedingungen die Lösung entarten kann. Der 
Korrelationsfaktor o ist dem Betrage nach kleiner als 
1, also 0®< 1. Daß dann der Nenner in (3.11) nor- 
malerweise nicht verschwinden kann, sieht man fol- 
gendermaßen ein. Ist etwa ®+ 0, so wird der Nenner 
D°{1 — 20 (Y]D) + (W/D)®%}. Die geschweifte Klam- 
mer ist eine quadratische Funktion von Y/®, die an 
der Stelle Y/® = o das Minimum 1 — 0? > 0 besitzt. 
Für Y=0 schließt man entsprechend. Nur wenn ® 
und Y gleichzeitig verschwinden, entartet die Lösung. 
Schließt man den trivialen Fall verschwindender 
Streuungen 0, =0 und oy=0 (bei dem kein Aus- 
gleichen notwendig ist) aus, so kann der Nenner 
wegen (3.10) nur 0 werden, wenn gleichzeitig @, H;—= 0 
N; Hy, @ =0ist. Dann muß entweder gleichzeitig 
(Gx2 =0 und H, = 0) oder (G, = 0 und 6; = 0) oder 
(H, —=(0 und H = 0) oder (H; =0 und G = 0) 
sein. Alle vier Fälle sind aber unmöglich, wenn 
G(z,t)=0, H(y,t)=0 und F(x, y)=0 keine singu- 
lären Punkte enthalten. 

4. Die Ausgleichbedingung. Nach den 
vorbereitenden Erörterungen des vorausgehenden Ab- 
schnitts ist die Lösung unserer eigentlichen Aufgabe 
einfach. Gegeben ist jetzt eine ganze Meßreihe P,. 
Gesucht werden die festen Parameter A, B. Als Aus- 

leichprinzip wählen wir die folgende Forderung: 
ir bestimmen die Parameter A, Bso, daß die Wahr- 
scheinlichkeitsdichte an den gegebenen Meßpunkten 
P, (2,, y,) in bezug auf die gesuchten ausgeglichenen 
Punkte P,(%, Y,) für die gesamte Meßreihe mög- 
lichst groß wird. Dann ist offenbar der Widerspruch 
zwischen der Meßreihe P, und der ihr zugeordneten 


Reihe P, am kleinsten. Ist die Streumatrix 3 unab- 
hängig von den Parametern A, B und sind die Meß- 


(6HK— HG, (3.12). 
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fehler normal verteilt, so führt die eben gestellte For- 
derung zu 


n 
ı J/A2\2 A 2 
EEE Hun 
al ®, %/v 9% /,\0y )» oY .) 
(4.1). 
Da nun zwischen den Koordinaten der gemessenen 
und der ausgeglichenen Punkte die Beziehungen (3.11) 
und (3.12) bestehen, so setzen wir Ax und Ay aus die- 


sen Gleichungen in (4.1) ein. Nach kurzer Rechnung 
nimmt dann die Ausgleichbedingung die Gestalt 


Tal (@, H:y er H,G; »)° 
rc -le200% +7 


| -Min (4.2). 


an. Bei der Bildung der partiellen Ableitungen nach 
den Parametern A, B haben wir G, = @, (2, {»| A, B) 
= @,(A, B) und sämtliche anderen Größen als Funk- 
tionen dieser Parameter aufzufassen. Die weitere 
Rechnung verläuft in bekannter Weise und führt 
schließlich zu den Normalgleichungen zur Berechnung 
von A und B. In verwickelten Fällen muß man natür- 
lich die Ausgleichbedingung durch Einführung von 
Näherungswerten A, und B, in bezug auf die Para- 
metera = A — A, und b = .B — B,linearisieren. 
Ist insbesondere @=x2—1=0, so wird H (y, x|A, B) 
=(. Mit@,=1,% =—1lund A; = H,nimmt dann 
die Ausgleichbedingung (4.2) die gewohnte Form 
für das Ausgleichen der Funktion H (x, y|A, B) = 0 


an, nämlich 
Hr ; 
a a li =. Min 
(Hz 0%)? +29 (H202)(Hyoy)» + (Hyoy)} (4.3) 


5. Beispiel. Als Beispiel sei die Ausgleich- 
bedingung aufgeschrieben, wenn man aus einer Meß- 
reihe P, den Halbmesser A des Rollkreises einer ge- 
wöhnlichen Zykloide 


G=x—Alt—sint) und H=y—A(1-—cost) (5.1) 
bestimmen soll. Die Meßgenauigkeit sei in der ganzen 
x, y-Ebene die gleiche und die Meßfehler seien unab- 
hängig voneinander, also 

O0 = 0, — konst., Op = oy = konst, » =0 (5.2). 


Dann wird ®= — Ao,sintund Y= — Aoy,(1— cost). 
Die Ausgleichbedingung (4.2) lautet nach einfachen 
trigonometrischen Umformungen 


2 
n I" sin ni — 1, C08 2) +4 (" cos 2 — 2sin 3) 


v=1 07, c08? nu + °, sin? u 
— Min (5.3). 
Sie liefert für 0, = 0, = o als besten Wert für den 
Halbmesser A des Rollkreises . 
ty RR» ty . ty 
(# 0085 —yrsin a (" 0085 — 2 sin 2) 
PARSE Seen N — —_—______ 2 (8.4). 


iv . W\2 
((e 0085 — 2sin 3) | 


Karlsruhe. K. Stange. 


Zur Theorie der Profile geringer Dicke und 
Wölbung 
Einleitung 


In der Theorie dünner, schwach gewölbter Trag- 
flügelprofile kennt man seit langem Verfahren, die die 
Strömung um ein solches Profil mit Hilfe einer im 
Innern gelegenen Singularitätenverteilung darstellen. 
Nach Ackermann-Birnbaum!) läßt sich das 


ı) W.Birnbaum, Die tragende Wirbelfläche .... Z. angew. 
Math. Mech. Bd.3 (1928) 8. 290. 
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Profil in erster Näherung durch eine Wirbelschicht 
längs der Mittellinie ersetzen; in zweiter Näherung 
kann man durch Hinzunahme einer Quellbelegung 
nach Föttinger?) auch die Dicke berücksichti- 
gen. Die Handhabung ist später verbessert worden 
von Keune?°) und Helmboldt) und insbeson- 
dere von Riegels?°) für die Praxis zu einem Ab- 
schluß gebracht worden. Im Folgenden will ich ver- 
suchen, die Darstellung der Theorie zu ver- 
einheitlichen. Ausgehend von den funktionentheore- 
tischen Ansätzen, unter Beschränkung auf Glieder 
erster Ordnung in Dicke und Wölbung, läßt sich bei 
Benutzung elliptischer Koordinaten in der Tat eine 
übersichtliche Darstellung finden, die sich in Einfach- 
heit und Form an die von Ackermann-Birn- 
baum und Föttinger eng anschließt. 


Grundgleichung und ihre verein- 
fachte Form für Profile geringer 
Dicke und Wölbung 


Das komplexe Strömungspotential einer Quell- und 
Wirbelbelegung auf einem schwach gewölbten glatten 
Bogen An (£) durch die Abszissenpunkte &= — 1 und 
= --1, in Verbindung mit einer Parallelströmung 
mit den Geschwindigkeitskomponenten U und V, ist 
gegeben durch 


+1 
En 


’ 


wobei g(£) die Quelldichte, y(£) die Wirbeldichte auf 
dem Bogen an der Stelle & bedeuten, die nur im Innern 
des Intervalls -—1<&<+ 1lals von Null verschie- 
den angenommen werden. 

Schreiben wir 


y-ınla)+Y . (2), 


wobei wir Y. von der KurveAn(z) ab zählen, so können 
wir folgende Abspaltung vornehmen: 


D=(U —iV)(e+iin(e)+iY) 
+1 


a IK +iy(H)n(e +iY — Has 
—l 
+1 

+ [a +iy) In 


= 


2—L 


Ra rer 


Beachten wir dabei, daß im letzten Ausdruck der Lo- 
garithmus auf die Form gebracht werden kann 


mil N share) 
() eli4 ee 


die sich nach Potenzen von A entwickeln läßt, wobei 
die mit den Reihengliedern zu bildenden Integrale auch 
für Y—0 existieren, so ergibt sich für ihn in linearen 
Gliedern in A: 


Sp 


| ee E-iNdet.. 


—1 


®2) H. Föttinger, Fortschritte der Strömungslehre.... 
Jahrb. d. Schiffbautechn. Gesellschaft 1924, S. 295. 
®) F.Keune, Die ebene Potentialströmung’ um allgemeine 
dicke Profile. Jahrb. d. deutschen Luftfahrtforschung 1938, 8. 13. 
%)H.B.Helmbold-F.Keune, Mehrere Aufsätzein Luft- 
fahrtforschung Bd. 20 (1943). 
*», F. Riegels, Das Umströmungsproblem bei inkompres- 
siblen Potentialströmungen. I. Mitteilung Ing.-Arch. Bd. XVI 
(1948). II. Mitteilung Ing.-Arch. Bd. XVII (1949). 
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Für die Stromfunktion liefert dieses Integral die An- 
teile 


die jedoch für die Nachbarschaft der Belegungslinie 
als Größen zweiter Ordnung fortzulassen sind, denn 
für dünne Profile ist die Quelldichte g(x) und die 
Ordinate Y als klein von erster Ordnung anzusehen. 
Für die Bestimmung der Profilstromlinie genügt es 
daher, die Darstellung 


B=(U—-iV)(r +iAn(e)+iY) 
+1 
EN („+iy)n(e—E+iY)de... (3) 


zu betrachten, die mit A = 0 für eine Streckenbele- 
gung genau gilt. Das hierbei auftretende Integral 
können wir nun mit Hilfe einer bereitgestellten Rei- 
henentwicklung für den Kern in elliptischen Koordi- 
naten auswerten ®). 


Auswertung durch Reihen- 
entwicklung 
Setzen wir 


2= —chs-c890; &= — c08%; 


Y=she-sin®; ch = 6]; | (4) 
0<O<n; 0<9<a; ei. 
OEL sh = ©in; 


und machen den Ansatz 


5 % 
We 
re m an 
Ay, = Ayr + tAyni 
so finden wir unter Berücksichtigungder Darstellung ®) 


In(e —E+iY)=e—i0+ const 


cand )...46) 


x PB nei), a 


n=1 

die Entwicklung 

®= (U —- iV)(- che- cos @ + iAn(x(e, 9)) 
+ishe-sin o)+ Ula(e — i +. const) 


An ) 

n 

Daraus ergeben sich die Stromlinien, unter ihnen das 
Profil e= (©), aus dem konstant gesetzten Imagi- 


närteil: 


U (An(z(s, 0)) + sh e-sin @)-+ Veh e: cos © 


(M). 


>‘ En e mrle—io), 
n=1 


gt — ne j 
77 (a E— AO — e —._ (@nicosn © (8). 


n=1 
+ Ayr sinn o)) =Im® = /(s, ©) 


Man erkennt aus dieser Gleichung, daß 
dgr = 0 . 


. (9) 


°%, K. Jaeckel, Ermittlung einer Reihendarstellung.... 
Z. angew. Math. Mech. Bd. 30 (1950) $. 184. 
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vorausgesetzt werden muß, wenn das Profil geschlos- 
sen sein soll. Soll weiterhin im Punkte & = + 1 Pro- 
filhinterkante und Staupunkt zusammenfallen, so 
müssen dort die Ableitungen der Stromfunktion nach 
© Er y oder auch nach e und © verschwinden. 

it 


— = U(—2n'(z(e, ©)) sh ecos © 
+ chhesin 0) + Vshecos® + Ula,; 
u > e"e(qn; cosn © + Ay, sinn 9)) » 
ov KB ‚ 
rn U (An’ (© (&,)) ch esin@+ sh e cos ©) 


_ Vchsin@- UN ot. 


n=1 
-(— 4y;5inn © -+ ay,cosn ©) } 


folgen daraus die Abfluß-Bedingungen (e= 0,0 =n) 


2 Ayr(— 1)* OR Gyr = O wegen (9) 
=i 


Z mi(—1)"=0 . (10). 
n=0 


Für genügend kleine Werte von eläßt sich die Strom- 
funktion nach e entwickeln. Es ergibt sich in Gliedern 


erster Ordnung 


An(& (0, 9)) 477 cos (0) 


- I (m osno + m im no) 
n n 
n=1 

: Y(e, ©) 
N re 


Stil) 


wobei angenommen wird, daß auch der Anstellwinkel 
V/Uklein sein soll. Dann nämlich sind alle a,,; und ay, 
von der ersten Ordnung klein und sorgen dafür, daß 
die Korrekturglieder vernachlässigbar klein werden. 

Innerhalb unserer Näherung können wir weiterhin 
An (x) als Mittellinie des dünnen Profils annehmen, da 
kleine Abweichungen davon auf das Ergebnis keinen 
Einfluß haben. Daher läßt sich e(©) für das Profil 
als gerade Funktion in © ansetzen. 

Dann wird die Kontur bestimmt durch die Glei- 
chung (11) mit 


WW 
m (<(9), ©) > 10 ’ 


die in die folgenden beiden zerfällt 


An(x(0, ©))+ (7 ee a) 008 @ 


oo 
Ani 
— > —, cos n®=(\,, 


— (12), 
a 
e(©)sin O — — sinnO=0, 
x(0, 9) = — cos © 


die besagen, daß die Wirbelschicht nach der Acker- 
mann-Birnbaumschen Theorie berechnet wer- 
den kann und die Quelldichte nach dem Föttin- 
gerschen Verfahren, bei dem man das Profil in ellip- 
tischen Koordinaten darzustellen hat. 


Bestimmung der Verteilungen und 
Geschwindigkeiten für ein 
gegebenes Profil 
Für ein vorgelegtes Profil geringer Dicke und Wöl- 
bung bestimmt man zunächst die Mittellinie und da- 
durch das symmetrische Profil, dessen Hinterkante 
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dem Punkte & = + lentspricht. Der Punkt& = — 1 
liegt im Innern dieses symmetrischen Profils und ist 
zweckmäßig so zu wählen, daß sich die Profilnase mit 
Hilfe des Krümmungskreises einer konfokalen Ellipse 
mit den Brennpunkten (— 1,0) und (+ 1,0) am besten 
anschmiegt, womit in den Normalfällen dem Ansatz (5) 
Rechnung getragen ist. Man entnimmt dann dem 
elliptischen Koordinatennetz, in das man das symme- 
trische Profil eingezeichnet denkt, den Zusammen- 
hang &e(©) und bestimmt nach (12b) die verlangten 
Fourierkoeffizienten der @Quellverteilung. Die 
Wirbelverteilung ist nach Gl. (12a) zu berechnen, 
indem man die Funktion 7(x) mit = — cos Odem 
Verlauf der Mittellinie entnimmt undderen Fourier- 
koeffizienten ermittelt. Wie beim Birnbaum- 
schen Verfahren kann man auch die differentiierte 
Gleichung benutzen, bei der die Ableitung n’(«) 
Verwendung findet. 

Zur Bestimmung der Strömungsgeschwindigkeit W 
am Profilrand hat man 


ov oy\: 
(+) 


auf die Koordinaten e und © umzurechnen, mit 


av. 
de 
Kraft und Moment ergeben sich aus den Blasius- 
schen Formeln. 
Hannover. 


e’ (9)-+ (9,9) —=(, 


K. Jaeckel. 


Die erste Kleinfunktionsrechenmaschine. 


Dank der Unterstützung durch die Deutsche 
Forschungsgemeinschaft konnte der Verfasser in der 
Werkstatt des Geodätischen Instituis der Tech- 
nischen Hochschule Stuttgart eine Funktions- 
rechenmaschine nach dem in Band 31 (1951), S. 301 
bis 309 dieser Zeitschrift beschriebenen Prinzip der 
einstufigen Interpolation herstellen. Als Grund- 
maschine wurde eine Sprossenradrechenmaschine 
Fabrikat Schubert mit Handantrieb verwendet. Ihr 
Einstellwerk wurde durch eine besondere Über- 
tragungseinrichtung mit einem Speicherwerk ver- 
bunden, dessen Wertträger aus 42 kreisförmigen 
„Stufenscheiben‘‘ besteht, auf deren Umfang die 
Ziffern der gespeicherten Grundwerte und Inter- 
polationsfaktoren durch verschieden tiefe Schlitze 
dargestellt sind. Der Wertträger (Durchmesser 
100 mm, Länge 108 mm) ist auswechselbar und ent- 
hält jeweils drei Funktionen mit je 100 siebenstelligen 
Grundwerten und Interpolationsfaktoren. Es wurden 
2 verschiedene Wertträger hergestellt, einer für die 
Funktionen sin &, cos © und arc tg x und ein zweiter 


für tg x, sec tg x (Secans aus Tangens) und / a2. 
Die Bedienung der Versuchsmaschine sei an folgen- 


dem Beispiel erläutert: Es soll sin 236°4357 be- 
rechnet werden. Man stellt den Funktionsschalter 
auf sin und Quadrant III und den Wertträger durch 
Drehen auf das Grundargument 236°. Dann wird 
durch einen Hebelzug der Grundwert (= |sin 2368 |) 
ausdem Speicherwerk in das Einstellwerk und von dort 
durch eine positive Kurbeldrehung in das Resultat- 
werk gebracht. Durch einen zweiten Hebelzug wird 
das Einstellwerk auf den Interpolationsfaktor ein- 
gestellt. Nun wird das Restargument (0,4357) in das 
Umdrehungszählwerk eingekurbelt, dann steht im 
Resultatwerk der Absolutbetrag des gesuchten sin- 
Wertes, dessen Vorzeichen durch die Stellung des 
Funktionsschalters angezeigt wird. Der berechnete 
Wert kann im Bedarfsfalle aus dem Resultatwerk 
mechanisch in das Einstellwerk übertragen werden. 

Soll aus einem sin-Wert der Winkel berechnet 
werden, so wird dieser Wert ins Resultatwerk ein- 
gestellt und um den dem Absolutbetrag nach nächst 


16* 
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kleineren Grundwert vermindert. Dann wird die 
Differenz durch den zum übertragenen Grundwert 
gehörigen Interpolationsfaktor dividiert, worauf an 
der Grundargumenteinstellung des Speicherwerks 
die vollen Grade und im Umdrehungszählwerk die 
Minuten und Sekunden des spitzen Winkels abgelesen 
werden können. Auf ähnliche Weise lassen sich die 
übrigen mechanischen Funktionstafeln für direkte 
und indirekte Interpolation verwenden. 

Das Quadratwurzelziehen erfolgt nach der Glei- 
chung z = + Az = 21 + (2? — x2)/2 x, wiefolgt: 
Man erzeugt im Resultatwerk x? und stellt das 
Speicherwerk auf den an einer Trommel ablesbaren 
nächst kleineren Wert z3 ein. Dann wird durch 
einen Hebelzug der zugehörige Wert x, in das Ein- 
stellwerk übertragen und durch negative Kurbel- 
drehungen mit sich selbst multipliziert. Dann steht 
im Umdrehungszählwerk 2, und im Resultatwerk 
x? — x@?. Hierauf wird, nachdem das Einstellwerk 
durch einen weiteren Hebelzug auf 2 x, eingestellt 
wurde, das Resultatwerk auf Null gekurbelt, dann 
kann im Umdrehungszählwerk x, + Az abgelesen 
werden. 

Ein Vergleich zwischen dem Rechnen mit der 
Funktionsrechenmaschine und dem Rechnen mit ge- 
wöhnlicher Rechenmaschine und gedruckter Funk- 
tionstafel (je 7 verschiedene Aufgaben mit je 10---20 
Zahlenbeispielen und 20---50 Versuchspersonen, ins- 
gesamt je etwa 3000 Rechnungen) zeigte, daß die 
neue Rechenart wesentlich einfacher, bequemer und 
sicherer ist als die alte. Die Zahl der durch Unauf- 
merksamkeit des Rechners verursachten Fehler war 
nur noch etwa halb so groß wie bei der bisher üblichen 
Berechnung. Sie gliedert sich wie folgt: 


Zahl der Fehler bei Argumenteinstellung. 35% 

Zabl der Fehler außerhalb Funktionsberech- 
ne ae 

Zahl der Fehler beim Bilden von Komplement- 
zählen #u..1,, So re a a 30 
Zahl der Bedienungsfehler 12% 
Zahl der Quadrantenfehler . sorr AR 9A 
Zahl der sonstigen Fehler. . . k 4% 
100% 


Das Rechnen mit der Funktionsrechenmaschine 
bringt auch bereits in der Versuchsausführung in den 
meisten Fällen eine beträchtliche Zeitersparnis 
(17---30%). Nur das Quadratwurzelziehen und die 
Berechnung von tg-Werten > 1 erfordern wegen der 
zeitraubenden Division vorläufig noch etwas mehr 
Zeit als die bisher übliche Berechnung. 

Die Genauigkeit der mechanischen Tafeln ist bei 
ben trigonometrischen Funktionen etwa fünfstellig. 
Bei der mechanischen Quadratwurzeltafel ist der 
Fehler < 2,6 - 10-2im Bereich 1 <x < 10. Weitere 
Einzelheiten können dem Bericht K. Ramsayer: 
Die geodätische Funktionsrechenmaschine, Zeit- 
schrift für Vermessungswesen 1952, Heft 4, ent- 
nommen werden. 


Stuttgart. K. Ramsayer, 


Über einige Hilfsfunktionen mit zwei Eingän- 
gen, die bei Lösung der ‚Gleichung der ein- 
dimensionalen gedämpften Welle‘ auftreten. 
Die Behandlung des linearisierten Systems in der 
Arbeit von E. Pieruschka ‚Die mathematischen 
Grundlagen zu einer Meßmethode des Schubmoduls 
zäher Flüssigkeiten“, Z. angew. Math. Mech. 31 (1951), 
S. 83—92, veranlaßt mich, über Untersuchungen zu 
berichten, die ich vor längerer Zeit ebenfalls über die 
„Gleichung der eindimensionalen gedämpften Welle“, 
die inder Elektrotechnik auch als ‚„Heavisidesche“ 
Gleichung bekannt ist, 

av _ 
028 

angestellt habe. 


oV 0:V 
dr a a EL) 
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Bei Vorgabe eines harmonischen Zeitablaufs führt 
die einfachere ‚Gleichung eindimensionaler Diffu- 
sions- und Wärmeleitvorgänge‘“ 


Er re) 
bekanntlich auf Lösungen 


Y(2,t) = e*!y0;V;i(e) ee h(2a): 


wobei die ‚Partikularlösungen‘“ V;(x) bestimmte 
Kombinationen von trigometrischen und Hyperbel- 
funktionen sind, die von nureine m Argument a ab- 
hängen, wobei a im wesentlichen der Veränderlichen x 
entspricht (in Verbindung mit bestimmten Ausgangs- 
daten der Gleichung). 

Die allgemeinere Gleichung (1) führt — ebenfalls 
bei Vorgabe eines harmonischen Zeitablaufes — auf 
Lösungen 


Va)=ertg Vila)... . (le), 
[7 


wobei diesmal die ‚„Partikularlösungen‘ V;(x) wie- 
derum bestimmte Kombinationen von trigometrischen 
und Hyperbelfunktionen sind, die aber jetzt vonz wei 
Argumenten, & und , abhängen. Dabei entspricht 
wiederum & im wesentlichen der Veränderlichen z (in 
Verbindung mit bestimmten Ausgangsdaten der Glei- 
chung), während p im wesentlichen der Größe b? von 
Gleichung (1) entspricht. — Im Falle 5? 0 strebt 9 
dem Grenzwert 1 zu. 

So lautet z. B. in einem bestimmten Beispiele !) die 
Endlösung im einfachen Falle (d.h. mit 5? = 0): 


3 Sina —sina 


ar © Coja — cosa a 
während sie im allgemeinen Falle !) lautet: 
3 „Ein &p—gpsin = 

EP re 2... (38). 


& 
Goiap — cos — 
jap 2 


Während bei g =1 [was dem Werte 5b? = 0, also 
der einfacheren Gleichung (2), entspricht] die 
Kurve monoton fällt, tritt mit kleiner werden- 


Fl(&,9) 


7 2 J 4 


dem Parameter (9<1) mindestens ein” Maximum 
auf, das aber zunächst (bis etwa 9 = 0,8) den Gesamt- 
verlauf der Kurve nicht wesentlich ändert, während 
bei noch kleiner werdenden Werten (z.B. bei 
p = 0,4) der Charakter der Kurve ein grundsätzlich 
ganz anderer wird. Besonders wichtig ist die Tat- 
sache, daß bei solchen kleinen Werten @ sehr steile 
Äste vorliegen, d.h. die Werte F gegen kleine Ände- 
rungen von & stark empfindlich werden. 


ı) Dieses Beispiel ist von Herrn Dipl.-Ing. Hennigerin 
seiner Arbeit: ‚„‚Wirbelstromverluste inferromagnetischen Stof- 
fen komplexerLeitfähigkeit‘‘, Deutsche Elektrotechnik 6 (1952) 
Heft 9 u. 10 behandelt worden. Herrn Dipl.-Ing. Henniger 
verdanke ich auch den Hinweis auf diese.'Problemstellung, wofür 
ihm auch an dieser Stelle nochmals gedankt sei. 
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In jedem Einzelfalle wird es also notwendig sein, 
den Verlaufsolcher vorkommender Funktionen F (x, p) 
für verschiedene Parameter p zu diskutieren. 

‚Eine wesentliche Hilfe bietet dabei die Kenntnis der 
Eigenschaften einiger Hilfsfunktionen, über die ich 
einige Ausführungen bringen will, 

Im einfacheren Falle kommen vorwiegend folgende 
Hilfsfunktionen vor: 


1 

5 [Sofa + cosa] = $(a) ; 

1 

Fo [Sin +sina]= T (a) ; 

1 

7 [Cojx —cosa]=U(o) ; 

l 

2 [Sina —sina] = V (a) 
Die hierfür angegebenen Bezeichnungen 8, T, U, V 
sind die im sowjetischen Schrifttum gebräuchlichen; 
sie wurden von A. N. Kryloff, Die Schwingungen von 
Schiffen, Gesammelte Werke, Bd. X, 1948, S. 231 in 
dieser Form eingeführt. 


Im allgemeineren Falle treten folgende analogen 
Hilfsfunktionen auf: 


:- Ir" Cof(ap) + kai cos )l = 8,(6, 9); 


N (6) 
u Ir" Sin(ap+ I sin (5) = T,(&, 9) | 
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7. Ir Eoi(ap) — () cos 6) =U,„(&, ) | 
1 n f 
= Ir" Sin (ag) — | sin 6) = 1.9 | 


Ähnlich wie für die Hilfsfunktionen (5) besonders 
einfache Differentialbezeichnungen gelten, nämlich: 
d Ä d m 8 
te) = rl); | 


(6) 


d d (Da) 
ZIM)=To); VW) =Ul) | 
gelten auch für die Hilfsfunktionen (6) einfache Diffe- 
rentialbezeichnungen, nämlich: 
[) 
0x Sn(&, p) > Vn+1(6, p) 3 
(7) 
35 Tn(&, 9) = 8n+1(0, 9) ; 


(] 
Em U„(&,9) = Tn+ı(0, ) ; 


(08). 


0) 
EM Vn(a, p) ZZ Un+1(a, 9) 


In der Deutschen Akademie der Wissenschaften, In- 
stitut für Mathematik, Abt. Angewandte Mathematik 
werden diese Hilfsfunktionen zur Zeit vertafelt. 


Berlin-Lichtenberg. K. Borkmann. 


BUCHBESPRECHUNGEN 


Dr. Milan Vidmar (Prof. a. d. Techn. Hochschule 
Ljublana), DieGestaltderelektrischen 
Freileitungen. (Lehr- und Handbücher der 
Ingenieurwissenschaften Bd. 21.) 199 S. mit 49 Abb. 
Basel 1952. Verlag Birkhäuser. Preis brosch. 16,65 Fr. 
gena1l9NTo Er 

Es ist immer begrüßenswert, wenn ein an sich be- 
kannter Gegenstand von einer anderen Seite her als 
bislang beleuchtet wird. Er lenkt damit von neuem 
auch das Augenmerk derer auf sich, die ihn längst 
vollkommen zu kennen glauben und regt zu erneuter 
Betrachtung an. Nicht selten zeigt er auch in dieser 
neuen Beleuchtung Seiten, die im bisherigen Licht 
nicht zu erkennen waren. In diesem Sinne verfährt 
der Verfasser, wenn er ‚‚in reifen Jahren, sozusagen 
von vorn anfangend‘“ sich mit der „Gestalt der 
elektrischen Freileitungen‘‘ beschäftigt. Er knüpft 
dabei an seine jahrzehntelange Beschäftigung mit 
dem Transformator, seinem Bau und seiner Be- 
ziehung zur Energieübertragungstechnik an und 
bringt tatsächlich überraschende Einblicke in das 
Wesen elektrischer Freileitungen, die den praktischen 
Freileitungsbauer interessieren und zur Überprüfung 
seiner eigenen, auf Erfahrung und Durchrechnung 
zahlreicher Einzelprojekte gegründeten Ergebnisse 
anregen. Solche Überprüfungen, wie sie z. B. in be- 
zug auf die wirtschaftliche Spannweite bereits statt- 
fanden, führen zwar keineswegs in allen Fällen zu 
einer restlosen Übereinstimmung mit der von Vidmar 
in den Mittelpunkt gestellten Formel wie in dem 
Kapitel 2.20 ‚„‚Bestätigung der Lösung des Spann- 
weitenproblems durch Erfahrungstatsachen‘‘. Immer- 
hin wird dem Leser vor Augen geführt, daß man das 
Problem des Freileitungsbaues mit seinen zahl- 
reichen Einflußfaktoren auch durch allgemeingültige 
Betrachtungen überblicken und daß man sich damit 
viel Mühe sparen kann, die man bisher mit einer Un- 


zahl von Durchrechnungen von Einzelfällen hatte. 


Das Eingeständnis des Verfassers am Schluß des 
Buches, daß er sich mit Problemen beschäftigt, ‚‚die 
sehr lebendig sind‘ und ‚die ganz gewiß keinen 


Ewigkeitswert beanspruchen dürfen“, ist sicher zu- 
treffend. Das Buch wird, wenn es in die Hände der 
Leitungsbauer kommt, aufrüttelnd wirken und durch 
die erneut eingeleitete Diskussion zu besseren Lö- 
sungen führen. Leider enthält es keinerlei Literatur- 
hinweise, so daß es für den Leser schwer zu erkennen 
ist, auf welchem Stand der Technik die Erörterungen 
aufgebaut sind. 


Dresden. F.Obenaus. 


Dr.Alfred Vogel, Klassische Grundlagen 
der Analysis. X-- 194 Seiten mit 17 Abb. 
Leipzig 1952. S. Hirzel-Verlag. Preis kart. 8,50 DM 


Das Buch stellt den ersten Teil einer zweisemestri- 
gen Ergänzungsvorlesung über Analysis dar, wie 
sie der Verfasser an der Technischen Hochschule 
Stuttgart vor Studenten der Mathematik und Physik, 
die daneben die entsprechende !Kursvorlesung 
„Höhere Mathematik“ hören, mehrfach gehalten hat. 

Es handelt sich im wesentlichen um den ‚‚gene- 
tischen Aufbau des Systems der reellen Zahlen aus 
dem axiomatisch begründeten System der natür- 
lichen Zahlen‘ mit dem Ziel, die archimedische 
Ordnung, die Abgeschlossenheit und die Vollständig- 
keit der reellen Zahlen zu erhalten. Eine kurze 
Einführung in die Mengenlehre führt dann noch bis 
zur transfiniten Kardinalzahl. 

Da es meistens in den Vorlesungen über Diffe- 
rential- und Integralrechnung aus Zeitmangel nicht 
möglich ist, diese Dinge ausführlich zu behandeln, 
ist es um so mehr zu begrüßen, nunmehr dem an- 
gehenden Studenten das vorliegende Buch emp- 
fehlen zu können, das, ohne Vorkenntnisse zu er- 
fordern, die Grundlagen der Analysis mit aller 
Strenge, dabei didaktisch sehr geschickt und le- 
bendig darbietet. Es eignet sich darüber hinaus 
auch für jeden gebildeten Laien, der eine gewisse 
Übung im abstrakten Denken besitzt. 

Die einzelnen Kapitel behandeln: 1. Die na- 
türlichen Zahlen, 2. die Brüche und die absolut- 
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rationalen Zahlen, 3. die Unterschiede und die ra- 
tionalen Zahlen, 4. die Schachtelungen und die 
reellen Zahlen, 5. die Mengen und die Kardinal- 
zahlen. — Diesem ersten Teil soll ein zweiter über 
die ‚‚Grundbegriffe der Analysis‘ (Funktionsbegriff 
bis zum Integralbegriff) folgen. 


Dresden. Maruhn. 


R.Rothe, HöhereMathematikfürMa- 

thematiker, Physiker und Inge- 
nieure (Teubners mathematische Leitfäden Bd. 21 
bis23).Teill:Differentialrechnung und 
GrundformelnderIntegralrechnung 
nebstAnwendungen. 10. Auflage, herausgeg. 
von Dr. W. Schmeidler (Prof. a.d. T. U. Berlin). 2118. 
mit 161 Abb. Leipzig 1952. B. G. Teubner Verlags- 
gesellschaft. Preis geb. 6,30 DM. 
Teilll:Integralrechnung. Unoendliche 
Reihen. Vektorrechnung nebst An- 
wendungen. 38, Auflage herausgeg. von Dr. 
W. Schmeidler (Prof. a. d. T. U. Berlin), 210 S. mit 
98 Abb. Leipzig 1952, B. G. Teubner Verlagsgesell- 
schaft. Preis geb. 5,60 DM. 
Teillll. Flächen im Raume. Linien- 
integrale und mehrfache Integrale. 
Gewöhnliche und partielle Diffe- 
rentialgleichungen nebst Anwend- 
dungen. 5. Auflage, herausgeg. von Dr.W.Schmeid- 
ler (Prof. a. d, T. U, Berlin), 236 S. mit 167Abb. Leip- 
zig 1952. B.G.Teubner Verlagsgesellschaft. Preis geb. 
5,80 DM. 

In dieser Neuauflage des bekannten, von Prof. Dr- 
W. Schmeidler herausgegebenen Werkes wurden 
gegenüber den früheren Auflagen neben einer Reihe 
kleinerer Verbesserungen in TeilI und Teil II die 
Abschnitte ‚‚Integrationsgrenzen bei neuen Ver- 
änderlichen‘ und ‚Eindeutigkeitssatz für Fourier- 
sche Reihen‘ neu aufgenommen. 


Dresden. Maruhn. 


K. Oswatitsch, (Dozent a. d. TH Stockholm), Gas - 
dynamik. XIII +4568. mit 300 Abb. und 
3 Tafeln. Wien 1952. Springer-Verlag. Preis geb. 
78,— DM. 

Das Hauptgewicht in diesem vorzüglichen Werk, 
in dem das in den letzten 50 Jahren angesammelte 
Wissen über die Aerodynamik der kompressiblen 
Medien verarbeitet ist, liegt nicht auf der Darstellung 
der mathematischen Methoden und der numerischen 
Berechnungsverfahren, sondern auf der anschau- 
lichen Behandlung der physikalischen und technischen 
Probleme. 

Das Buch ist in 12 Kapitel gegliedert: Grund- 
begriffe der Thermodynamik; stationäre Stromfaden- 
theorie (‚„hydraulische‘‘ Behandlung der stationären 
Gasströmung); nichtstationäre Stromfadentheorie 
(einschließend die ebenen, zylindrischen und sphä- 
rischen akustischen Wellen); Gleichungen und Sätze 
der allgemeinen Strömung (einschließlich Verdichtungs- 
stoß); spezielle Anwendungen der Integralsätze (z. B. 
Schaufelgitter, Düsenschub, - Raketenantrieb); spe- 
zielle strenge Lösungen der reibungsfreien stationären 
Strömung (Prandtl-Meyersche Eckenströmung, Trans- 
formationen von Legendre, Molenbroek 
und Tschapligin); allgemeine Behandlung der 
reibungsfreien, stationären Unterschallströmung im 
ebenen und im drehsymmetrischen Fall; dasselbe für 
die Überschallströmung (einschließlich Verdichtungs- 
stoß, kegelige Felder, Strömungen mit veränderlicher 
Entropie, Charakteristikenmethoden, Hyperschall- 
strömung); schallnahe Strömungen (Ähnlichkeits- 
gesetze, stoßfreie lokale Überschallgebiete); spezielle 
stationäre und nichtetationäre räumliche Strömungen 
(z. B. tragendes Dreieck in Überschallstrom); Strö- 
mungen mit Reibung (Grenzschichttheorie); Über- 
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blick über die Versuchstechnik (einschließlich elek- 
trische und andere Analogien). Den einzelnen Ka- 
piteln sind ausführliche Literaturverzeichnisse an- 
gefügt. 

Nicht nur wegen seines umfassenden Inhalts, son- 
dern ebensosehr wegen der ausgezeichneten Dar- 
stellung wird dieses Werk sowohl für den Theoretiker 
wie für den in der Praxis stehenden Ingenieur von 
größtem Nutzen sein. Es enthält viele eigene For- 
schungsergebnisse des Verfassers aus seiner früheren 
Tätigkeit bei L. Prandt| und aus der nachfolgen- 
den Zeit. 


München. Sauer. 


Dr. Ing. Bruno Eck, Ventilatoren. Ent- 
wurfund Betrieb der Schleuder- und 
Schraubengebläse. 3048. mit 344 Abb. 
2. verb. und erweiterte Auflage. Berlin-Göttin- 
gen-Heidelberg 1952. Springer-Verlag. Preis geb. 
36,— DM. 

Seit der 1937 erschienenen ersten Auflage des wohl- 
bekannten Werkes hat der Ventilatorenbau eine 
außerordentlich starke Entwicklung erlebt. Die an- 
fängliche Überlegenheit des Axialgebläses ist von dem 
Radiallüfter fast vollkommen aufgeholt worden. 
Diese Fortschritte in der wissenschaftlichen Erkennt- 
nis und konstruktiven Gestaltung haben ihren Nieder- 
schlag in der vorliegenden Neuauflage gefunden. Der 
technischen Entwicklungstendenz entsprechend hat 
der Verfasser besondere Sorgfalt den Abschnitten 
über die Theorie der Strömung in Laufkanal und 
Diffusor und die beste Gestaltung dieser Bauelemente 
zur Erzielung hohen Wirkungsgrades unter Heran- 
ziehung der neuesten Forschungsergebnisse gewidmet. 
Beachtung verdient die Herausstellung der besonderen 
konstruktiven Erfordernisse der Lüfter hinsichtlich 
Geräuscharmut, Platzbedarf, Regulierfähigkeit und 
Parallelbetrieb. Die vom Verfasser vorgeschlagenen 
Kennziffern lassen die geeignete Grundlage zur Be- 
urteilung dieser Verhältnisse erhoffen. Zu begrüßen 
ist die Ergänzung der früheren Auflage durch die 
Aufnahme der neuen Entwicklungen der Radial- und 
Axialgebläse in Form der meridian beschleunigten 
Axialgebläse, der Radialgebläse mit Vorläufern, der 
Querstromgebläse, Kleinstgebläse u. a. Das Werk be- 
schäftigt sich zwar ausschließlich mit den einstufigen 
Gebläsen, doch bilden die Ausführungen ohne weiteres 
auch die Grundlagen für die Gestaltung mehrstufiger 
Kreiselverdichter. Angesichts des Mangels an hoch- 
wertiger Literatur auf dem Gebiete der Kreiselver- 
dichter muß das Erscheinen dieses Buches sowohl von 
dem Konstrukteur als auch dem Betriebsingenieur 
dankbar begrüßt werden. 


Dresden. HoEaltin 


Helmut Hasse, Mathematik als Wissen- 
schaft, Kunst und Macht. 368. Wies- 
baden 1952. Verlag für angewandte Wissenschaften. 
Brosch. DM 4,80. 


Die Schrift will nicht wissenschaftliche Erkennt- 
nisse, sondern persönliche Bekenntnisse darstellen. 
Dieser betonte Einsatz der Mathematiker-Persönlich- 
keit und ihrer seelischen Schaffensimpulse im In- 
teresse der rechten Bewertung ihrer Wissenschaft ist 
lebhaft zu begrüßen. Auch ist wärmester Beherzigung 
zu empfehlen, was der Vf. über die an die Darstellung 
der Mathematik zu stellenden Anforderungen aus- 
führt. Forderungen, die gemeinhin als nur didaktisch 
angesehen und wegen dieses Signums vom „Wissen- 
schaftler‘“ oft abgewiesen werden, erhebt Hasse in 
die Sphäre des ‚„Lebensgesetzes der Mathematik“ 
selbst. Er spricht von der ‚„Tragik des Mathematikers 
durch Streben nach Formalisierung‘ in dem Sinne, 
daß durch die Überbetonung der äußeren Form „der 
eigentliche mathematische Inhalt übertönt wird‘, 
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Gerade aus der Feder des modernen Algebraikers und 
Zahlentheoretikers wiegt ein solcher Hinweis besonders 
schwer. 

Angesichts dieser Vorzüge brauchen gewisse Mängel 
und Unklarheiten nicht verschwiegen zu werden. Auf 
zweierlei sei kurz hingewiesen: 

Wenn die Mathematik, weil sie nicht in der „Fähig- 
keit des Rechnens und Vermessens‘“ besteht, als 
„‚Geisteswissenschaft“ bezeichnet wird, so stellt diese 
These genau denselben fruchtlosen Streit um Worte 
dar wie 1925 bei ihrer erstmaligen Formulierung durch 
den Vf. Das der Schrift gestellte Ziel, über ‚den 
tieferen Sinn der Mathematik und ihre Stellung im 
Ganzen der Kultur“ Verständnis zu wecken, wird 
dadurch nicht erreicht — ganz abgesehen davon, 
inwieweit das durch persönliche Bekenntnisse über- 
haupt möglich ist. 

Zum anderen möchte ich ergänzen, was über das 
Ansehen der Mathematik in Deutschland gesagt wird. 
Es erreichte seinen Tiefstand zur Zeit des Nazismus, 
und die spätere Anerkennung der Mathematik als 
„kriegswichtig‘ kann im Blick auf das Ganze der Kul- 
tur nicht gerade als grundsätzliche Meinungsänderung 
bezeichnet werden. Falls solche Erinnerungen im Vf. 
nachklingen sollten — er selbst erwähnt nur die 
Romantik als Ansehen vermindernd —, muß doch er- 
gänzend hinzugefügt werden, daß darin heute ein ent- 
scheidender erfreulicher Wandel eingetreten ist. 


Dresden. M. Draeger. 


Dr. R.M. El Badrawy, Ebene Plangitter 
bei Überschallgeschwindigkeit. (Mit- 
teilungen aus dem Institut für Aerodynamik an der 
Eidgenössischen Technichen Hochschule in Zürich, 
Nr. 19.) 90 S. mit 71 Abb. und 11 Tab. Zürich 1952. 
Verlag Leemann. 

Die unter Leitung von Prof. Ackeret durch- 
geführte wertvolle Arbeit behandelt als Grundlage 
der Gitterströmung bei Überschall das ebene Plan- 
gitter. Zunächst wird für Ma = 1,2 bis 10 und ver- 
schiedene Anstellwinkel Auftrieb und Widerstand 
einer unendlichen dünnen ebenen Platte exakt be- 
rechnet, in ausführlichen Tabellen und Kurven- 
darstellungen verzeichnet, mit der linearisierten 
Theorie nach Ackeret verglichen und dabei gute 
Übereinstimmung festgestellt. Es folgt eine analy- 
tische Behandlung des für Überschall-Gitterströmung 
wichtigen Einholens, Durchkreuzens und Reflek- 
tierens von Verdichtungsstößen und Expansions- 
wellen. Auf diesen Grundlagen beruht die anschließen- 
de Behandlung des Gitterproblems, wieder in exakter 
und linearisierter Form. Ein Zahlenbeispiel liefert 
wieder gute Übereinstimmung. Zur Illustration 
dienen Schlierenaufnahmen zwischen zwei parallelen 
Platten für mehrere Anstellwinkel. Als ein instatio- 
närer Fall wird die plötzliche Veränderung des Anstell- 
winkels behandelt. Zum Schluß wird noch Reibungs- 
und Dickeneinfluß auf den wie bei einem Propeller de- 
finierten Wirkungsgrad untersucht und der Wirkungs- 
grad eines Überschallpropellers abgeschätzt. 


Weilburg. L. Schiller. 


HansReichard, Gesetzmäßigkeiten der 
freien Turbulenz. (VDI-Forschungsheft 414) 
2. erweiterte Aufl. II +30S. mit 18 Abb. Düssel- 
dorf 1951. Deutscher Ingenieur-Verlag G. m. b. H. 
Preis brosch. 12,50 DM. 

Die früheren rechnerischen Arbeiten über freie 
turbulente Strömungen fußten auf versuchsweisen 
Ansätzen über die turbulente Schubspannung 
(Mischungsweg) und leiteten „‚deduktiv‘ hieraus ge- 
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wisse Eigenschaften der Strömung ab. Hierbei auf- 
tretenden Schwierigkeiten begegnet Reichardt 
durch eine mehr ‚‚induktive“ Betrachtungsweise, 
indem er umgekehrt die Theorie aus den Versuchs- 
ergebnissen entwickelt, wodurch er in die Lage kommt, 
nur schwer oder gar nicht meßbare Größen zu be- 
rechnen: Quergeschwindigkeit, die turbulente Impuls- 
übertragung (Scheinreibung) und die gesamte Impuls- 
übertragung. Erste Ergebnisse dieser Behandlungs- 
weise sind in dieser Zeitschrift Bd. 21 (1941), S. 257 
veröffentlicht, weitere Ergebnisse in einem nur wenig 
zugänglichen Heft in Bd. 24 (1944), S. 268. Eine 
breitere zusammenfassende Darstellung bietet das 
jetzt in 2. Auflage vorliegende VDI-Forschungsheft. 
Über den Inhalt der 1. Auflage hinaus enthält es einen 
Abschnitt über die Geschwindigkeitsverteilung in der 
x-Richtung, einen über Impuls- und Wärmeaustausch 
(nach Z. angew. Math. Mech. 24) und einen völlig 
neuen über das Problem eines allgemeinen Schub- 
spannungsansatzes. In diesem setzt sich der Verfasser 
grundsätzlich und kritisch mit den verschiedenen hier- 
zu gemachten Versuchen auseinander, ohne selbst 
einen befriedigenden Ausdruck entwickeln zu können. 
Im Sinn seiner „induktiven‘ Betrachtungsweise for- 
dert er für den theoretischen Fortschritt als not- 
wendig die Weiterbildung der Versuchstechnik, um 
eine genauere Kenntnis der Profile und ihrer Ab- 
leitungen zu erhalten. 

Allen mit den Turbulenzfragen Beschäftigten wird 
die neue erweiterte Auflage des Forschungsheftes sehr 
willkommen sein. 


Weilburg. IesSchrller, 


Prof. Dr. Fritz Neiß, Einführung in die 
Zahlentheorie. VIII+113S. Leipzig 1952. 
S. Hirzel-Verlag. Preis kart. 4,80 DM. 


Was jeder Mathematiker von der Zahlentheorie 
wissen muß, soll nach der Absicht des Vf. in seinem 
Buche dargestellt werden. Es dürfte kaum zwei 
Autoren geben, die bei dieser Zielsetzung denselben 
Inhalt vorlegen würden. Deshalb wäre es müßig und 
objektiv kaum möglich, die betroffene Stoffauswahl 
zu kritisieren. 

Teilbarkeitseigenschaften, Kongruenzen, quadra- 
tische Reste, quadratische Formen und einige zahlen- 
theoretische Funktionen bilden den wesentlichen In- 
halt. Viel könnte daran auch kaum geändert wer- 
den. Daß einige grundlegende Begriffe der modernen 
Algebra eingearbeitet sind, ist nur zu begrüßen. 

Ein anderes ist die Art der Darstellung. Ihr fehlt — 
und zwar gerade für den Mathematiker mit „ein 
wenig Zahlentheorie zur allgemeinen Bildung‘ — das 
Wesentliche: der systematische Zusammenhang, der 
zur Wissenschaft konstitutiv gehört. In der vorliegen- 
den Form gewinnt der Leser, an den der Vf. sich vor- 
nehmlich wendet, den Eindruck, daß Zahlentheorie 
aus einer Reihe zwar interessanter, aber einzelner 
Probleme besteht. Und dieser — objektiv falsche — 
Eindruck wird verstärkt, wenn ein Kapitel die aus- 
drückliche Überschrift „Einzelne Sätze‘ trägt. Unter 
diesen erscheint als erster der kleine Fermatsche Satz, 
der gerade alles andere als ein „einzelner Satz‘ ist. 
Ein in sich geschlossenes, plastisches und darum un- 
verlierbares Bild von der Zahlentheorie und ihrem 
Wesen ist leider nicht gestaltet. 

Schließlich noch ein Wunsch: Das Buch enthält 
viele, z. T. sinnstörende Druckfehler. Sie in einem 
Nachtrage zu berichtigen, würde das Verständnis 
nicht unwesentlich erleichtern. 


Dresden. Draeger. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen, 
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EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.R.Reissig, Die pandiagonalen Qua- 
drate vierter Ordnung. (Berichte über 
die Verhandlungen der Sächsischen Akademie der 
Wissenschaften zu Leipzig. Mathematisch-natur- 
wissenschaftliche Klasse. Band 100, Heft 6) 54 8. 
Berlin 1952. Akademie-Verlag. Preis brosch. 5,— DM. 


Dr.-Ing.W.Fries,Fachwerk undRahmen- 
werk. Ein systematischer Grundriß der Statik des 
ebenen Tragwerkes. X + 3688. mit 365 Abb. 
Berlin-Göttingen-Heidelberg 1953. Springer-Verlag. 
Preis geb. 42,— DM. 


Dr. K. Mayrhofer (Prof. a. d. Universität Wien). 
Inhalt und Maß. VIII + 269 S. mit 17 Abb. 
Wien 1952.Springer- Verlag. Preis brosch. 36,— DM., 
geb. 39,— DM. 


Prof. Dr. Ulrich Graf und Dr. H.-J. Henning, 
Formeln und Tabellen der mathe- 
matischen Statistik. VI + 102S. mit 
9Abb. Berlin-Göttingen-Heidelberg 1953. Springer- 
Verlag. Preis brosch. 9,— DM. 


Wissenschaftliche Zeitschrift der Universität Leip- 
zig. Mathematisch-naturwissenschaftliche Reihen 
Nr. 1/2. 168 S. Selbstverlag der Universität. 


Ing. Simon Gershanik, Contribucion al 
estudio de cargas sismicas sobre 
edifictios. (Publicaciones del Observatoro 
Astronomico de la Universidad Nacional de la Plata. 
Serie Geofisica VIII) 49 S. mit 3 Abb. La Plata 1951. 


Deuxi&me Colloque de G6omö6trie Algöbrique. (Tenu 
a Liege les 9, 10, 11 et 12 juin 1952.) Centre Belge de 
Recherches Mathematiques. 244 S. Liege-Paris 1952. 
Georges Thone — Masson & Cie. Preis brosch. 
375 belg. Francs oder 2.625 franz. Francs. 


Dr. D. Küchemann — Dr. Johanna Weber, Aero- 
dynamics of Propulsion. IX+3408. 
mit Abb. New York, London, Toronto. 1953. 
McGraw - Hill Publishing Company Ltd. Preis 
64/6 Schilling. 


Dr.-Ing. W. Matz, Berechnung der Aus- 
mauerung stählerner Gefäße. (Ver- 
fahrenstechnik in Einzeldarstellungen. Heft 1.) 
VIII +75S. mit 15 Abb. Berlin-Göttingen-Heidel- 
berg 1953. Springer-Verlag. Preis brosch. 10.50 DM. 


M. Landolt, Größe, Maßzahl und Ein- 
heit. Band 1. Zweite Auflage. 127 Seiten. Zürich 
1952. Rascher-Verlag. DM 8.—. 


M. Bergsträsser, Evolventengeometrie 
für Stirnradgetriebe. (VDI-Forschungs- 
heft 436.) 24 +38S. mit 41 Abb. und 3 Zahlen- 
tafeln. Düsseldorf 1952. Preis 15.— DM, für VDI- 
Mitglieder 13.50 DM. 


NACHRICHTEN 


Arpad Nadai zum 70. Geburtstage 


Am 3. April dieses Jahres begeht ArpadNadai 
in Pittsburgh seinen 70. Geburtstag, In mehr als 
50 Veröffentlichungen über die Mechanik der elasti- 
schen und der plastischen Verformung und über die 
Mechanik ihres Zwischenbereiches hat er in seiner 
stets eigenständigen Weise die einen Probleme gelöst 
und andere gefördert und geklärt. 

Geboren in Budapest, hat er in Zürich, haupt- 
sächlich bei Stodola, Maschinenbau studiert und 
später in Berlin bei Eugen Meyer mit einer Disser- 
tation über die Temperatureffekte verformten Eisens 
promoviert. In Berlin hat er auch begonnen, sich mit 
den elastischen Platten und Schalen zu beschäftigen 
und die „Ersatzscherkräfte‘‘ als gedankliches Hilfs- 
mittel auszunutzen. 

Nach militärtechnischem Kriegsdienst in der 
österreich-ungarischen Armee ist er 1919 in das 
von Felix Klein gegründete Institut für ange- 
wandte Mechanik in Göttingen gekommen, wo er 
sich 1923 habilitierte und 1926 zum a. o. Professor 
ernannt wurde, In den Jahren bis 1927 hat er hier 
in enger Fühlung mit Ludwig Prandtl, unterstützt 
durch die einschlägigen wissenschaftlichen und tech- 
nischen Gesellschaften, mit seinen Doktoranden ge- 
arbeitet u. a. über dicke Platten unter konzentrierter 
Belastung und über die qualitativen und quantita- 
tiven Effekte (Fließfiguren und Spannungsgesetze) 
an der Fließgrenze zäher Metalle. Aus diesen Jahre 
stammt auch sein Buch, ‚‚Der bildsame Zustand der 
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Werkstoffe‘, in dem wohl erstmalig die Plastizität 
der Metalle in der eine mathematische Theorie ge- 
stattenden Klarheit behandelt ist. 

Nach einer Gasttätigkeit in der Westinghouse 
Comp. siedelte Nadai 1929 endgültig nach Pitts- 
burgh über. Hier beschäftigte er sich eingehend mit 
dem Kriechen der Metalle, mit der Gültigkeit der 
modernen Ansätze hierzu (Sinusgesetz von Prandtl; 
Oktaeder-Schubspannungen gleich Funktion der 
Oktaeder-Schiebungsgeschwindigkeit) und mit deren 
Anwendung in der Theorie des Walzens. Diese Ar- 
beiten sind in Deutschland nicht so bekannt und aus- 
genutzt worden, wie sie es verdienen. 

Möge Herr Na dai noch lange am Gebäude seiner 
Wissenschaft weiterbauen können und der Früchte 
seines Schaffens sich erfreuen. 


Leverkusen. W. Lode. 


Dresden: Zu seinem 70. Geburtstag wurde 
Professor Dr. Willers der Doktor der Natur- 
wissenschaften Ehrenhalber von der Technischen 
Hochschule Darmstadt verliehen. 


Hannover: Am9. Januar 1953 starb Professor 
Dr. Conrad Müller. 


Frankfurt/M.: Herrn Prof. Dr. Lore y wurde 
aus Anlaß seines 80. Geburtstages am 23. Jan. 1953 
das Verdienstkreuz verliehen. 
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